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ABSTRACT 

We argue that discrete dynamics has natural hnks to the theory of analytic functions. Most important, 
bifurcations and chaotic dynamical properties are related to intersections of algebraic varieties. This paves the 
way to identification of boundaries of Mandelbrot sets with discriminant varieties in moduli spaces, which are 
the central objects in the worlds of chaos and order (integrability) respectively. To understand and exploit this 
relation one needs first to develop the theory of discrete dynamics as a solid branch of algebraic geometry, which 
so far did not pay enough attention to iterated maps. The basic object to study in this context is Julia sheaf 
over the universal Mandelbrot set. The base has a charateristic combinatorial structure, which can be revealed 
by resultant analysis and represented by a basic graph. Sections (Julia sets) are contractions of a unit disc, 
related to the action of Abelian Z group on the unit circle. Their singularities (bifurcations) are located at the 
points of the universal discriminant variety. 
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1 Introduction 



This paper is devoted to the structure of Mandelbrot set a remarkable and important personage of the modern 
theoretical physics, related to chaos and fractals and simultaneously to analytical functions, Riemann surfaces 
and phase transitions. This makes Mandelbrot set one of the bridges connecting the worlds of chaos and order 
(integrability) . At the same time the Mandelbrot set is a very simple object, allowing detailed experimental 
investigation with the help of the easily available computer programs and profound theoretical analysis with 
the help of undergraduate-level mathematics. This makes Mandelbrot set a wonderful subject of physical 
investigation: one needs to look for theoretical explanations of experimentally observable properties of the 
object and theoretical predictions can be immediately verified by the easily affordable new experiments. In this 
paper we just lift the cover of Pandora box: if interested, the reader can easily continue along any of the research 
lines, which are mentioned below, or, digging a little deeper into the experimental data, find and resolve a lot 
of new puzzles, - all this without any special mathematical background. Remarkably, when thinking about this 
simple subject one appears in the close vicinity of modern problems of the "high theory": the subject is not 
only simple, it is also deep. 

Mandelbrot set shows up in the study of trajectories (world lines) behavior change under a change of the 
laws of motion. Generally speaking, trajectory depends on two types of data: on the motion law (e.g. on the 
choice of the Hamiltonian) and on the initial (boundary) conditions. The pattern of motion (the phase portrait 
of the system) can be rather sophisticated: trajectories can tend to fixed points, to limit cycles, to strange 
attractors and go away from the similar types of structures, some trajectories can go to infinity; trajectories 
can wrap around each other, knot and unknot etc. In this paper the information about the phase portrait will 
be collected in Julia sets, the structure of these sets encodes information about the world lines dependence on 
the initial conditions. If the law of motion is changed, the phase portrait = Julia set is also deformed. These 
changes can be pure quantitative (due to the variations of the shape of trajectories, of foci and limit cycles 
positions) and - sometime - qualitative (looking as reshuffling of the phase portrait, creation or elimination 
of the fixed points, cycles, attractors etc). The Mandelbrot set or its boundary, to be precise, is exactly 
the set of all the points in the laws-of-motion space (say in the space of Hamiltonians) where the qualitative 
reshufhings (phase transitions) take place. Therefore the problems about the Mandelbrot set (or, better, about 
the set of the phase portraits, associated with various laws of motion, i.e. about the "sheaf of Julia sets 
over the Mandelbrot set") are the typical problems about the structure of the space of theories (laws 
of motion), about what changes and what remains intact in transition from one theory to another. In other 
words, these are typical problems from the string theory^ - the science which studies the various sets of similar 
physical models and studies the change of the experimentally measurable characteristics with the change of the 
model. Moreover, in the case of the Mandelbrot set one can immediately address the problems from a difficult 
branch of string theory, which nowadays attracts a good deal of attention, - from landscape theory, which 
wonders how often is a given property encountered in the given set of models (i.e. associates with the various 
physical quantities, say, scattering amplitudes or mass spectra, the measures on the space of theories, which 
specify how frequent are the amplitudes of the given type or masses of the given values in the models of the 
given type). One of the lessons, taught by the study of Mandelbrot sets, is that even if the phase transitions 
are relatively rare in the space of theories, the phase transition points are not distributed in a uniform way, 
they tend to condense and form clusters of higher dimensions, up to codimension one, which actually separate 
- the naively un-separated - domains in the space of theories and make this space disconnected. This lesson 
is important for the theory of effective actions, for problems of their analytical continuations, for construction 
of special functions (r-functions) which can describe these effective actions, and for many other - conceptual 
and practical - applications of string theory. Dynamics, essential for this type of applications of Mandelbrot-set 
theory, is not obligatory the ordinary dynamics in physical space-time, often the relevant one is dynamics in 
the space of coupling constants, i.e. the (generalized) renormalization group flow. 

In this paper we restrict consideration to discrete dynamics of a single variable. This restriction preserves 
most essential properties of the subject, but drastically simplifies computer simulations and mathematical 
formalism (for example, substitutes the study of arbitrary analytic functions by that of polynomials or the 
study of arbitrary spectral surfaces of arbitrary dimension by that of the ordinary Riemann surfaces). However, 

^As often happens, the name is pure historical and refers to a concrete set of models, which were first successfully approached 
from such direction. Better name could be "the theory of theories" (sometime "abbreviated" to pretentious "theory of everything") 
or, most adequate, "generic quantum field theory". Unfortunately, these names are already associated with narrower fields of 
research, where accents are put on somewhat different issues. 
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throughout the text we use every chance to show directions of possible generahzations. Getting rid of the infinite- 
dimensional algebras (loop algebras), which would be unavoidable in considerations of continuous dynamics, 
we concentrate instead on another infinity: unification of various one-dimensional Mandelbrot sets, shown in 
numerous pictures below and associated with particular 1-dimensional families of evolution operators, into a total 
(infinite-dimensional) universal Mandelbrot set, shown symbolically in Fig^ni^^t the very end of the paper. 
Understanding the structure of this set (which has a pure cathegorial nature of the universal discriminant 
variety) can be one of the intermediate targets of the theory. The notions of discriminant and resultant are 
assumed familiar from the courses on general algebra and are actively used already in the Introduction. Still, 
if necessary, all the necessary definitions and properties can be found in s.4.12 of the present paper (the section 
itself is devoted to generalizations from polynomials to arbitrary analytic functions). 
The theory of deterministic dynamical systems J^, 

i'=(3'{x), (1) 

is one of the main chapters of theoretical and mathematical physics. One of the goals of the theory is to classify 
the types of motion - from completely integrable to fully ergodic, - and understand their dependence on the 
choice of the vector field /3(x) and initial conditions. Discrete dynamics, in the avatar of Poincare map, 

— > fix), 

P'{f{x))^^P\x), (2) 

naturally arises in attempts to develop such a theory and plays the principal role in computer experiments 
(which - in the absence of adequate theory - remain the main source of information about dynamical systems). 

In quantum field/string theory (QFT) the most important (but by no means unique) appearance of the 
problem |^ is in the role of renormalization group (RG) flow. As usual, QFT relates to a problem, 
formulated in terms of effective actions, this time to a Callan-Symanzik equation for a function !F{x^ (p) of the 
coupling constants x^ and additional variable if (having the meaning of logarithm of background field) 

In other contexts solutions to this equation are known as characteristics ,2^ of the problem Q : in the case of one 
variable the solution is arbitrary function of ip — s{x) or of x(s = if), its ip dependence is induced by the action of 
operator exp {—f^iP^di). Until recently only simple types of motions (attraction and repulsion points) were 
considered in the context of RG theory, but today this restriction is no longer fashionable (see, for example, 
PI for the first examples of RG with periodic behaviour and for further generalizations). In this context 

the discrete dynamics is associated with the original Kadanoff-style formulation of renormalization group, and 
Callan-Symanzik equation (PJ) becomes a finite-difference equation.^ Bifurcations of dynamical behaviour of 
trajectories with the change of /3(a;) in ^ are associated with transitions between different branches of effective 
actions 0. One of the most intriguing points here is that the theory of effective actions is actually the one of 

^This discrete equation looks like 

Its generic solution is arbitrary function of any of its particular solutions. Formally, such particular solutions are given by the 
differences (/p =p n{x), where h{x) is the "time" , needed to reach the point x, i.e. a solution of the equation /°" (xq) = x for given / 
and XQ. In other words, particular solutions are formally provided by the ip-th image and Lp-fa preimage of /(x), 

.F±(x,¥>) = r(Tv'(xo 

For example, if /(x) = x'^ and /°(x) = x'^ , then the generic solution to is an arbitrary function of the variable cT*" logx. 

Consideration of the present paper provides peculiar 93-independent solutions (zero- modes), associated with periodic /-orbits. 
For example, the "orbit (5-function" 

T±{x,^) = F^{x)S{F^{x)) 
is a solution of such type (see eq. 1231 below) - a generalization of continuous-case zero-mode 



n^,^) = Aet(^^ <5(/3(x)). 



At least in this sense periodic /-orbits are discrete analogues of non-trivial zeroes of /3-function, in continuum limit they can turn 
into limit cycles and strange attractors. 
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integrahle systems - effective actions usually are r-functions, RG flows are multi-directional (since the boundary 
of integration domain in functional integral can be varied in many different ways) and related to Whitham 
dynamics jH|- Obvious existence of non-trivial RG flows (describing, for example, how a fractal variety changes 
with the change of scale) ~ along with many other pieces of evidence - implies existence of a deep relation between 
the worlds of order and chaos (represented by r-functions and dynamical systems respectively and related by 
the - yet underdeveloped - theory of effective actions). Description of effective actions (prepotentials) JF(<y9,x) 
in situations when |^ exhibits chaotic motion remains an important open problem. Of course, one does not 
expect to express such T through ordinary elementary or special functions, but we are going to claim in this 
paper that the adequate language can still be searched for in the framework of algebraic geometry, in which the 
theory of dynamical systems can be naturally embedded, at least in the discrete case. 

The goal of this paper is to overview the connection between the discrete dynamics of one complex variable 
which studies the properties of the iterated map 

X - I{x) ^ r\x) : = fU{x)) ^ . . . ^ r{x) ^ . . . , (4) 

and the algebraic properties of the (locally) analytic function f{x) = a^a;'^, considered as an element of the 
(infinite-dimensional) space M of all analytic functions C ^ C on the complex plane C. We attribute the 
apparent complexity of the Julia and Mandelbrot sets, associated with the pattern of orbit bifurcations, to the 
properties of well defined (but not yet well-studied) algebraic subspaces in moduli space: discriminant variety 
V and spaces of shifted n-iterated maps Mon 

F^{x-J) = n{x)-x. (5) 

In other words, we claim that Julia and Mandelbrot sets can be defined as objects in algebraic geometry and 
studied by theoretical methods, not only by computer experiments. 

We remind, that the boundary of Julia set dJ{f) for a given / is the union of all unstable periodic orbits 
of / in C, i.e. a subset in the union 0{f ) off all periodic orbits, perhaps, complex, which, in turn, is nothing 
but the set of all roots of the iterated maps Fn{x; f) for all n. The boundary of algebraic Mandelbrot set 
dM{n) C n C M consists of all functions / from a given family ii C M, where the stability of orbits changes. 
However, there is no one-to-one correspondence between Mandelbrot sets and bifurcations of Julia sets, since 
the latter can also occur when something happens to unstable orbits and some of such events do not necessarily 
involve stable orbits, controlled by the Mandelbrot set. In order to get rid of the reference to somewhat subtle 
notion of stability (which also involves real instead of complex- analytic constraints), we suggest to reveal the 
algebraic structure of Julia and Mandelbrot sets and then use it is a definition of their algebraic counterparts. It 
is obvious that boundaries of Mandelbrot sets are related to discriminants, dM = 1) <Z M and can actually be 
defined without reference to stability. The same is true also for Julia sets, which can be defined in terms of grand 
orbits of /, which characterize not only the images of points after the action of /, but also their pre-images. The 
bifurcations of so defined algebraic Julia sets occur at the boundary of what we call the grand Mandelbrot set, 
a straightforward extension of Mandelbrot set, also representable as discriminant variety. Its boundary consists 
of points in the space of maps {moduli space), where the structure of periodic orbits (not obligatory stable) and 
their pre-orbits changes. Hypothetically algebraic Julia sets coincide with the ordinary ones (at least for some 
appropriate choice of stability criterium). 

To investigate the algebraic structure of Julia and Mandelbrot sets we suggest to reformulate the problem 
in terms of representation theory. 

An analytic function as a map / : C — s- C defines an action of Z on C: n / o . . . o /. For / with 

n 

real coefficients this action has two invariant subsets R and C — R. Often in discrete dynamics one considers 
only the action on R. However, some peculiar properties of this action (like period doubling, stability etc.) 
become transparent only if we consider it as a special case of the action of / with complex coefficients and no 
distinguished subsets (like R) in C. 

The set of periodic orbits of this action can be considered as a representation of Abelian group Z, which 
depends on the shape of /. One can study these representations, starting from generic representation, which 
corresponds to generic position of / in the space M . Other representations result from merging of the orbits of 
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generic representation, what happens when / belongs to certain discriminant subsets = {V O Ain)* C Ai. 
These discriminant sets are algebraic and have grading, related to the order of the corresponding merging orbits. 
The boundary of Mandelbrot set 9M(/i) for a one-parametric family of maps fi C A4 is related to the union 
of the countable family of algebraic sets of increasing degree, UJ^q {V n Hon), and this results in the "fractal" 
structure of Mandelbrot set. This approach allows us to introduce the notion of discriminant in the (infinite- 
dimensional) space A4 of analytic functions as the inverse limit of a sequence of finite-dimensional expressions 
corresponding to discriminants in the spaces of polynomials with increasing degree. 

On the preimages of periodic orbits (on bounded grand orbits) the above Z-action is not invertible. The 
structure of the set of preimages of a given order k has special bifurcations which are not expressible in terms 
of bifurcations of periodic orbits. This implies the existence of a hierarchy of discriminant sets in the space of 
coefhcients (secondary Mandelbrot sets). The union of secondary Mandelbrot sets for all k presumably defines 
all the bifurcations of Julia sets. We call this union the Grand Mandelbrot set. 

Mentioned above is just one of many possible " definitions" of Julia and Mandelbrot sets [HII ■ Their interesting 
(fractal) properties look universal and not sensitive to particular way they are introduced, what clearly implies 
that some canonical universal structure stands behind. The claim of the present paper is that this canonical 
structure is algebraic and is nothing but the intersection of iterated maps with the universal discriminant and 
resultant varieties T) C Ai and TZ C M x M, i.e. the appropriately defined pull-back TZ* C A4. Different 
definitions of Mandelbrot sets are just different vies/projections/sections of TZ* . The study of Mandelbrot sets 
is actually the study of TZ* , which is a complicated but well-defined problem of algebraic geometry. Following 
this line of thinking we suggest to substitute intuitive notions of Julia and Mandelbrot sets by their much 
better defined algebraic counterparts. The details of these definitions can require modifications when one starts 
proving, say, existence- or equivalence theorems, but the advantage of such approach is the very possibility to 
formulate and prove theorems. 

Our main conclusions about the structure of Mandelbrot and Julia sets are collected in section 13 There 
are no properties of these fractal sets which could not be discovered and explained by pure algebraic methods. 
However, the adequate part of algebraic theory is not enough developed and requires more attention and work. 
In the absence of developed theory - and even developed language - a lot of facts can be better understood 
through examples than through general theorems. We present both parts of the story - theory and examples 
- in parallel.^ For the sake of convenience we begin in the next section from listing the main relevant notions 
and their inter-relations. 



^To investigate examples, we used powerful discriminant and resultant facilities of MAPLE and a wonderful Fractal Explorer 
(FE) program '11'. In particular, the pictures of Mandelbrot and Julia sets for 1-parameter families of maps in this paper are 
generated with the help of this program. 
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2 Notions and notation 



2.1 Objects, associated with the space X 

• Phase space X. In principle, in order to define the Mandelbrot space only X only needs to be a 
topological space. However, to make powerful algebraic machinery applicable, X should be an algebraically 
closed field: the complex plane C suits all needs, while the real line R does not. We usually assume that 
X = C, but algebraic construction is easily extendable to p-adic fields and to Galois field. In fact, further 
generalizations are straightforward: what is really important is the ring structure (to define maps from X into 
itself as polynomials and series) and a kind of Bezout theorem (allowing to parameterize the maps by collections 
of their roots). ^ 

• A map / : X ^ X can be represented as finite (polynomial) or infinite series 

f{x)=J2'^kx''. (6) 

fe=0 

The series can be assumed convergent at least in some domain of X (if the notion of convergency is defined on 
X). Some statements in this paper are extendable even to formal series. 

• /-orbit 0{x; /) C X of a point x € X is a set of all images, 

{xj{x),r^x),...,mx),...}. 

It is convenient to agree that for n = /°°(a;) = x. The orbit is periodic of order n if n is the smallest positive 
integer for which /°"(a;) = x. 

• Grand /-orbit GO{x; /) C X of x includes also all the pre-images: all points x' , such that f°^{x') — x 
for some k. There can be many different x' for given x and k. The orhit is finite, if for some fc > and n > 1 
/"("+'=) (a;) = f°^{x). If k and n are the smallest for which this property holds, we say that x belongs to the 
fc-th preimage of a periodic orbit of order n (which consists of the n points {/"'^(cc), . . . , /°("+'^^^) (a;)}). The 
corresponding grand orhit (i.e. the one which ends up in a finite orbit) is called bounded grand orbit (BGO), 
it looks like a graph, with infinite trees attached to a loop of length n. Grand orbit defines the ramification 
structure of functional-inverse of the map / and of associated Riemann surface. Different branches^ of the 
grand-orbit tree define different branches of the prepotential (the discrete analogue of T in eq.©). 

• The union £'„(/) C X of all periodic /-orbits of order n, the union On^kif) C X of all their fc-preimages 
and the unions 0{f) = U^^iO„(/) C X of aU periodic orbits and gO{f) = ^t=o,n=i^n,kU) C X of aU bounded 
grand orbits. The sets 0{f ) C GO{f ) can be smaller than X (since there are also non-periodic orbits and thus 
unbounded grand orbits) . 

• The sets Sn{f ) C X and Sn.kif) C X of all roots of the functions 
Fn{x; f) = n(x) - X (7) 

■'In order to preserve these properties, multidimensional discrete dynamics can be introduced on phase space X™ in the following 
way. The x- variables are substituted by m-componcnt vectors xi, . . . , Xm (or even m + 1-component if the space is CP™ and afHne 
parameterization is used). The relevant maps / ; X"^ X'" are defined by tensor coefficients a^.j:, 

oc 

fi{x)= ^ . . .x^", i = l,...,m 

fcl,...,fcm=0 

instead of Then one can introduce and study discriminants, resolvents, Mandelbrot and Julia sets just in the same way as we are 
going to do in one-dimensional situation (m = 1). For relevant generalization of discriminant see 1121 . Concrete multidimensional 
examples should be examined by this technique, but this is beyond the scope of the present paper. 

^For a rooted tree we call branch any path, connecting the root and some highest-level vertex. For infinite trees, like most 
grand orbits, the branches have infinite length. 
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'■ills' 



Figure 1: Typical view of a bounded grand orbit of order n = 2 (thus the loop is of length 2) for f{x) which is a quadratic 
polynomial (thus two arrows enter every vertex), (a) Example of generic case: f{x) = — 3, shown together with its embedding 
into the complex- a; plane, (b) The same example, f{x) = a;^ — 3, only internal tree structure is shown, (c) Example of degenerate 
case: f{x) = — 1 (only one arrow enters —1, its second preimage is infinitesimally close to 0). 




Figure 2; A non-degenerate bounded grand orbit of order n = 5 for cubic f(x). Only internal structure is shown, without 
embedding into X. 
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and 



Fn.ki^; f) - r^^+'Hx) - r\x) = F„+kix; f) - Fk{x; f) = F„ (x; f"") . (8) 

All periodic orbits Om{f) C whenever m is a divisor of n (in particular, when m = n). Similarly, 

Om,k{f) C Sn,k{f)- This property allows to substitute the study of (bounded) periodic (grand) orbits by the 
study of roots of maps F„_fc(x;/). Reshufflings (mergings and decompositions) of (grand) orbits occur when 
the roots coincide, i.e. when the map becomes degenerate. This property allows to substitute the study of the 
phase space X (where the orbits and roots are living) by the study of the space of functions, where the 
varieties of degenerate functions are known as discriminants T). 



• On{f) C X is a set of zeroes of function G'„(x; /), which is an irreducible-over-any-field divisor of reducible 
function Fn(x] /), 

Fn{x-J)^X{GUx-J). (9) 

m|n 

where the product is over all divisors m of n, including m = 1 and m = n. The functions Gn are more adequate 
for our purposes than Fn, but the Fn are much easier to define and deal with in sample calculations. 



• Stable periodic orbit of order n is defined by conditions 

|(r")'(a;)| < 1 i.e. |i^^(x) + 1| < 1 



G„(.) = 

for all its points. This definition is self-consistent because 

(r")'(^) = n /'(/°'(^)) = n /'W' (11) 

k—O js^^orbit 

i.e. the l.h.s. is actually independent of the choice of x E On{f)- Other periodic orbits are called unstable. 
Bounded grand orbits are called stable and unstable when periodic orbits at their ends arc stable and unstable 
respectively. 



• The set of stable periodic orbits C+(/) and its complement: the state of unstable periodic orbits 0^{f). 
By definition 0+{f) U 0_(/) = 0{f). 

• Julia set J(/) G X is attraction domain of stable periodic orbits of /, i.e. the set of points a; G X with 
orbits 0{x;f), approaching the set of stable orbits: Ve > 3fc and 3zk G ©+(/): |/°'^(a:) — Zk\ < e. This 
definition refers to non-algebraic structures like continuity. However, the boundary dJ{f) is an almost algebraic 
object, since it actually consists of all the unstable periodic orbits and their grand orbits, which are dense in 
dJif): 

dJif) ^o^if) = go^T) (12) 

Moreover, each grand orbit "originates in the vicinity" of 0-{f): generically, every branch of the grand-orbit 
tree is associated with - "originates at" - a particular orbit from 0_(/), and when C+(/) = this is a one- 
to-one correspondence. If 0+(/) ^ 0, the "future" of the grand orbit can be three-fold: it can terminate in 
a periodic orbit, stable or unstable, approach (tend to) some stable orbit or go to infinity (which can actually 
be considered as a "reservoir" of additional stable orbits). If C+(/) = 0, the situation is different: strange 
attractors can also occur. Description of dJ{f) in terms of orbits works even when 0+{f) — 0, but then it is 
not a boundary of anything, the Julia set J(/) itself is not defined. 

For a better description of Julia sets see summary in sIHland s l5. 51 below. 



10 



2.2 Objects, associated with the space M. 

• The moduli space M of all maps. In principle it can be as big as the set of formal series 0, i.e. can be 

considered as the space of the coefficients {a^}. Some pieces of theory, however, require additional structures, 
accordingly one can restrict M to sets of continuous, smooth, locally analytic or any other convenient sets of 
maps. We assume that the maps are single- valued, and their functional inverses can be described in terms 
of trees, perhaps, of infinite valence. The safest (but clearly over-restricted) choice is a subspace V C A4 oi 
all polynomials with complex coefficients, which, in turn, can be decomposed into spaces of polynomials of 
particular degrees d: V = W^Vd- One can also consider much smaller subspaces: families fi{c) C M, with ak{c) 
in © parameterized by several complex parameters c = (ci, C2, . . .). Most examples in the literature deal with 
one-parametric families. 

• The moduli subspace Mon C Al of shifted n-iterated maps, consisting of all maps F : X — > X which 
have the form lO for some f{x). The shift by x is important. Eq.(|Zj) defines canonical mappings 

/„ : M — > Mon, 

which can be also written in terms of the coefficients: for f{x) parameterized as in eq.®, and 

oo 

Fr.{x-J)^n{x)-x = Y,ai"^x'^ 

k=0 

we can define /„ as an algebraic map 

L : {ak} {4"^} (13) 

where all a^"-* are polynomials of {ai}.^ Inverse map 1^^ is defined on Mon only, not on entire M, but is not 
single- valued. In what follows we denote the pull-back of /„ by double star, thus [Inf)** = f (while I^^{Inf) 
can contain a set of other functions in addition to /, all mapped into one and the same point /„/ of Mon)- 
The map /„ can be restricted to polynomials of given degree, then 

4 : Vd—* Vd^ (14) 

embeds the {d + l)c-dimensional space into the (d" -I- l)c one as an algebraic variety.^ Like in (|14() . often 
Mon = Inip) 't- One can, however, use a pull-back to cure this situation: \B n /^on]** C /i for any B d M. 

In the studies of grand orbits more general moduli subspaces are involved: Mo{n,k) C M, consisting of all 
maps of the form © for some /. 

• Similarly, the functions Gn{x; f) define another canonical set of mappings, 

Jn- M — > Mn- 

into subspaces Mn C M, consisting of all functions which have the form Gn{x; f) for some /. These are, 
however, somewhat less explicit varieties than Mom because such are the functions Gn{x; f) - defined as 
irreducible constituents of Fn{x; /). Still, they are quite explicit, say, when / are polynomials of definite degree, 
and they are more adequate to describe the structures, relevant for discrete dynamics. The J-puUback will be 
denoted by a single star. 

^In order to avoid confusion we state explicitly that a^"' = and aj.^' = aj, — Also Moi = M, but in general analogous 

relation is not true for subsets fi C M: fiol does not obligatory coincide with fi, it differs from it by a shift by x. 
^For example, (7^2)02 is embedded into ■P4 as algebraic variety: 

SA^{D + 1) - AABC + = 0, 

QAA^{AAB'^E + ~ iA^{D + 1)2)3 = B^{1QA^{D + 1)- B^f 

This follows from 

Ax'^ + Bx^ + Cx^ + Dx + E = {ax^ + bx + c)°^ - x = 
= a^x* + 2a'^bx^ + {ah^ + ah + 2a^c)x^ + {2ahc + fe^ - l)x + (ac^ + hc + c). 
Note, that because of the shift by x one should be careful with projectivization of this variety. 
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Figure 3: Schematic view of inter-relations between the subsets Mon, Mon and 'Von- 

• Discriminant variety T) C M consists of all degenerate maps /, i.e. such that f{x) and its derivative 
f'{x) have a common zero. It is defined by the equation 

in A4, where D{f) is the square of the Van-der-Monde product of roots differences: 

for fix) ^ l[{x - rk)D{f) ^ \{{rk - n). 

k k^l 

Restrictions of discriminant onto the spaces of polynomials of given degree, ViVd) = T) C\ Vd, are algebraic 
varieties in Vd, because D{f) is a polynomial (of degree 2c? — 1) of the coefficients {a^} in representation © 
of / (while roots themselves are not polynomial in {ofc} the square of Van-der-Monde product is, basically, 
by Vieta's theorem). Discriminant variety V has singularities of various codimensions k in V, associated with 
mergings of fc -I- 2 roots of /. 

• Resultant variety TZ C A4 x M consists of pairs of maps f{x), g{x) which have common zero. It is defined 
by the equation 

^(/,5)-0, 

and 

for f{x) ^ JJ(x - rk) and g{x) ^ ]^(a; - s/) 

k I 

R{f,g)-l[{rk-si). 

k,l 

Again, if / and g are polynomials, resultant is a polynomial in coefficients of / and g. Resultant variety 
has singularities when more that one pair of roots coincide. Discriminant can be considered as appropriate 
resultants' derivative at diagonal f — g. Higher resultants varieties m Ai x . . . x A4 are also important for our 
purposes, but will not be considered in the present paper. 

• Intersections of 2? and TZ with subspaces like Mon and Mn, i-e. Von = VnMon, TZo(m,n) = T^f^ {-Mom x 
■A^on); T^n = Vr] Mn, T^m,n = TZ C] {Mm X Mm) cousist of disconnected components (see Fig|2J. /-puUbacks 
'^on, Klm,n) ^ and J-pullbacks 

!?: = {/: DiGnif)) = 0}cM, 

n*m^n = {/ : RiGmif), G„(/)) = 0} C M 

also consist of numerous components, disjoint and touching. Particular families ^ G M intersect these singular 
varieties and provide particular sections dM (/i) of this generic structure: the universal discriminant (or 
resultant) variety 

cx) C30 oo 

v* = n* = \Jnn^ U n*^,n U 

n— 1 n,m— 1 n— 1 
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which is the boundary of the Universal Mandelbrot set dM{A4). 



• For a given family of maps fi C A4 the boundary of algebraic Mandelbrot set dM{p) C /i is defined 
as a J-puUback 

aM(M) = r^^i)^{vr^ (u- i^™))* = 

= A*U(U-i{/: D{GM)) = Q}) (15) 



• Varieties I?**, 7?,**, I?*, 7?,* - and thus the boundary of algebraic Mandelbrot space dM - can be instead 
considered as pure topological objects in M. independent of any additional algebraic structures, needed to define 
Fm G„, discriminant and resultant varieties T) and TZ. For example, V^* — Vl consists of all maps f E M with 
two coincident fixed points. Higher components P* C I?** consist of maps /, with two coincident fixed points 
of their n-th iteration (for V**) and such that any lower iteration does not have coincident points (for I?*). 
Resultants consist of maps / with coinciding fixed points of their different iterations. Fixed points and iterated 
maps are pure categorial notions, while to define coincident points one can make use of topological structure: 
two different fixed points merge under continuous deformation of / in A^. 

• Stability domain of the periodic order-n orbits 5'„ C is defined by the system (|1(J|I : if a root x of 
Gn{f) is substituted in to inequality, it becomes a restriction on the shape of /, i.e. defines a domain in A4. 
This domain is highly singular and disconnected. For a family ^ C A4 we get a section Sn{^) = 5„ n /x, also 
singular and disconnected. 

• Mandelbrot set M{A4) C is a union of all stability domains with different n: 

M = Ujf^i5„, Mip) = U^=i5„(Ai) = M n fx. (16) 

Definitions (|16|l and (|15(l leave obscure most of the structure of Mandelbrot set and its boundary. Moreover, 
even consistency of these two definitions is not obvious.^ For better description of Mandelbrot set see s|3|below. 

• Above definitions, at least the algebraic part of the story, can be straightforwardly extended from orbits 
to grand orbits and from maps Fn{f) and (?„(/) to Fnk{f) and Gn,k{f)- The maps Gn,k{x', f), which define 
fc-th pre-orbits of periodic order-n orbits are not obligatory irreducible constituents of Fn^kix] f) (what is the 
case for fc = 0). They are instead related to peculiar map 

z: X^X 

associating with every f E Ai the points a; e X with degenerate pre-image: 

Zf^{z: D{f{x)^z)^0}. 

Generalizations of Mandelbrot sets to fc > 1 are called secondary Mandelbrot sets. The Grand Mandelbrot 
set is the union of secondary sets with all k > 0. Bifurcations of Julia sets with the variation of / inside fi d M. 
are captured by the structure of Grand Mandelbrot set. 

We see that, the future theory of dynamical systems should include the study of two purely algebro-geometric 
objects: 

(i) the universal discriminant and resultant varieties V <Z M and TZ C Ai x A4, 
(n) the map /x — > ji = UJ^;^^^ and the pull-backs V* C M and TZ* C M. 

It should investigate intersections of 2? and p, and TZ with jl x jl and consider further generalizations: to 
multi-dimensional phase spaces X and to continuous iteration numbers n. It is also interesting to understand 
how far one can move with description of TZ* and universal Mandelbrot space by pure topological methods, 
without the ring and other auxiliary algebraic structures. 

*It deserves emphasizing that from algebraic perspective there is an essential difference between the Mandelbrot and Julia sets 
themselves and their boundaries. Boundaries are pure algebraic (or, alternatively, pure topological) objects, while entire spaces 
depend on stability criteria, which, for example, break complex analyticity and other nice properties present in the description of 
the boundaries. 
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Figure 4: Divisors tree t[2A]. It contains divisor trees for all its divisors: t[12], t[6], t[4], t[3], t[2], t[l], some in several copies. 

\ \ 




Figure 5: A piece of the forest of multipliers trees. Trails are shown by "steps". 



2.3 Combinatorial objects 

• Divisors tree t[n] is a finite rooted tree. Number n stands at the root, and it is connected to all T(n) — 1 
of its divisors k\n (including fc = 1, but excluding k = n). Each vertex k is further connected to all t(A;) — 1 of 
its divisors and so on. The links are labeled by ratios m = n/k > 2 and = ki/ki+i > 2. Number n is equal 
to product of m's along every branch. See examples in Fig01 The number T{n) of divisors of n = . . .p'^'' is 
equal to T{n) ~ (ai + 1) . . . (ofe + 1). A generating function (Dirchlet function) has asymptotic 

D{x) = ^T{n) ^ xlogx + {2C - 1) + 0{V^), (17) 

C - Euler constant. 

• Multipliers tree r„ is a rooted tree of infinite valence. At the root stands the number n. The branches 
at every level p are labeled by positive integers = 1,2,3,4, . . .. Every vertex at p-th level is connected to 
the root by a path mi, m2, . . . , mp and the product nmim2 ■ ■ ■ mp stands in it. The basic forest T = U„=iT„. 
A number of times B[n) the number n occurs in the forest T is equal to the number of branches in its divisors 
tree. 
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Fi gure 6; A fragment of the basic graph. Only the generating arrows, associated with primes, are shown (to avoid overloading 
the picture). All compositions of arrows should be added as separate arrows (for example, arrows connect 1 to 4, 6, 8, 9 and all 
other integers). 



The numbers Tp{n) of ways to represent n as a product of p integers, i.e. the numbers of times n appears at 
the p — th level of multipliers tree, are described by the generating function 



C('S) = Sm=i ™ is Riemann's zeta-function. Note that H17|l describes r(n) = T2{n). 

• Basic graph B is obtained from the multipliers tree Ti by identification of all vertices with the same 
numbers, so that the vertices of B are in one-to-one correspondence with all natural numbers. Since Ti was a 
rooted tree, the graph B is directed: all links are arrows. There are rin) — 1 arrows entering the vertex n and 
infinitely many arrows which exit it and lead to points mn with all natural m. 

2.4 Relations between the notions 

Relations are summarized in the following table. 
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objects in X 



objects in M 



/ 



unstable periodic orbits = 

= 9(Julia set ,/(/)) = 
^(algebraic Julia set JaU)) 
I 

reshuffling of JA{f) 



map / 

I 

iterated maps Fn{f) 

I 

I 

n — periodic / — orbits = 
= roots of Gn{f) 

I 

coincident roots of Gn{f) 
= reshuffling of orbits 



maps Fn^kif) 

bounded grand / — orbits 

= roots of F„^k{f) 
infinite preimage of BGO 

I 

coincident roots of -Fn,fe(/) 
resuffling of BGO = 



moduli space M 
U 

spaces Mon 
spaces jVlfc 



{D{G„{/)} = 0} 
discriminant varieties = 
= ^(Mandelbrot set M) 



spaces Mo(„,k) 



= {D{F„,kU)} = 0} 
9(grand Mandelbrot set) 
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3 Summary 



In this section we briefly summarize our main claims, concerning the structure of Mandelbrot and Julia sets. 
They are naturally splited in three topics. 

3.1 Orbits and grand orbits 

This part of the story includes: 

- The theory of orbit and pre-orbit functions /) and G'„^fe(a;; /). 

- Classification of orbits, pre-orbits and grand orbits. 

- Intersections (bifurcations) of orbits and degenerations of grand orbits. 

- Discriminant and resultant analysis, reduced discriminants dn and resultants rm,n intersections of resultant 
and discriminant varieties. 

All these subjects - to different level of depth - are considered in s0] below and illustrated by examples in 
sslHlandEl Systematic theory is still lacking. 

3.2 Mandelbrot sets 

Figs|7|a-b show three hierarchical levels of Mandelbrot set: the structure of individual component {Mi in 
Fig[7|a); existence of infinitely many other components (like M3.Q in Fig|7|b), which look practically the same 
after appropriate rescaling; trails connecting different components Mka with Mi, which are densely populated 
by other M^p (well seen in the same FigI3b). These structures (except for the self-similarity property) are 
universal: they reflect the structure of the universal Mandelbrot set, of which the boundary is the universal 
resultant variety TZ* . This means that they do not depend on particular choice of the family /i (above figures 
are drawn for /c = 2;* + c, which is in no way distinguished, as is obvious from similar pictures for some other 
families, presented in ssElandEl). The only characteristics which depend on the family are the sets of indices a, 
labeling degeneracies: different domains with the same place in algebraic structures. Algebraic structure reflects 
intersection properties of different subvarieties 7?,J^„, while a's parameterize the intersection subvarieties. To 
classify the a-parameters one should introduce and study higher resultants and multiparametric families fi d M, 
what is straightforward, but (except for introductory examples) is left beyond the scope of the present paper. 

The main property of the resultant varieties 7?.*j „, responsible for the structure of Mandelbrot set, is that 
they are non-trivial in real codimension two only if n is divisible by m or vice versa. All the rest follows from 
above-described relation between the resultant varieties and the boundary of Mandelbrot set. Stability domains 
Sn and Mandelbrot sets M are made out of the same disk-like^ building blocks, but connect these blocks in two 
different ways. Mandelbrot set and its "forest with trails" structure reflect the universal (/i-independent) labeling 
(ordering) of these elementary domains (denoted by a below), while their sizes, locations are self-similarity 
properties (which depend on the choice of ^i) are dictated by stability equations. 

3.2.1 Forest structure 

For any family of maps /i Mandelbrot set M{fj,) has the following hierarchical structure, see Fig[7|a (for particular 
(U some components - or, better to say, their intersections with fi - can be empty): 

M{n) consists of infinitely many disconnected components (of which only one is clearly seen in Fig[7|a and 
some others are revealed by zooming in Figdb), 

M {p) = U Mkc.ifJ-),Mkc. n Mii3 =0, a k^l OT P, 

labeled by natural number k and an extra index a, belonging to the /^-dependent set Vk{n). For polynomial / 
the set is finite, its size will be denoted by |z/fc(/i)|. 

Every component MkaifJ-) is a union of elementary domains, each with topology of a disc (interior of a 
{2D — l)R-dimensional sphere if /i is a Z?c-parametric family), which form a tree-like structure. Vertices of the 
tree Tka - the skeleton of Mka ~ are associated with "centers" of elementary domains, a link connects two 

^In the case of multi-parametric families fi C A4, when dimension of components of Mandelbrot set l\I{fi) is equal to dim fi , 
"disc" may be not a very adequate name. In fact, these components look more like interiors of cylinders and tori rather than balls. 
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Fi gure 7'. a) General view of the Mandelbrot set for the family fc{^) = + c. Actually, only the Mi component is well 
seen. Arrow points at another component, Ma.^, shown unlarged in b). The structure of the trail between Mi and Ms,^ (almost 
unobservable in this picture) is also seen in b) . Zs-symmetry of the picture is due to invariance of equations like + c = a:: under 
the transformation x —^ e'^'^^^^x, c — > e^'^'^^c. Associated Julia set has Z4 symmetry, see Figl^l b) Enlarged component M^^a, 
included in the box in a). It looks like exact copy of Mi in a). In fact similarity is approximate, but the deviations are damped by 
the ratio of sizes, r^/ri <^ 1. In this unlarged picture also the trail is seen between Mi and M^^a and other Mf^^ component in it. 



vertices whenever the two corresponding domains have a common point - and there can be at most one such 
point, where the two domains actuaUy touch each other. The skeleton is a rooted tree, moreover, it is actuaUy 
a covering of the multiphers tree Tk, introduced in s l2.3l The difference is that every branch at the level p 



carrying the label rup has multiplicity |^'mp(/i)| and thus an extra label ap G Vmlil^)- The union of skeleton 
trees is the skeleton of the Mandelbrot set - the Mandelbrot forest 

T{M) = |jTfe„. 



(p)/ 



Thus 



M. 



and, finally, M (n) = M D n, 



k.OL 



Mk. = U ^fc^ 

p=0 



(p) _ 



u - 

111 ,--.,mp = l 

ai ,...,ap 



(p) 



Ull 

ai 



00 ( 00 

M = uu U 

k,ap=0 \ mi,. .,mp = l 
\ Qi,...,ap 



mi 
ai 



The elementary domains are ct^^-* 



mi 
ai 



ka 



they are connected by a link in the powerful tree Tka, i.e. when p' 
i=l,...,p. 



and two of them touch at exactly one point, whenever 
p + I, and m,^ = m^ , = for all 



3.2.2 Relation to resultants and discriminants 

The touching point (map / € /i) belongs to the resultant variety 7?.* „, , i.e. is a root of i?(G„(/), G„'(/)) — 0, 



with n — fcmi , 



. TUp and n' 



fcmi . 



'p+i- 



The boundaries of all elementary domains cr^^' with 



p > are smooth (unless fj, crosses a singularity of TZ), while Scr'-*'-' [fca] has a cusp, located at discriminant 
variety I?^, i.e. at / € /i, satisfying D{Gk{f )) = (o-^*^^ [/ca] itself has a peculiar cardioid form). It is clear that 
the bigger the family /j,, the more intersections it has with the resultant and discriminant varieties, thus the 
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bigger are the sizes of the sets v^P^di). For entire Ai the indices ap get continuous and parameterize the entire 
puUbacks of the resultant and discriminant varieties. 

Every elementary domain a^^^ touches a single domain of the lower level p — 1 and infinitely many domains 
of the next level p+ 1 (which are labeled by all integer m = mp+i > 2 and a — ap+i € i'm'^^\fJ-)) ■ The touching 
points - belonging to zeroes of i?(G„,G'„,„) with all possible m — mp+i are actually dense in the boundary 
da'^P^ i.e. the boundary can be considered as a closure of the sets of zeroes. 



3.2.3 Relation to stability domains 



As already mentioned, every ct'^' 



ka 



is characterized by an integer n = kmi . . . rUp. Some- 



mi . . . nip 
ai ... Q!p 

time, when other parameters are inessential, we even denote it by cr^f''. Stability domain SnifJ-), defined by the 
conditions (|10|l is a union of all (t(/i) with the same n. This implies, that Sn is actually a sum over all branches 
of divisor tree t[n\, introduced in s l2.3l Then p is the length of the branch at k is its end-link. Other links carry 



numbers mi, 



The only new ingredient is addition of extra parameters a at every step. This means that 



we actually need a-decorated divisor trees (for every family /z), which we denote i[n](/i) and imply that links 
at p-th level carry pairs nip, ap with ap £ VmlilJ-), and sum over decorated trees imply summation over a's. In 
other words, the elementary domains a are actually labeled by branches of the decorated trees, Bt, and 



Sn= [j an[Bi] 



Bth 



( 



u u 

p=0 



u 



mi 
ai 



ka 



(18) 



It remains an interesting problem to prove directly that eas. (|10() imply decomposition (|18|l . 

3.2.4 Critical points and locations of elementary domains 

Let Wf denote critical points of the map f{x): 

f{wf) = 0. 

Then (/°")'('w^/) = 0, i.e. the critical points always belong to some stable orbit, see H10() . This orbit is periodic 
of order n provided 



Gn{wf;f)^0. 



(19) 



This is an equation on the map / and its solutions define points / in the family /i C AI, which lie inside 



stability domain Sn, and, actually, inside certain elementary domains cr^P^ 



. In fact this 



mi . . . mp 
ai ... ttp 

is a one-to-one correspondence: solutions to ea. (|19|l enumerate all the elementary domains. 

This implies the following pure algebraic description of the powerful forest of M(/i) (formed by the coverings 
Tka of the multipliers trees Tfc). All vertices of the graph are labeled by solutions of H19I) and accordingly carry 
indices n. Links are labeled by roots of the resultants R{Gm{f),Gn{f)), with appropriate m and n (standing 
at the ends of the link) . a-parameters serve to enumerate different solutions of (|19|) and different roots of the 
same resultants. 

Critical points Wf also play a special role in the study of bifurcations of grand orbits and Julia sets, see s l4.6l 



3.2.5 Perturbation theory and approximate self-similarity of Mandelbrot set 

For concrete families ^ G M a, kind of approximated perturbation theory can be developed in the vicinity of 
solutions / = /o to (|19|l . Namely, one can expand equations 110|l . which define the shape of stability domain 
(and thus of the elementary domains CT[/ca]), in the small vicinity of the point {x, f) = {wf^, /o) and, assuming 
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that the deviation is small, substitute original exact equations by their approximation for small deviations. 
Approximate equations have a universal form, depending only on the symmetry of the problem. This explains 
why in many examples (where this method is accurate enough) all the components Mka of Mandelbrot set look 
approximately the same, i.e. why the Mandelbrot set is approximately self-similar (see Fig|Z|b) and help to 
classify the components of stability domain. The same method can be used to investigate the shape (not just 
structure) of Julia sets. See in sec l4.9.4l below an example of application of this procedure to the map families 
fc = x'' + c. 



3.2.6 Trails in the forest 

The last level of hierarchy in the structure of Mandelbrot set is represented by trails, see Fig[7|b. Despite the 
components Mka do not intersect (do not have common points), they are linked by a tree-like system of trails: 
each Mka is connected to Mi by a single trail Tfca, which is densely populated by some other components Mip. 
The trail structure is exhaustively described by triangle embedding matrix, r = {T[ka,lf3]} with unit entry 
when one trail is inside another, ti/s C Tka, and zero entry otherwise. If locations of all Mka in parameter-c 
space are known (for example, evaluated with the help of the perturbation theory from s l3.2.5(l . then embedding 
matrix fully describes the trails. It is unclear whether any equations in parameter space /i can be written, which 
define the shape of trails. 

Embedding matrix seems to be universal, i.e. does not depend on the choice of the family /x. If true, this 
means that the trail structure is indeed a pertinent characteristic of Mandelbrot /resultant variety, not of its 
section by fi. However, as usual, the universal structure is partly hidden because of the presence of non-universal 
degeneracies, labeled by a-parameters (at our level of consideration, ignoring higher resultants): concrete trail 
systems T(/i) are coverings of presumably- universal system. 

Trail structure requires further investigation and does not get much attention in the present paper. 



3.3 Sheaf of Julia sets over moduli space 

The structure of Julia set J{f) is also hierarchical and it depends on the position of the map / in Mandelbrot 
set, especially on its location in Mandelbrot forest, i.e. on the set of parameters {fca,p, m^ai}, labeling the 



elementary domain a 



which contains /. Since in this paper we do not classify a- 



ip) ™i ••■ "^P 
a\ ... Qfp 

dependencies, only the numbers fc (labeling connected component), p (labeling the level) and mi, . . . , mp will 
be interpreted in terms of J(/) structure. 

Julia set is a continuous deformation of a set of discs (balls), which - the set - depends on fc; increase of p by 
every unit causes gluing of infinitely many points at the boundary, at every point exactly m^ points are glued 
together. Every elementary domain is associated with a set of stable orbits (often there is just one), contained 
inside the Julia set together with their grand orbits. All other (i.e. unstable) periodic orbits and their grand 
orbits belong to the boundary of Julia set and almost each particular grand orbit fills this boundary densely. At 
the touching point between two adjacent elementary domains the stable orbit approaches the boundary dJ from 
inside J (perhaps, it is better to say that the boundary deforms and some of its points - by groups - approach 
the orbit which lies inside J), intersect some unstable orbit and "exchange stability" with it. Their grand orbits 
also intersect. The ratio of orders of these two orbits is rup and every nip points of unstable grand orbit merge 
with every point of the stable one, thus strapping the disc at infinitely many places (at all points of the merging 
grand orbits) and pushing its sectors into bubble-like shoots, see FiglHl As / transfers to the elementary domain 
at level p + 1, the new stable orbit (the one of bigger order) quits the boundary and immerses into Julia set, 
while the orbit of smaller order, which is now unstable, remains at its boundary: the singularities, created when 
the grand orbits crossed, can not disappear. 

One can say that every link of the forest describes a phase transition of Julia set between neighboring 
elementary domains. An order parameter of this transition is the distance between approaching stable and 
unstable orbits on one side, looking as the contact length between emerging components of Julia set, and the 
angle between these components on the other side, when they are already separated except for a single common 
point. Both quantities vanish at transition point (i.e. when / is at the touching point between two domains 
in the Mandelbrot set). If we move around the transition point in A4, there is a non-trivial monodromy (the 
Julia set gets twisted), but on the way one should obligatory pass through the complement of Mandelbrot set. 
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where the discs - which form JuUa set - disappear, only the boundary remains, JuUa set has no "body", only 
boundary, and exact definition of monodromy is hard to give. 

Juha sets can change structure not only when stable orbit intersects with unstable one - what happens when 
the corresponding elementary domains a of Mandelbrot set touch each other, - but also when unstable grand 
orbits at the boundary dJ{f) cross or degenerate. Such events are not reflected in the Mandelbrot set itself: 
a bigger Grand Mandelbrot set, which contains all zeroes of D{Fn.k) and R{Fn,k, Fm,i) on its boundary, 
should be introduced to capture all the bifurcations of Julia set. 

See s l5.5l for some more details about Julia sets and their bifurcations. 
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Figure 8: The typical forms of the JuHa set for the map /c = + c and c £ Mi. a) c in the center of tr^" (c = 0). Z4 symmetry 
of the picture is due to the symmetry of the map fc- Associated Mandelbrot set has symmetry Z3, see Fig|7|a. b) c at the touching 
point (tJ°' U o-^^' {c= ^21/3(1 + 4^3)). Infinitely many pairs of points on the unit circle (a) are identified to provide this pattern, 
c) c at another touching point ctj"' U CTj^' (c = g(9 + « lllx/S)^''^)- Infinitely many triples of points on the unit circle (a) are 

identified to provide this pattern, d) c inside CTj^' (c = 0.45 + iO.8). This is a further deformation of (a) in the direction of (b). 
Once appeared in (b) the merging points at the boundary and straps of the disc caused by this merging of boundary points do not 
disappear, but angles between the external bubbles increase from the zero value which they have in (b). They will turn into 2it/2 

and bridges become needles when c reaches the boundary of Cj^' , say an intersection point U <j^^^ (e) or cri^^ U <jj?^ (f ). e) c at 

the touching point cr!^^ U (c 0.56 + 0.97i). f) c at the touching point ai^^ U (c 0.309 — 0.862j). A mixture of merging 
pairs and triples is clearly seen. The bridges for pairs became needles, while triplets are touching. Angles between triplets will 
increase as c goes inside and reach 27r/3 at its boundary, and so on. Arrows point the positions of Julia sets shown in Figlslin 
the Mandelbrot set. The "Julia sheaf is obtained by "hanging" the corresponding Julia set over each point of the Mandelbrot set. 
At the boundary of Mandelbrot set the Julia sets are reshuffled and the task of the theory is to describe this entire variety (the 
sheaf) and all its properties, both for the universal Mandelbrot set and multi-dimensional Julia sets associated with multicomponent 
maps, and for the particular sections, like the one-parametric family of single-component quartic maps, f{x; c) = x'^ + c shown in 
this particular picture. 
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4 Fragments of theory 



4.1 Orbits and reduction theory of iterated maps 

In discrete dynamics trajectories are substituted by orbits and closed trajectories - by periodic orbits 0„ of the 
finite order n. Different orbits have no common points, and / acts on each 0„ by cyclic permutation of points. 
All points of every 0„ belong to the set 5„(/) of roots of the function 

F„{x;f) ^ n{x) - X. 

The map / generates the action of cyclic group f{n) on iS„(/), and this set is decomposed into orbits of /(n) 
of the orders k, which are divisors of n (all On with k = n are among them). If n is divisible by fc, then F„(x) 
is divisible by -F'fe(x), and in general -F„ decomposes into irreducible constituents: 

T(n) 

Fn{x-f)=\{Gk{x-f) (20) 

k\n 

with a product over all possible rin) divisors k oin (including k = \ and k = n). The number of periodic orbits 
of order n is equal to Nn{f)/n, where Nn{f) is the number of roots of Gn{x; /). Irreducibility means that Gk 
are not decomposable into simpler constituents in an X- and /-independent way (of course, when Gk{x) are 
polynomials they can be decomposed into products of monomials over C, but this decomposition will not take 
place over R). 

Among the tasks of the theory is the study of reducibility of the sets iS„(/), the ways they decompose into 
orbits and the ways this decomposition changes under the deformations of / within ^ <Z M.. Since the story is 
essentially about the roots of functions, it gets much more transparent in the complex setting than in the real 
one. Also, the entire theory is naturally a generalization from the case of the polynomial functions f E V C A4. 

Obviously, for f EVd CV, Fn € T^d", and Gk{x) in l^njl is a polynomial of certain degree Nk{d). If d = p is 
prime, then the group f{n) — Z„ and is isomorphic to the Galois group over Fp of the cyclotomic polynomial 

— X (while the Galois group of / itself over C is trivial). If instead n = pis prime, then d out of the roots 
of Fp{x) are also the roots of Fi{x), i.e. are invariant points (orbits of order one) of /, while the remaining 
dP — d roots decompose into ^-y^ orbits of order n = p. According to 120|l the numbers Nnid) for f € Vd can 
be defined recursively: from 

T(n) 

^A^fe(rf)=degF„ = d" 

k\n 

it follows, that 

iV„(d) = rf"-^7Vfe(d). 

k<n 

The lowest numbers Nn{d) are: 

Ni{d) = d, 
N2id) = d^ -d = d{d-l), 
Nsid) = d^-d={d- l)d{d + 1), 
Niid) ^d^ -d^ ^ d^{d^ - 1) = {d-l)d^{d+l), 

N^id) ^d^ -d, 
Ne{d) = (d - l)d{d +l){d^ + d- 1), 
Njid) = d'^ -d, 
Ns{d) = d'^ -d^ = d\d'^ -I), 
Ng{d) = d° - d^ = d^(d^ - 1), 
iVio(d) = rf(d-*-l)(d^ + d-l), 
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For n prime, n = p, Np{d) — dP ~ d — d{dP — 1), and small Fermat theorem guarantees that Np{d) is divisible 
by p. Further, 

Np.{d) = dP" ~ dP"'" = dP"'" (d(f-i)p'"' - l) , 
Np,p,{d) = dP'P^ ~ dP' - dP^ + d, 

in particular, 

N2p{d) = {dP - d){dP + d~l), 
Nzp[d) ^ {dP - d){d^P + dP+^ +d^- 1). 

4.2 Bifurcations and discriminants: from real to complex 

Usually considerable part of discrete dynamics concentrates on the real part of above story, i.e. deals with the 
functions /, which map the real line R C C into itself. Such functions form a subspace A4^^^ in A4. The 
subset 7^(^) c A^*^-^-* consists of all polynomials with real coefficients. In variance with the complex situation, 
different polynomials Fn in the same V^d can have different numbers of real roots, thus even the size of the 

set si^\f) can change when / is varying inside A^^^^. In fact, the orbits of / G A4^^'^ are either entirely 
real, O C R or entirely complex, O C C — R. In the latter case the orbit is either self-conjugate or there is a 
complex conjugate orbit O ^ O. The number of real roots of Fn{x), i.e. the size of the set Sn (/), can change 
whenever a self-conjugate orbit or a pair of conjugate orbits, that was complex, becomes real or vice versa, a 
real orbit turns into a complex self-adjoint or a pair of conjugate orbits. This is generalization of the well known 
phenomenon, when a pair of complex-conjugate roots of a polynomial from becomes real, or vice versa, 

when a pair of real roots collide and go away into the complex domain - just in the case of peculiar polynomials 
Fn{x) from Mon this happens at once with entire orbits, not just with pairs of roots. In the case of polynomial's 
roots, the points in V^^\ where some roots can migrate between the real line and complex domain, belong to 
discriminant varieties C Vd, consisting of polynomials with at least a pair of coincident roots (i.e. such, that 
P{x) and P'{x) have at least one common root).^° The same is true for the orbits of /: all kinds of reshuffling 
of orbits take place when Fn G "D C A4. Note, that in the last statement neither the discriminant varieties nor 
the maps (and their orbits) are restricted to polynomials with real coefficients, everything is equally well defined 
at least for entire space of complex polynomials V. Moreover, a generalization should exist to entire space M. 
of analytic (not obligatory polynomial) functions, see s l4. 101 below. 
Coefficients a{,"^ of 

oo 

i^„(a;) = r"(x)-x = ^ai")x'= 

k=0 

are polynomials of the coefficients ak of f{x) — X]fc°=o ^■kx'' , thus we have an algebraic map (|13|l between the 
spaces of coefficients /„ : {a} {a(")(a)}. Then the space A4 of coefficients {a} of / is divided into connected 
components by the /„-pullback of the discriminant variety Von- Given a path a{t) in the space M, its image 
In{a(t)) can cross Von = V C\ Mon at a generic (non-singular) point /„(a(to))- Fn{x;t) will have different 
number of real roots before and after the crossing (i.e. for t' < to < t" in a neighborhood of to)- This change in 

^"The simplest example of non-trivial discriminant variety is T>2 - the quadric ~ iac = in the space 1^2 of quadratic 
polynomials ax^ -\-bx -\- c. In accordance with 1531 below 



4ac = - i det ( 2a b 



b 
b 

2a 6 



If quadratic polynomial is considered as a quadric in CP-"-, Q(x, y) = ax^ + bxy + cy^ , discriminant is a determinant of 2 X 2 matrix 

f 2a b \ 
{ b 2c J 



with two lines formed by the coefficients of and : 

^ ox ay 



■ 4ac = — det I 
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the number of real roots is called bifurcation in the theory of dynamical systems. Let the number of real roots 
m' for t' < to be less then the number m" for t" > ta. If a new root Xi{t") belongs to an /-orbit O, then all the 
elements of O belong to the /-invariant domain C — R for t' < to and " simultaneously come" to the /-invariant 
domain R at t = to- Since elements of O and its conjugate O come to the real axis R simultaneously, one can 
say, that real roots "are born in pairs" at t = to- 

4.3 Discriminants and resultants for iterated maps 

General comments on the definitions of discriminants and resultants (including the non-polynomials case) are 
collected in s l4. Ill below. However, for iterated maps these quantities are highly reducible. 
Discriminant variety is defined by the equation 

DU) = 

where, for polynomial /, D{f) is the polynomial of the coefficients {a}, see eq.ljHU below. Similarly, D{Fn) is 
a polynomial of the coefficients {a^"^(a)} and the resultant R{Fn, Fm) is a polynomial of coefficients |a("^(a)} 
and {a(™)(a)}. 

Since function Fn{x) is reducible over any field (e.g. over R as well as over C), see (j2()|> . the resultant 
factorization rule H54|l implies: 

D{Fn)^l[D{Gk)l[R\Gk,Gi)- (21) 

k\n k,l\n 

k>l 

However, this is only the beginning of the story. Despite Gk{x) are irreducible constituents of Fn{x), the 
resultants R{Gk, Gk/m) and discriminants D{Gk) are still reducible: 

R{Gk,Gi) = r\Gk,Gi), l< k, 

n 

D{Gk)=d\Gk) n Ii"''\Gk,Gk/m) = d\Gk)Y[r''-'{Gk,Gi), 

m>l l\k 

thus D{Fr,)='[[id\Gk)Ylr''+'{Gk,Gi)]- (22) 

k\n \ l\k J 

Indeed, whenever a resultant vanishes, a root of some G^./™ coincides with that of Gk and then - since roots of 
Gk form orbits of order k - the k roots should merge by groups of m into / = k/m roots of Gfc/m, and there are 
exactly / such groups. In other words, whenever a group of m roots of Gk merges with a root of Gk/m, so do 
I = k/m other groups, and R{Gk, Gk/m) is an l-ih power of an irreducible quantity, named r{Gk, Gk/m) is (|22|l . 

In a little more detail, if ai, . . . , with k = Im are roots of Gk and . . . , are those of Gk/m, then 
ai — f3i implies, say, that 

ai = .- 

C^m+l = - 



Then in the vicinity of this point 

R{Gk,Gi)^Y[{a,-(3s) 

i,s 

(factors, which do not vanish, are omitted). Internal products are polynomials of the coefficients of / with 
first-oidei zeroes (the roots themselves are not polynomial in these coefficients!), and only external product 
enters our calculus and provides power I in the first line of ea. l|^ . 



• = "m = Pi , 
- - = Cl2m = 02, 

= • ■ • aim = Pl- 



l I sra 



n n 

s=l \i={s-l)m+l 
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Similarly, 



n 



n =ni n ("-/^^ 



■ .,3 = (B-l)r, 



J 



what implies that zero of D{Gk) is of order l{m — 1) = fc — to. 

The remaining constituent d{Gk) of discriminant D{Gk) describes intersections among order-fc orbits. In 
this case, whenever two points of two orbits coincide, so do - pairwise - the k — \ other points. Thus the 
corresponding zero is of order fc, and irreducible quantity d{Gk) enters DiGk) in fc-th power. 

For examination of examples in ssElandlHlit is convenient to have explicit versions of (|22|l for a few lowest 
n and k. It is also convenient to use condensed notation: 



dk^d{Gk), ru ^ r{Gk,Gi), 
Tki = 1 unless I is divisor of fc, l\k (or vice versa, fc is divisor of I). Then 



D{Fk) 



D{Gk) 



di 
didlrl^ 
didl 



3' 31 



«l"2"4'4l'42'21 



5' 51 



did2 d'^dgrli r^^ J^f i ''31 



di 
dlr2i 

'^3'''31 



"4' 41 '42 

dirt, 

"6 '61 '62 '63 



Decomposition of discriminant imply that, as / is varied, new roots can emerge in different ways, when 
different components of I?* are crossed. If the new roots are born when the component of i?(G'„, Gk) is crossed, 
where the orbits of orders n and fc intersect, then they appear at positions of the previously existing roots of 
Gk- If this happens at real line, than the phenomenon is known as period doubling (it is doubling^ since when 
more than two new roots occur at the place of one, they are necessarily complex). 

4.4 Period-doubling and beyond 

The simplest example of orbit reshuffling with the change of the map / is the period-doubling bifurcation |13| , 
which can be described as follows. Let xq be an invariant point of /, i.e. Fi{xq) = f{xo) — xq ~ 0. Let us see 
what happens to an infinitesimally close point xq + e. 

f{xo + e) = xo + f'{xo)e + 
r^{x„+e)^xo + [r{xo)fe + ..., 

We see that the necessary condition for xq + e with infinitesimally small but non-vanishing e to be a root of 
Fi or F2 is f'{xo) = 1 or f'{xo) = ±1 respectively (i.e. Fi or F2 should be degenerate at the point xq). The 
period-doubling bifurcation corresponds to the case f'{xo) — —1, i.e. the map / is such that ^2(2;), but not 
-^1(2;), becomes degenerate, then in the vicinity of the corresponding stable point (the common zero xq of F2 
and F2) a new orbit of order 2 can emerge. If, more generally, xq is a root of some other F„, Fn^xo) — 0, 
i.e. describes some (perhaps, reducible, if n is not a simple number) /-orbit of order (period) n, then the same 
reasoning can be repeated for /°" instead of /: 

Fn{xo + e) = F^{xo)e + . . . 
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and 

F2n{xo + e) ^ F^^{xo)e + . . . 

Derivative ^2)1(2^0) = Fn{^o){Fn{xo) + 2). Period-doubling bifurcation occurs for / with the property that 
(r")'(xo) = -1 or F:,ixo) = -2. 

The period-doubhng bifurcation, though very important, is not the only one possible: new orbits can emerge 
in other ways as well. 

First of all, doubling is relevant only in the case of maps with real coefficients. In fully complex situation one 
can encounter the situations when a higher power /'(xq)'"' = 1, k > 2 (and all lower powers of f'{xo) 7^ 1): then 
we have the bifurcation when the period increases by a factor of k and the new orbit emerges in the vicinity of 
original one (which loses stability, but survives). 

Second, the new orbits can emerge "sporadically" at "empty places", with no relation to the previously 
existing orbits and no obvious criterium to warn about their appearance. The only reason for them to occur is 
the crossing between discriminant variety V and moduli spaces Mon, Mn of iterated maps or their irreducible 
constituents. 

4.5 Stability and Mandelbrot set 

Above consideration of period-doubling bifurcation implies introduction of the notion of "stable orbits" in the 
following way: the point x of the orbit of order n, i.e. satisfying Gn{x; /) = 0, is called "stable" if 

i(r")'(^)i < 1 

and "unstable" otherwise. Since 

n 

(r")'(2^)=n/'(/°'(^))= n /'w (23) 

k—l all z G orbit 

all points of the orbit are simultaneously either stable or unstable. 

The Julia set J(/) C X is attraction domain of stable periodic orbits in X. Its boundary dJ{f) consists of 
all unstable periodic orbits and their grand orbits. 

Excluding x from the pair of stability conditions 

r G„(a:;/)=0, 

\ \F;,ix-j) + 1\ < I, 

we obtain a stability domain 5„ of the order-n orbits in moduli space A4 . All zeroes of reduced discriminants 
d{Gn) and resultants r(G„, Gmn) with arbitrary to > 1 lie at the boundary of Sn- This property is used in our 
description of Mandelbrot set in sOl 

Mandelbrot set describes exchanges of stability between pairs of orbits: at the boundary of Mandelbrot set 
stable orbits intersect with unstable ones, stable become unstable, while unstable become stable. Since stable 
orbits lie inside Julia set, and unstable ones - on its boundary, this causes reshuffling between interior and the 
boundary of J(/). Since all this actually happens with entire grand orbits, reshuffling involves infinitely many 
points and looks like a fractalization of the boundary. In fact, this does not exhaust all possible bifurcations of 
Julia set J(/): they can be also caused by crossings and degenerations of unstable orbits (with no reference to 
the stable ones). In order to study bifurcations of Julia sets one should include into consideration the pre-images 
On,s of periodic orbits of /, associated with the roots of the functions 

Fr^A^) : = r("+^)(a;) - r{x) = F„ (/-(x)) (24) 

and study their reducibility properties. Examples in secsEl and El below demonstrate that consideration of 
orbit's pre-images, discriminants D(Fn_s) and resultants R{Fn^s, Fk^r) can indeed capture the bifurcations of 
J(/), which are overlooked by consideration of orbits alone. Presumably, it provides complete theory of Julia 
sets in terms Grand Mandelbrot set, which characterizes reshufflings of all orbits, stable and unstable, and is 
fully algebraic, does not refer to additional stability structure. 
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4.6 Towards the theory of JuUa sets 

4.6.1 Grand orbits and algebreiic Julia sets 

For a given (algebraically closed) field X each series / = X]n=o '^n^;", a„ G X, defines a map X ^ X and thus 
a (pre-)order f{x) ;^ a; on X. Then the set of points of X splits into connected components with respect to >- 
("grand orbits"): we say that xi,X2 belong to the same grand orbit, iff xi y x,X2 y x for some x G X., i.e. 
/°"(.Ti) = /°™(.T2) for some n and m. In particular, for n = 0, m = 1 the set of points X2 satisfying xi = /(.T2) 
is the /-preimage of a given xi. So a grand orbit for generic point x G X can be represented by an oriented tree 
with the valency of each given vertex x equal to d+ 1, where d is the number of roots of the equation f{z) = x. 
For a given x the set a;+ of points x' y x will be called the orbit of x. 

However, if x G X is a root of some Fn{x; /), then the orbit becomes a closed loop of finite length (order) 
n, and the grand orbit GO is a + 1-valent graph (see Fig.), obtained by gluing vertices at a distance n on 
a certain totally ordered chain of the general position tree. For a pair of elements Xa,xp on a periodic orbit 
of order n we have 0„ = a^a =2;^, moreover both Xa >- and x^ >- Xa is true, so for the points of 0„ the 
relation >- is just a pre-order. But for each x G On there is a subset x~ C GO, called the pre-orbit of x, being 
a well-ordered tree rooted at x. The set of points z G for which f°^{z) = x will be denoted (not to be 
confused with a negative power!), so that x~ = UgX"*. Each set belongs to the roots of 

FnA^) = Fn+s{x) - F,{X) = r^^+'\x) - f^x). 

Then for each n we have a finite set of grand orbits with n-periodic "bases" 0„. Then we may characterize 
the initial map / by the structure of this discrete set of data. 

In particular, taking z G x~ we can regard z~ C x~ as a sequence of maps we can take its inverse limit, i.e. 
the set of sequences {zi} where zq = z,Zi = f{zi+i) (or Zi G x~*). For different z',z" (which may in general 
belong to different grand orbits as well) their preorbits {z')~, {z")~ are not only isomorphic as ordered sets, but 
have isomorphic inverse limits. In particular, the set of limit points for lim(z')~ and liin(z")^ in F coincide. 

We call them the algebraic Julia set of /: 

J^(/) = lim(x)- VxgX. 

Hypothesis: JaU) coincides with J(/). 

4.6.2 Prom algebraic to ordinary Julia set 

In order to explain our expectation that algebraic Julia sets are related to conventional ones, we formulate two 
more hypotheses. 

A map / defines in X two important subsets: the unions 0+(/) and 0-(/) of all stable and unstable 
periodic orbits of all orders - both are countable sets of points in X. 

Hypothesis: For almost all / G the /-orbit of almost any point x G X approaches 0+{f). Moreover, it 
approaches exactly one concrete stable orbit of particular order n(x), which is a characteristic of the point x. 

The pre-orbit tree of almost any point x G X has infinitely many branches. Going backwards along particular 
branch we approach its inverse limit, or origin. 

Hypothesis: For almost all / G and x G X the origins of almost all periodic branches of grand /-orbit 

of X belong to 0-{f). Moreover, different periodic branches usually originate at different unstable orbits from 
0-{f) and almost every orbit in 0-{f) is an origin of some branches of /-pre-orbit of x. 

In other words, almost each grand orbit originates at the closure of entire 0-{f) and terminates at (tends 
to) a particular orbit in £)+(/). If these hypotheses are true, one and the same set C-(/) is an origin of almost 
all grand orbits (not obligatory bounded) . It is this set (or its closure, to be precise) that we call the boundary 
of the algebraic Julia set, 9Ja(/) = C_(/). Periodic branches and unstable periodic orbits play the same role 
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as periodic sequences (rational numbers) in the space of all sequences (real numbers), see s l5.5l for some more 
details. 

Bounded grand orbits reach (not just tend to) a periodic orbit, however, in this case it can belong to ©_(/), 
not obligatory to ©+(/): bounded grand orbits are not in generic position in what concerns the "future". Still, 
they are not distinguished from the point of view of the "past" and can be used to study Ja(/). Bounded grand 
orbits are convenient to deal with, because they consist of roots of Fn^s (with all s) and can be studied in pure 
algebraic terms. 

4.6.3 Bifurcations of Julia set 

If we deform / (change the coefficients of the corresponding series) then the set of grand orbits move in F and 
may undergo the following two structural changes: 

(i) merging of two distinct periodic orbits O', O", which happens as soon as any pair of elements from these 
orbits merge, 

(ii) merging of elements of the same pre-orbit . 

This merging of components of grand orbits results into splitting of JaIJ) into disjoint components. 

Case (i) corresponds to merging of the roots of F„ and ii) corresponds to merging of the roots of F„ s- 
This means that / hits the discriminants Von = 'D{Fn) and 'D{Fn^s) with s > 1 respectively, i.e. that some 
functions Fn{x; f) or Fn.s{x; f) become degenerate. If a; is a multiple root of Fn then Fn{x) = and F'^{x) = 
(/°"(x) — x)' — so (/°"(x))' = 1. Since n-periodic a; is a stationary point for the map /°", then the last 
equation means that the points of the merging orbits and O'^ change their stability type. Thus in the 
case when we can speak about convergence in X, the language of "stability" is related (equivalent?) to that 
of "discriminants", but we will stick to the latter since it allows to handle also the case ii) and the case of 
the arbitrary field X as well (which does not need to be full and have continuous and differentiable functions 
defined). 

In the case (ii) some point of the grand orbit has degenerate preimage (two or more preimages coincide). A 
point S X has degenerate /-preimage when discriminant 

D{f{x) ~ zf) = 0. (25) 

This equation defines an important (multi- valued) map 

z: M — » X, 

associating a set of points {zj} - solutions to (|25(l - with every map f E A4. Bifurcations of the type (ii) occur 
whenever some bounded grand orbit crosses this set, i.e. when 

K,.(^/;/) = o 

for some n and k. 

The fact, that all bifurcations of bounded grand orbits are either of the type (i) or of the type (ii), implies 
that discriminants D(Fn^s) sue products of discriminants D{Fn) and the new canonical functions Fn^r{z) on M, 
with r < s. Similarly, the resultants R{Fn^s, Fk^r) are made from R{Fn, Fk) and various Fm,t{z) with m\n,k 
and t < s,r. Actual expressions are somewhat more involved, because the functions Fn.s{z) are highly reducible 
(see s 14. 71 below 1 and their different components enter differently into formulas for particular resultants. 

For description of these reductions it is important that the map z f is intimately related to the critical points 
Wf of /. Indeed, ea. (|25|l implies that f{x) — Zf and f'{x) have a common zero. The zeroes of derivative f'{x) 
are critical points Wf: 

/'K) = o, 

and therefore 

i.e. there is a one-to-one correspondence between the points Zf and Wf. Moreover, from the definition of Fn^s 
it follows that 

F„^s{zf) = Fn,s+iiwf). (26) 
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4.7 On discriminant analysis for grand orbits 
4.7.1 Decomposition formula for Fn.s{x; f) 

Each zero of Fk^rix) satisfies f°^''~^^\x) = f°^{x). Applying f°'' to both sides of this equation, we obtain 
jo(2fc+r)^^) ^ jo(/c+r)(^-) ^ Repeating this procedure several times we get /° (2;) = /"'■(x) and, 

finally, applying f°('^-^) we obtain (x) — f°'^{x) for any s>r and any m. This means that Fk^rix) is a 

divisor of any i^n,s, provided k\n and s > r. Consequently, similarly to H2Q(I we have the following decomposition 
of Fn,s(x) into irreducible (for generic field X and map /) components: 

r(n) s 

FnA^\f)=\{\{GkA^;f) (27) 

k\n r=0 

The irreducible function Gkfl{x; f) = Gk{x; f) appeared already in H2()|l . In particular, H27|l states that 

-^1,1 = G1G14, 
Fi^2 = GiGisGi^2, 

Fi^s = GiGi^iGi^2 ■ ■ ■ Gi^s, 



F2,s — {GlGi^iGi^2 ■ ■ ■ Gi^s){G2G2,lG2,2 ■ ■ ■ G2,s), 

F^^s — (GiGi^iGi^2 ■ • ■ Gi_s)(G3G3_iG3_2 • • ■ G^^s), 

F4,s — (GiGi_iGi_2 ■ • ■ Gi_s)(G2G2,lG2,2 ■ • ■ G2,s)(G4G4_iG4^2 ■ • ■ G4_s), 



4.7.2 Irreducible constituents of discriminants and resultants 

A direct analogue of H21I) expresses discriminants D{Fn.s) through DiGk.r) and the resultants R{Gk,r,Gk' ,r')- 
A less trivial thing is expression of these quantities through elementary constituents, which are, as already 
predicted, the familiar from ea. (|22|l irreducible discriminants and resultants dk — d{Gk) and rnk — r{Gn, Gk), 
as well as the new quantities, which are irreducible components Wk.rif) of 

WnAf) ^Y[GnAwfJ) (28) 

Wf 

- the products of all values of Gr.s{x; f) at all critical points Wf of /, f'{wf) — 0.^^ Note, that G, not F, enters 
the definition of W in (jJHl, still after substitution of peculiar values of x = Wf this quantity often gets further 
reducible, the first few reductions are: 



Wnif) = wAf) (irreducible), 

WiAf)^Mf)wiAf), 
( W2,2r{.f) ^ W2^2r{f) (irrcduciblc) , 

\ VF2,2r+l(/) = W2(/)w2,2r+l(/), 

W3,3r(/) = W3,3r(/) (irreducible), 

W3,3,.+l(/) = W3(/)w3_3r+i(/), 

W3^3r+2if ) = W3,3r+2(/) (irrcduciblc) , 

^^We hope that the use of the same letter w for critical points Wf £ X and irreducible components «>fc^r(/) G will not cause 
too much confusion. 
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W4,4r{f) = W4Ar{f) (irrcducible) , 

W4,4r+l(/) = W4(/)w4,4r+l(/), 

W4,4r+2(/) =W4,4r+2(/) (irrcducible) , 
Wi^r+sif) ^ W4Ar+3{f) (irreducible), 

(we actually checked most of these statements only for n + s < 6). Presumably, in general 

Wn,s — Wn^s for s 7^ 1 mod n (irreducible). 



(29) 



4.7.3 Discriminant analysis at the level (n, s) = (1,1): basic example 

Let us begin the proof of ()29fl from the simplest case of Wi,i- The key point is that 

Gi,i{x) = fix) mod Gi{x). (30) 

Then, according to H55() . the residue R{Gi,i, Gi) can be expressed as a product over all critical points Wf, i.e. 
the roots of f'{x): 

ri|i,i = i?(Gi, Gi,i) ^ i?(Gi, /') ^ n^iH = = (31) 

w 

Ea. H31|l is the first result of grand orbit discriminant calculus. It relies upon the relation (|30|l . actually, on the 
fact that f'{x) is residual of Gi^ix) division by Gi{x): then, though Gi^i{x) is not divisible by Gi{x) as a 
function of x, Wi{f) = JJ^, Gi{w) divides VFi,i(/) = Hto Gi^iw) as a functional of /. In order to see that 

Gi^ijx)- f'jx) ^ /(^)-^ - / [X) 

Gi{x) f{x)-x ^ ' 

is non-singular at all zeroes of Gi{x) = Fi{x) = f{x) — x it is enough to consider an infinitesimal variation of 
such root, X ~ xi + X: fi^i) — ^ii arid expand all functions at xi up to the second order in x'- 

G,(x) = Fi(x) = f{x) - X = (/'(xi) - l)x + \f"{x4)x^ + ...= 



(/'(.r)-l)x(l + X^(^gi^ + 0(x^)) 



F^A^) = /'(.Ti)(/'(xi) - l)x + \nx,){nx4f + /'(xi) - l)x' + O(x^) = 



, , f"{x,){nx4? + r{x,)-i)x , ^ r(^i)x . 2 

2(m)-l) ^""^^ V (v'"2(/'(.r)^l)+^(^ ■ 



m) + ^x/"(a:i)(/'(:ri) + l) + 0(x'), 



so that H32|) becomes 



(f(.r)-l)x + 0(x-) ^ ^ ^"^^ + ""^^^ 
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and is finite at x = 0. This proves divisibility of Wi.i(/) and allows to introduce its irreducible constituent 

W,{f)^l[G,iwf)=w,{f), 

Wf 

Wf 

It enters expression for discriminant £)(Gi.i): 

D{Gi^,) ^ 4-^wi,i (33) 

i.e. is actually the irreducible part of this discriminant: 

di^i'^wi^i. (34) 

Power d — 1 in which di — d{Gi) enters ()33|l depends on the degree d of the map f{x) (i.e. f{x) is assumed to 
be a polynomial of degree d). 
Eqs.(El| and 

ri\i,i ^ wi, 
dis ~ Wis 

together with 

i?(Gi,i,G'„) 1, forn>l 

are the outcome of discriminant/resultant analysis in the sector of Gi and Gi^i, responsible for fixed points and 
their first pre-images. 

4.7.4 Sector (n, s) = (l,s) 

For a fixed point xi = f{xi) we denote /' = f'{xi), /" = f" {xi). Then for x = xi + x '^6 have: 

Fn{x) = x(/'" "1)1 + ^xf'^jr^ + O(x') 



and 



Gk{x) = gkif) [l + \xf"hk{n + 0(x')) 



where QkiP) are irreducible circular polynomials (which will appear again in sEJ and hki/3) are more sophisti- 
cated (with the single exception of hi they are also polynomials). The first several polynomials are: 

<?!(/?)= /3-1 higiip)^! 

ff2(/3)=/3 + l h2{f3) = l 

53(/3) +1 ft3(/3)=/3+l 

54(/3)=/3' + l /i4(/3) =/3(/3 + l) 

55(/3) =/9' + /3' + /3'+/3 + l /i5(/3) =/3' + + 1 

56(/3) =/3'-/3 + l /i6(/3) = /3* + /3' + /3'-l 

Making use of Fn^s — F71+S — Fg and of decomposition formula H27|l . it is straightforward to deduce: 

Gi.A^) = - /' + ix/"(/' + 1) + O(x^), 



Gi,i(x)-f(x) 1 /' + 1 - 2 , ^f'^nt \ 
G^ix) ^2-f3^^ +0(X)=2/ +0(X), 
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as we already know. Here and below f'{x) = /' + xf" + 0{x^). Further, recursively, 

= ^ = r + \xnr + i)r"' + o{x% 

i.e. Gi^s(x) — f'{x) mod (71(2;) and, as generalization of H3U|I and (|31|l . we obtain for all s 

ri|i,, = i?(Gi, Gi,.) ~ i?(Gi, /') ^ n /) = 

Wf 

4.7.5 Sector (n, s) = (2,s) 

For consideration of the sector (n, s) = (2, s) we need to consider the vicinity of another stable point X2, which 
belongs to the orbit of order two, i.e. satisfies f{f{x2)) = X2, but X2 = f{x2) 7^ X2- Denote /' — /'(X2), 
/' = /' (i2) = nf{x2)), f" = nx2), /" = /"(i2) = /"(/(X2)). For X = X2+ X we have: 

F^{x) = G,{x) = {X2 - X2) + xif - 1) + \x^f" + O(x'), 



(a;) = G^{x)G2{x) = x (/'/' - l) + {f"f + ff') + O(x'), 
3(2;) = Gi(x)G3(x) = (i2 - X2) + X (r/' - 1) + (/"/'/' + /"/'' + /"(/'/')') + O(X'), 



Then 



G2 i(a;) — ^2,1 _ Fi^F^ — Fi) _ 
G1G2G14 (F2 — Fi)F2 

^ {i2~x2)+x{f-l) + h'f" + Oix') 

(i2 - X2) - xif - 1)/' - (/"(/' - 1) + /"/'') + Oix') 

X {f'f - 1) f + \x^ (/"./'' + /" ((/'/')' + ./'/' - 1)) + O(x^) 
X {f'f' - 1) + (/"/' + /"/'') + 0{x^) 

' Xif'f'-l) , ^.,2^^ , 1 f , 1^ , 



so that 



- + 4 _^^ + o(x') j (^/' + 2X/" (/'/' + 1) + O(x^) ) , 

G2i(x) + fix) = -xif'f - 1) (h" + + oix': 

\2 X2-X2J 
G2.ix) + f'ix)_ ^l;. + _^Uo(x) 



and 

G21 

G2{x) \2'' ' X2 — X2 

is finite at x = 0. Therefore 

?-2|24 - i?(G2,G24) - ~R{G2,f') -Y[G2{Wf) = W2{f) = W2{f). 

Wf 
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4.7.6 Summary 

In the same way one can consider other quantities, complete the derivation of (|29|l and deduce the formulas for 
resultants and discriminants. Like in H3,S|I . these decomposition formulas depend explicitly on degree d of the 
map f{x). In what follows 

#,s= (d-l)d^-i 

denotes the number of s-level pre-images of a point on a periodic orbit, which do not belong to the orbit. Also, 
Wn{f) enter all formulas multiplied by peculiar d-dependent factors: 



= d'i'('i'-l)u-4. 



22 = 4, 


33 = ^ 


22 = 4, 


36 = 7 


2^ = 64, 


324, 


12 = 4096, 


3^2 


230, 


3240 


2^4, 


37176 



(the last two columns in this table contain the values of numerical factors for d — 2 and d — 3, question marks 
label the cases which were not verified by explicit MAPLE simulations) . Obviously, these factors are made from 
degree Nd{n) of the map G„(a;), which was considered in s l4.1l 

WniD^d^'-^^^^Wnif). (35) 

Non-trivial (i.e. not identically unit) resultants are: 



urr' n ^ _ / if m = n, 



if m ^ n 



and 



n(n n ^ - / ^{^k, G„) = r*^ if s' = s, k <n (actually, k\n) 

[ W*'' if S < s' 



or, in the form of a table: 
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11 

12 
13 
14 

15 
16 



11 



1 
1 
1 

1 



12 



■iiii 
1 
1 
1 
1 



,r,#2 
1 

,r,#2 
1 

#2 
1 



w 



13 



Wl 

1 
1 
1 

1 



"^1 

-#2 



,r.#3 



14 



1 
1 
1 

1 



"'I 



-#4 
-#4 



21 



1 
Ul2 
1 
1 
1 



22 



1 
W2 
1 
1 

1 



23 



1 

W2 
1 
1 

1 



24 



1 

W2 
1 
1 
1 



^12 


1 


1 


1 


„#1 


1 


'l2 


1 


1 


1 


1 


1 


'l2 


1 


1 


1 


1 


1 


^12 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 
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1 
1 

W3 
1 
1 



32 



1 

1 

1 

1 



1 

r#2 
1 
1 
1 
1 
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1 

1 

1 

1 



1 
1 

r#3 
1 
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1 
1 

1 

1 



1 
1 

1 

'^13 
1 



21 

22 
23 
24 
25 
26 



,#1 
' 12 

1 
1 
1 
1 
1 



1 


1 


1 


#2 
12 
1 

1 


1 

#3 

n2 
1 


1 
1 

"^12 


1 


1 


1 


1 


1 


1 



wf' 
wf' 





wf' 


wr 


1 


1 


1 


1 




~%2 

wf 


-#2 

w1t 


1 


1 


1 


1 


-#2 
wJT 




wf" 


1 


1 


1 


1 


-#2 

wj 


w*' 




1 


1 


1 


1 


-#2 
Wo 


w*^ 


-#4 
Ulo 


1 


1 


1 


1 


-%2 

w^ 




wf 


1 


1 


1 


1 



31 


,.#1 

^13 


1 


32 


1 


„#2 
^^13 


33 


1 


1 


34 


1 


1 


35 


1 


1 



1 
1 

r#3 
1 



1 

1 

1 
r#4 

1 





wr 




w*' 
w*' 




W^ 


-#2 






w!; 




wf' 


-*2 
wZ 


-#3 

wZ 


wr 


'#■2 
W^ 


wr 



wf 

-#3 



-#4 
"3 



41 


'l4 


1 


42 


1 


„#2 


43 


1 


1 


44 


1 


1 


45 


1 


1 



1 

,#3 

1 
1 



1 
1 
1 
„#3 



{rl,)*' 



{rid*' 
1 

1 
1 



{rid*' 



1 
1 

{rid*' 



51 


'^15 


1 


1 


52 


1 


„#2 

rib 


1 


53 


1 


1 


„#3 

rib 


54 


1 


1 


1 


55 


1 


1 


1 



1 
1 
1 

r#4 
1 



61 
62 
63 
64 
65 



,#1 

'^16 
1 

1 
1 
1 



1 
,#2 

1 
1 
1 



1 

1 
,#3 

1 
1 



1 

1 

1 
,#4 

1 



{rh)*' 
1 

1 
1 
1 



{rh)*' 
1 

1 
1 



{rid*' 



{rid*' 
1 

1 
1 
1 



{rid*' 



{rie)*' 



Similarly, for discriminants: 



r=2 
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<]1 



,n—k 
k.7i 



W. 



(5(n,s) 



n 



(36) 



k<n 



nr' s — 1 r—1 

Exponents in this expression are #s = [d — l)d*~^, 4j=s/jj=r = d^"'', and (5(n, s) = c?'^~^4f^^T" ~ — ^'^'^ 

s = nr' + r, r', r > (if such r' and r do not exist, i.e. s < n, 6(n, s) = 0). The first few values of 5{n, s) are 
listed in the table (stars substitute too lengthy expressions): 



/ 



n \ s 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 





1 


d+1 


d'^ + d + l 


+ d2 + d + 1 


* 


* 




* 


* 


2 








1 


d 


rf2 + 1 


d^ + d 


d^+d'^ + 1 


d^ +d^ +d 




* 


3 











1 


d 


d^ 


d^ + 1 


d'^ + d 


d^ +d^ 


d^ + d^ + 1 


4 














1 


d 






d^ + l 


d^ +d 


5 

















1 


d 




d^ 


d* 


6 




















1 


d 




d3 


7 























1 


d 




8 


























1 


d 


9 





























1 


10 

































The first few discriminants are listed in the table (£)„ = D(Gn), not to be confused with irreducible c?„; the 
last two columns contain numerical factors for d = 2 and d = 3^^): 

^^Numerical factors are equal to d'i'C'"'— where the sequences #„(d) are 1, 1, 3, 6, 15, . . . for d = 2 and 1,2,... for d = 3. 
One can observe that #„{d) = deg^Gn{x; x'^ + c), but the reason for this, as well as the very origin of numerical factors, here and 
in eg. 1351 . remain obscure. 
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d -. 



11 

12 

13 
14 
15 



Df '^''-''>wf;wf.,Wi3 



D 



,d(d-l) 
*1 



Mill 



.d^(d-l),„d+l 



■'I 



UI 



rd''(d-l)^d'^+d^+d+l^ 



Il«'f2"'l3"'l4"'15 



1 

2" 
2l6 
248 
2128 



-3 
327 

3I35 

q567 



21 
22 
23 
24 
25 

31 
32 
33 
34 
35 



D. 



d(d-l) 



14^21^22 



d^(d-l),,.d 



d: 

.d^fd-i) 



^21 M^22 ^^231^24 

wiwiwiwi^w26 



d: 



d(d-l) 



1^11^32 



^-S "31 
-jd3(d-l)^,,rf3, 

^d^d-1) d4^d3 Wd= Wd W35 



D," ^" ^'Wg'^i W'3i;W3i5W34 



,2(d-l) d-l 

^12 '^21 
,2d(d-l) d(d-l) w 
<^2 ''12 "'21 '"22 

j2d^(d-l) d^rd-l) rf2 , 

"2 ''12 '"''2'(i'2l«'22'^'23 

,2d^(d-l) d^(d-l) d w-' d 



,2d«(d-l) d*(d-l)^ d^ + l 



'12 '"2 "^21 

,3(d-l) 2(d-l) 

,3d(d-l) 2d(d-l) rf 
3 13 t«3j^ty 

,3d2(d-l)^2d^(d-l),„d2,„d 



Ulo-, «)22 '"'23 "^24 '"'ZS 



13 



,3d3(d-l) 2d3(d-l) 



3 



,3d* (d- 1) ^2d* (d- 1) ,„d,„d\„d' 



"'31'"' 32 '"'33 
Ul3Ul|j MiIj ^"33 "'34 



'13 



«'3"'3l'^32«'33«'34"'3S 



1 

2" 
2I6 

248 

2128 

1 
2I2 

248 

2144 

2384 



36 
354 

? 

? 

? 

324 
? 
? 
? 
? 



41 
42 

43 
44 
45 



D 



Df'-''-^^wiwi.,W43 

d^(d-l)„.d3„,d2Tr^d ■ 



,d*(d-i) 



'vyfi w^^j wi'3W44 
^i? w£ wt, W45 



,4(d-l) 3(d-l) 2(d-l) 
"4 '14 ' 24 

,4d(d-l) 3d{d-l) 2d(d-l) 



1041 

W41 «'42 



,4d^(d-l) 3d^(d-l) 2d^(d-l) d^ rf 

,4d3 (d- 1) 3d3 (d- 1) 2d3 (d- 1) d3 d2 d 
«4 ^^14 ^24 M'4if42 "'43^44 

^d* (d- 1) 3d* (d- 1) 2dnd- 1) d4 d3 .2 ^ 

"4 "14 "24 '''4f'4iW42 ^43^^44^46 



1 

224 

296 
? 
? 



51 
52 



Dg-Vsi 

nd(d-i)^dw52 



,5(d-l) 4(d-l) 
'^5 ''is ■^Sl 

,5d(d-l) 4dtd-l)^„d 



"'51W62 



61 
62 



d(d-l) 



W^, W62 



,(i(d-l) 5(d-l) 4(d-l) 3(d-l) 

"^e '"16 ''26 ''36 '"'61 

Yd-1) 5d(d-l) 4d(d-l) 3d(d-l) d 

r^g rjg r3g tU2i^22 



' 26 
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4.7.7 On interpretation of Wn,k 

Responsible for degenerations of pre-orbits are intersections with the set {2;/} = f{{wf}), which is /-image of 
the set of critical points of /. Let us describe what happens when {zf} is crossed by pre-orbits of different 
levels. 

•{Gn(/)=O}n{^/}^0. 

Let A be the intersection point, denote its pre-images along the periodic orbit through A ^,A . . ,A " = 
A, and the pre- images of level s on the pre- image tree through B~'^ ^^^{A), where indices i can be further 
ordered according to the tree structure. Points on the level ,s of the tree, rooted at pre-image will be 

denoted through B^^'' If ^ G then two of its pre-images coincide. There are two possibilities: 

either these coinciding preimages belong to pre-orbit tree, say B^^{A) = B^^iA), or one of them lie on the 
orbit, say B^^{A) = A~^. In the first case Wn,i{f) = 0, in the second case Wn{f) = 0, and the fact that only 
two such possibilities exist is reflected in decomposition formula 

XlGn{Zf) = []G„,i(w/) = TK,i(/) = 

Zf Wf 

(products at the l.h.s. are needed to build up quantities, which depend on the coefficients of /(x), not on their 
irrational combinations, entering expressions for individual points Zf and Wf - the latter are roots of f'{x)). 
Let us analyze these two cases in more detail. 
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Fi gurc 9: A periodic orbit of order n and grand-orbit trees, rooted at the points of the orbit. Wn{f) vanishes when (i) a point A 
on the orbit coincides with an "irreducible point" from the set {zf}, which is /-image of the set of critical points {wf}, and (ii) 
when one of the coinciding first pre-images of A is also a point on the orbit. Arrows show the pre-image points of different levels, 
which coincide when A S {zf} and Wn{f) = 0, i.e. when A~^ = B. 



o 



w„{f) = 0. See FigEI 

Since B^^{A) = A^^ exactly one point from {G„.i — 0} coincides with exactly one point on {Gn = 0}, thus 
the resultant has simple zero, 

Numeric d-dependent factor, distinguishing between Wn and w„ requires separate explanation. 

Pre-images of B^^{A) - the set of d (different!) points B^^{A), . . . , B'^^^A) from {G„.2 = 0} should coincide 
pairwise with pre-images of A"^, which are d — 1 points B^^{A~^), . . . , B'^\{A~^) from {Gn,i = 0} and the 
point A~^ on the orbit, i.e. from {G„ = 0}. Thus the corresponding resultants will have zeroes of orders 
^1 = d — 1 and one respectively: 

i?(G„,i,G„,2)=rf,|„2^^i;f 

and 

R{Gn, Gn,2) = rn\n,2 ^ Wn 

Next pre-images of A^^ consist of A^^ e {G„ = 0}, #1 = d - 1 points B^^{A^^) e {G„,i = 0} and 
#2 = d{d - 1) points B-^{A-^) £ {G„,2 = 0}, while those of Bf^ are #3/#l = d^ points from {B-^{A)}. 
Each of d^ points from the A~^ second pre-image should coincide with one of the d^ points from that of B^^. 
This provides relations 

i?(G„,2,G„,3)=rf,, 3~u;f , 



R{Gn,i,Gn,3)^rt'-^t 



n,2|n,3 
' »i,l|n,3 



and 

R{Gn, Gn,3) = rn\n,3 ~ Wn 

Continuing along the same line, we deduce: 

RiGn,r, Gn,s) = rt'^r\n,s ~ for r < S 
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and 



W„ 



Actually, this exhausts the set of non-trivial (non-unit) resultants for orbits of coincident periods n, but 
does not exhaust the possible appearances of ?i;„; it will show up again in discriminants, see below. 

o wn,i{f) = 0. SeeFigCni 

In this case we deal with a single tree, rooted at A G Gn n {zf} and thus we can safely omit reference to 
A in B{A). Thus, the starting point is ^ = -62^^. (In particular, there is nothing to discuss in the case of 
d — 2, when ^1 = 1, there are no a priori different points at the level s — 1, which could occasionally coincide, 
and all Wn^i = 1.) We can conclude, that discriminant D{Gn,i) ~ Wn,i- 

This is not the full description of discriminant, because it also contains contributions coming from intersec- 
tions of different periodic orbits, which cause the corresponding points on pre-orbit trees to coincide as well. 
These factors were already evaluated in s 14. 31 above, in application to discriminant of Gn,s they should just be 
raised to the power #s, counting the number of pairs of points which are forced to coincide at level s when 
their roots on the orbits merge. Thus finally 



D{Gn,l) ^ D{Gn)*^WnS = 



< n ^fc,7' 
\ k<n 



and the remaining undetermined factor is just a constant on M, which depends on d, but not on /. The same 
argument explains the expressions for non-trivial resultants 

R{Gn.s, Gk^s) = i?(G„, Gk)*' = r'^f, k\n, k < n. 

Considering next pre-images of 1 and we obtain at level s the two coincident sets of = d** ^ 

points each and thus D{Gn,s) ^ wf^ i D{Gn)'^^ ■ However, many more factors enter in the expression for higher 
discriminants: all Wn.r with r < s contribute, and, remarkably, Wn also contributes when s > n - in accordance 
with decomposition formula H29(l . To see all this we should return to the very beginning and consider the 
intersections of pre- orbits with the set {z/}. 

• {G„,,_i(/) - 0} n{2/} 0, r > 1. See FigCH 

Let '"^ e {-2/}- This means that some two preimages of this point coincide, say B^^ — -B^*^, and at level 
s >r there will be two coinciding sets, of #s/#r = d^~^ points each. Since, Gn,r-i{zf; /) = Gn,r-i{f{wf); /) ~ 
Gn,r{wf; f), all this happens when Wn,r{f) — 0, and we conclude that D{Gn,s) ~ T^nr fo'" ^'^Y f < s. We 
can further use (|29|l to substitute W„,r(/) by w„^rif) for r 7^ 1 mod n and by Wn,rif)wn{f) for r = 1 mod n. 
Collecting everything together we finally reproduce (|36|1 : 



n.r 



I 



A;<n 



V 



Index S{n, s) is obtained by summation over all 1 < r < s, such that r ~ 1 mod n, of the weights 4j=s/^r = 
Each term in this sum corresponds to the tree, warped times around the periodic orbit, see Fig ll2l For 
s = nr' + r and r' , r > we have 

d(n,s] = d — — = ; ; 

If such r' and r do not exist, i.e. if s < rt, 6{n, s) — 0. 
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Figure 10: The same periodic orbit of order n and the same point A, intersecting {zf} as in the previous Fig|5| «'ri.i(/) vanishes 
when (i) a point A coincides with a point from {zf}, and (ii) when both coinciding first pre- images of A belong to pre-image tree, 
rooted at A. Arrows show the pre-image points of different levels, which coincide when A £ {zf} and ««n,l(/) = 0. 
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Fi gure 11: A point which intersects {zf} = /({'if/} belongs to the r — 1-th level of the grand orbit tree, r > 1. Arrows show the 
pre-image points of different levels, which coincide in this situation. 




4.8 Combinatorics of discriminants and resultants 

Systematic analysis of orbits of the map / : X ^ X includes determination of the following characteristics, 
used in our presentation in s|21 

• Number iV„(/)/n of periodic orbits of order n, which is n times smaller than the number Nn{f) of roots 
of Gn{x;f), 

7V„(/)=degJG„]. 

• Number S{n; /) of different combinations of periodic orbits of order n which can be stable for some /. (It 
deserves emphasizing, that some orbits are never stable: for f = + c only one of the two fixed points - orbits 
of order one ~ has non- vanishing stability domain in the plane of complex c.) 

• Number A/'„(/) of elementary components a in stability domain 5„. It is also useful to count separately 
the numbers Afn'\f) of elementary components u at level p, 

oo 

• Number Nn{f ) of solutions to the equation G„(x = 0, /) = 0, considered as an equation for /. Solutions 
of this equation for particular families of maps fj, G M are also important, because they define some points 
inside the cr-domains (and thus - approximately - the location of domains) . If the family /i is one-parametric 
and parameter is c, then 

AA„(/)=deg,[G„]. 

Presumably, 5(n; /) = Af^if) = K.{f). 

• Number of zeroes of reduced discriminant dn = d{Gn), introduced in s l4.3l For one-parametric families of 
maps / the number of zeroes is equal to degj,[d„]. Zeroes of dn are associated with self-intersections of periodic 
orbits and define the cusps on various components. Usually the number q of cusps depends only on the level of 
p of the component ai^^ (and not on n and other parameters like and a^). 

• Number of zeroes of reduced resultants Vn^k = f{Gn,Gk)- For one-parametric families of maps / the 
number of zeroes is equal to deg^[r„_/j]. Each zero describes some merging of the elementary cr-domains. The 
number of zeroes is related to the number of a-parameters. 

All these numbers are related by sum rules: 

oo 

Y^Af^^ =AegAGn{w.,c)], 

oo 

Y^qipWi^^ =deg,[dr.], 

OQ 

^AA(/)^5]degJr„,,]. (37) 

p=l k\n 

Example. For / which is a polynomial of degree d, f{x) = x'^ + c, we have: 

• n = 1 : 

• of orbits: d, 

• of "potentially-stable" orbits: 5(1; a;"* + c) = 1. 

Since Gi(0; /) = /(O) = ao = c, the number of solutions to Gi(0; c) = is 1. 

Discriminant d{Gi) — D{Gi) — D{Fi) vanishes when simultaneously x'^ + c ~ x and dx'^~^ = 1: there are 
d — 1 solutions, thus discriminant has d — 1 zeroes. 

• n = 2 : 

• of orbits: 
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# of "potentially-stable" orbits: S{2; x"^ + c) = d - 1. 

Since G2{0; c) = c''"^ + 1 = 0, the number of solutions to 6*2(0; c) = is 1. 

Discriminant d{G2) = yi^^^^jy has [d— l)(<i— 2) zeroes. Each of the d—1 elementary stability domains 
for orbits of order two has d — 2 cusps. 

We list much more data of this kind for particular families of maps in Tables in s|Hl below. Similar analysis 
can be performed for bounded grand orbits. 

4.9 Shapes of Julia and Mandelbrot sets 

4.9.1 Generalities 

The shape of Mandelbrot can be investigated with the help of stability constraints H10|l : 

|^^,;(x;/) + l| = l, 

G„(a;;/)=0. ^^^^ 

Exclusion of x from these equations provides a real-codimension-one subspace in C tM, which defines the 
boundary dSn of stability domain S'„. Its connected components are boundaries of the elementary domains a. 
Enumerating all n we can obtain in this way the entire boundary of Mandelbrot space M(fi). 

Zeroes of all discriminants D{Gn) belong to this boundary. Indeed, if D{Gn) = there is a point a;„ G X 
which is common root of Gn and GJ^: G„(a;„) — G'^{xn) = (actually, there is entire set of such points, labeled by 
additional a-index). Then, since F„ = ]\k\n = G„F„, we have F^{xn) = Gn{xn)F'n{xn) + G'^{xn)Fn{xn) = 
and thus both equations H38|l are satisfied. Furthermore, since according to (|22|l . a zero of every resultant 
i?(G„, Gn/m) is also a zero of D{Gn), we conclude that zeroes of all such resultants also belong to dM. In the 
world of iterated maps only resultants of this kind are non-trivial, i?(G„, Gk) = const unless k\n or n\k. 

One possibility to define the shape of M{ia) is just to plot many enough zeroes of discriminants and 
resultants and - since they are dense in dM{ii) - this provides approximation with any desired accuracy. 
Advantage of this approach is that it is pure algebraic and - once formulated - does not contain any reference 
to stability equations and to the notion of stability at all. The disadvantage is that - at least in presented form - 
it does not separate points from different elementary domains: zeroes of each particular resultant are distributed 
between many elementary domains, which can even belong to different disconnected components M^a of M . 
Also, it does not describe directly the boundaries of elementary domains, which - outside a few cusps, located at 
zeroes of reduced discriminants d{G„)~ are smooth curves of peculiar (multi-cusp) cardioid-like form. Finally, 
it does not explain the similarity of various components Mka , which - for reasonably chosen families fi - belong 
to a set of universality classes labeled by Zd-i-symmetries. 

Another option is to use equations H38(l more intensively and transform them into less transcendental form, 
making reasonable reparametrizations and approximations. The choice of parameterization, however, can im- 
pose restrictions on the families /i, but instead one can move much further in explicit description of constituents 
of Mandelbrot and even Julia sets. 

In the case of Julia sets the situation with non- algebraic approaches is somewhat more difficult. These sets 
do not have elementary smooth constituents, except for exactly at the bifurcation points, i.e. at the boundary 
points of Mandelbrot set, but even on dM the decomposition of J(/) depends discontinuously on the point at 
the boundary. As soon as one goes inside Mandelbrot set, infinitely many different smooth structures interfere 
and only some traces of them can be seen approximately. The structure of Julia sets is pure algebraic, not 
smooth (or, if one prefers, consistent with d-adic rather than with C-topology). The boundary of Julia set for 
given / is formed by solutions of all the equations Fn^s = for all n and s, with exclusion of a few orbits, which 
are stable for this /, - and this remains the best existing constructive definition of Julia set in general situation. 
Also, approximate methods can be used to analyze some features of Julia sets. 

4.9.2 Exact statements about 1-parametric families of polynomials of power-d 

Assuming that fj, d Vd and dimc/i = 1; we can label the maps / in the family by a single parameter c: fi = {/c}- 
For given c G Mi(^) define: 

F'^{x-c) + l^e'^'^^-^\ (39) 
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This relations maps the roots of Gn{x; c) (the periodic orbits of order n) into the unit circle, parameterized by 
the angle ip. Considering all possible n and taking the closure, one can extend this map to entire boundary of 
Julia set. For c G Mi this provides a one-to-one correspondence between dJc and the unit circle. If c e Mkoi, 
analogous expression, 

relates a part of dJc and unit circle. Other parts are mapped onto additional circles. Parameterization (|39() 
is adjusted to describe the boundary of Mandelbrot set: it explicitly solves the first equation in (|38|l and 
the function x(ip) can be substituted into G„(a;;c) = to obtain c{ip). The map c{ip) has several branches, 
associated with different elementary domains a from stability domain S'„. This program can be realized in 
certain approximation. 

Before going to approximations, let us give an example of exact statement. The map —> da is singular 
(there are cusps on the boundary da) whenever 



{d-l 
Gn = 



^' (40) 



Here (?„ = dGn/dc and 

Hn = {F^, Gn} = F^G'„ - F^Gn- 
Solutions to this system are all zeroes of the resultant 

R{Gn,G'J ~ D{Gn) = d"(G„) [] ^""'(Gn, Gk) (41) 

k\n 

which are not simultaneously zeroes of another resultant, R{Gn, Hn), or, what is much simpler to check, 
R{Gn, Hn). Actually excluded are roots of all reduced resultants r(G„,Gfe), and the cusps of da are at ze- 
roes of reduced discriminants d{Gn)- Indeed, whenever r(G„, G^) — for fc|n, there is a common zero x of the 
two functions G„(5;) = Gk{x) = 0. Also, because of (|^ . G'n{x) = 0, therefore Hn{x) = —Fl'Gn{x). Both G„ 
and Gk enter the product H2U|) for Fn, Fn ^ G„Gfc and two derivatives in F^' are not enough to eliminate G„, 
Gk and G'n from the product, therefore also F^'{x) — and Hn{x) — 0. 

For polynomial maps elementary domains of Mandelbrot set M(c) belong to universality classes, which are 
represented by multi-cusp cardioids (cycloids), described by the function e(<p) in complex e-plain 

Go: £ = e^'^, 

Cd-i : e^e"^- -e"^'^, for d>l. 
d 

The real curve Cd-i has discrete symmetry 

Zd-i : (fi^ (fi + e, £ ^ e'^e, 

it is singular, \de/d(p\ = 0, provided = fc = 0, 1, . . . , d — 2, i.e. the number of cusps is d — 1. 

Generically, the boundaries of elementary domains belong to the class Gi for level p = and to Gq for all 
other levels p > 1. However, for special families, like fc = x'^ + c, the situation will be different: at level p = 
the elementary domains cr^"-' can belong to the class Gd~i (has d — I cusps), while at higher levels it decreases 
to Cd-2 (has d — 2 cusps). If symmetry is broken by adding lower powers of x with small coefficients to fc, the 
small-size components M^a fall into generic Go © Gi class, and the bigger the symmetry-breaking coefficients, 
the bigger are components that switch from Cd-2 ® Cd-i to Go ® Gi class. 



4.9.3 Small-size approximation 

In order to explain how cardioids arise in description of elementary domains, we continue with the example of 
one-parameter family of polynomials of degree d and describe an approximation, which can be used to explain 
the above mentioned results. It is based on expansion of equations (|38|l around the point (x, c) — {w, c„), where 
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Wc is a critical point of fc{x) and c„ is the "center" of an elementary domain from S'„, defined from solution of 
the equations 

fci^c) = 0, , . 

G„K;c) = 0. ^^^> 

We now substitute into H38|l x — w + x ^^'^ c = c„ + e, expand all the functions in powers of e and Xi leave 
the first relevant approximation. Typically |x| < o?~" and |e| < actually there is much stronger damping 
for elementary domains belonging to " remote" M^a 7 thus approximation can be numerically very good in most 
situations. 

One can even promote this approach to alternative description of Mandelbrot set as a blow up of the system 
of points 142|l : each point gets surrounded by its own domain a. Some of these domains get big enough to 
touch each other - and form the connected components Mkoi - some remain disconnected topologically, but get 
instead connected by trails. Consistency of such approach and the one, based on resultant's zeroes, requires 
that the number of solutions to 142|l with given n coincides with the total number of zeroes of r„„i (c) with all 
m < n plus the number of components Mna growing from the elementary domain fT*-°^ [na], 

degj,G'„(wc, c) = l^-^"^ I + ^ degj,r„^m(c) V 1 — parametric famihes (43) 

where ^F{c) — deg^F(c) denotes the number of roots of F{c), which for polynomial F{c) is equal to its degree 
(for multiparametric families this is the relation between the degrees of algebraic varieties) . A similar sum rule 
related the number of roots of reduced determinant (i„(c) and the numbers Aflf'^ of solutions to H42I) with q{p) 
cusps: 

E'?(p)-^i"^=deg,d„(c). (44) 



4.9.4 Comments on the case of fc{x) ~ x'^ + c 

In this case there is a single critical point Wf = 0. 

Stability domain 5*1 for orbits of order one consists of the single elementary domain, 5*1 = a[l] and is defined 
by a system 

d\x^-'\ = 1, 

x'^-x + c = 0. (45) 
Its boundary can be parameterized as follows: x — li^^/'^'^^^^e*"^, and 

dS,: c = x-x' = ^^[e^^^^^e^^'^^ (46) 

This is a curve with a cusp at c = c? — l/d''/'^''"^' and the symmetry group Zd-i, the shift 0^0+ 
multiplies c by (so that there are actually d — I cusps). Julia sets have another symmetry, Zd. 

We use this example to show how approximate similarity of different components Mka of Mandelbrot set 
can be explained by "perturbation theory" from the previous subsection 14. 9. 31 We discuss just a very particular 
part of the story: the central domains cr(''^[fcQ;], obtained by the procedure from s l3.2.5l The procedure implies 
that we solve the equations ifTnji 

Fn{x] c) ~ 0, (A7\ 

|F^(x;c) + l| = l ^^'> 

by expansion near the point (x, c) = (0, co), where co is a root of the equation Fk{Q^ co) = 0. This means that 
we now write x = Q + x^ c = cq + £ and assume that x ^^'^ £ a-i'^ small - and this will be justified a posteriori. 
In this approximation 

F„(x, c) = c„(c) ~x + bnX"" + 0(x"'), 
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and 



K+1 = Kc't^d. 

At c = Co the Cfe(co) — and Cfc(co + e) = CkS (dot denotes derivative with respect to c). Substituting 
Fk — Cfc(co)£ — X + bk{co)x''' + 0(x^'^) into (|T7jl . we get in this approximation: 

d&fclx'-'l = 1, 
bkX'^ -X + CkS = 0, 

i.e. X = Pke'*, Pk = M-^'^'^-^) and 

a^, : £ = ^ ^ ^ f e'^ - 36'^^") (48) 

Cfc Cfe V " / 

Thus such 9(7^: are approximately similar to = dui , described by eq. 1461) . 

Approximation works because the radius pk is numerically small. Presumably it can be used to find Feigen- 
baum indices and their d-dependence. The same approximation can be used in the study of Julia sets for 
c S Mka with fc > 1. From this calculation it is also clear that the shape of Mka is dictated by the symmetry of 
the problem: for families p, of maps, which do not have any such symmetries one should expect all the higher 
Mka to be described by H48|l with d = 2. Obviously, Feigenbaum parameters depend only on effective value 
of d. Also, one can easily find families /x, for which above approximated scheme will not work: then one can 
expect the breakdown of self-similarity of the Mandelbrot set. 

However, if one tries to be more accurate, things get more sophisticated. Precise system (|38|l contains the 
condition G„(x;c) — rather than i^„(a;;c) = in (|47|l . In variance with i^„(x), which has a gap between x 
and x'', Gn{x) contains all powers of x, for example, for n — 2 



G2 = i + J2 c'^^'^^'x^ + J2^d -k- ly 



-k-2 d+k 
X 



k=0 k=0 



This makes analysis more subtle. We leave this issue, together with approximate description of Julia sets to the 
future work. 



4.10 Analytic case 

Before ending discussion of approaches to creation of theory of the universal Grand Mandelbrot (discriminant) 
variety and Julia sets and proceeding to examination of various examples, in the remaining two subsections we 
explain how the relevant notions look in the case of generic analytic functions, not obligatory polynomial. As 
already mentioned, the theory of discriminant, resultant and Mandelbrot varieties can be made pure topological 
and free of any algebraic structure, however in this case one can loose other interesting - pure algebraic - 
examples (like X — Fp, to begin with) as well as the skeleton machinery, based on the use of functions Fn 
and Gk (to define them one needs subtraction operation and entire ring structure respectively) and ordinary 
discriminant calculus (which is again algebraic). Thus in the remaining subsections we accept that these 
structures are important and explain instead, why we do not think that restrictions to polynomials is needed 
anywhere. Actually, polynomials play distinguished role in particular examples: this is because the vocabulary 
for iterations of other functions - even trigonometric ~ was never developed and we do not have any adequate 
language to discuss, say 

, . , sin(a;sina;a;) — a; 

G2(a;;sma;a;) = 

smujx — X 

or 

„ , . X sin(tj sinwa;) — sinwa; ^ , . 

Gi i(x; smwx) = = G2(x; smwx) — 1 

SVlYLJX — X 
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or there no-less-interesting trigonometric counterparts 

~ , . sinfwsinwcc) — sinx 

G2(a;; sinojx) = 

sin Lox — sin x 



- . . , sm(ujsinu)x) — smwx ^ , . n -, 

Gi i{X]smLOx) — : : — G2(x; smwxj — 1 

sm ojx — sm x 

While iterated polynomial is always a polynomial, just of another degree, even iterated exponential or iterated 
sine do not have names. This could be not a too important restriction if one agrees to rely more on computer 
experiments than on theoretical considerations, but even here existing symbolic-calculus programs are better 
adjusted to work with polynomials... Still, non-polynomial singularities play increasingly important role in 
modern theories, they appear already in the simplest examples of r-functions of KdV and KP equations, nothing 
to say about more general partition functions. Thus we find it important to emphasize that discussion of phase 
transitions and their hidden algebraic structure in the present paper is by no means restricted to polynomials. 

Take an arbitrary (locally) analytic function f{x). By the same argument as in the polynomial case / maps 
roots of Fn{x) into roots, giving us the a representation of Zn on the set of zeros of Fn{x),^'^ which is the 
analogue of the Galois group action for polynomials. The orders of orbits of this actions are also divisors of n, 
but the set of those orbits now is infinite for each given n. Although for series there is no notion of divisibility, 
still if n is divisible by fc, then Fk{x) ~ implies Fn{x) = 0. Indeed, if n = km and /°'^(xo) — xq, then 
/°"(a;o) = /°('="')(a;o) = r^^^'''~^Hr^{xo)) = f^'^^^'-^^^xo) = . . . = xq, i.e. every root of Fk{x) is obligatory 
the root of Fn{x). 

Example: Take f{x) — sin(2;) + cx + hx'^. FigllHIshows real roots of ^4(2;) and F2{x) for c = 2, h — —0.01. 
The root x w 2.39 is the root of F2{x) and its orbit has order 2. The root x k, 2.03 is the root of G4 (independent 
of F2{x)) and its orbit has order 4. 

If U is the definition domain (the set of regular points) of /, then f o f and F2 are defined in [/02 = Un f{U) 
and so on. Denote by 

f/°°(/) = n^=or"(f/) 

the inverse limit of the sequence 

U ^Uo2 ^ ■■■ ^ Uon ^ ■■■ 

For polynomial f{x) the set U°°{f) = C is the whole complex plane. 

The Z„ representation on the roots of Fn is also a representation of Z, which we denote hy V„. In n is 
divisible by k, then every root of Fk is the root of Fn and we have an embedding, which is actually a morphism 

-"^^This is because if xq is a root of Fn{x) = f°^{x) — x, then f{xo), f°^{xo) etc are also roots, and f°"{xo) = xq. 
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of Z-representations 14 — > y„ . This provides a directed set of Z-representations 

/ \ 

\ / 
/ 

Its direct limit KxjC/) = U„Ki is a representation of Z which is a pertinent characteristic of the map /. All 
roots of every F„ satisfy x = = /"^"'"^(x) for any m, thus they belong to ?7°°, i.e. each orbit of Kxj 

totally belongs to \J°° . Thus we have the following hierarchy of subsets, characterizing a holomorphic function 
/: 

Ko(/) c C . . . C (7o2(/) c C/(/) c C 

The closure of V"oo-space is called the algebraic Julia set J a of /. 

Hypothesis: This definition coincides with the one in sec 14. 61 

Example: For /(x) — x"^ the Juha set J(/) is the unit circle S*^ C C and the Vao space is an everywhere 
dense subset in this circle. 



In the general case, when / is a series with complex coefficients, there is no distinguished /-invariant 
subsets analogous to R C C. But the existence of /-orbits implies, that a deformation of coefficients of /(x), 
which makes any two roots of some Fn{x; f) lying on different orbits in T4o(/) C U°°{f) coincide, results in 
simultaneous pairwise merging of all the roots on those orbits, so the resulting orbit acquires multiplicity. 

For given n take the subset C„ in the space of coefficients a of f where some of the orbits in the set of zeros 
of Fn merge, i.e. / belongs to the /„-pullback of the discriminant component Von- The closure of Coo := U„C„ 
is the boundary dM of the Mandelbrot set of /. 

Hypothesis: For any finite-parameter family /i C of maps Coo(/i) is a subset of real codimension 1 in 
the family /i). 

Provided this is true, dM separates the space of complex coefficients into disjoint components. Note, that 
every particular C„ has real codimension 2 (since the two equations Fn{x) ~ and F^(x) — can be used to 
define a complex common root x and one complex relation between the coefficients of Fn (x) - and thus between 
the coefficients of /). 

4.11 Discriminant variety I) 

This section reminds the standard definitions of discriminants and resultants for the case of polynomial functions 
and briefly comments on the possible ways of generalization to the case of arbitrary analytic functions. 

4.11.1 Discriminants of polynomials 

In the polynomial case ^{Vn) is defined ^3] as an algebraic variety in Vn by an equation 

D{F) = (49) 
in the space of coefficients. Here D{F) is the square of the Van-der-Monde product, for F{x) = f „ Hfe^i {^~^k) 

n 

D{F)=vl--^l[{ak-air (50) 

k<l 
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Remarkably, it can be also expressed as a resultant of F{x) and its derivative F'{x) 

D{F)^{~r^-'y\-'R{F{x),F'{x)), 



(51) 



which is a polynomial of the coefficients of F[x) (while the individual roots can be sophisticated functions 
of these coefficients, not even expressible in radicals), so that the equation (|49|l is indeed algebraic. 

The resultant of two polynomials F{x) = Yll^^a'^kx'^ = ''^nnfc=i(^ ~ "^fc) ^"^^ ^{^) = SfcLo^fc-^'^ = 
w„i YViLii^ ~ A) is defined as a double product 



R{F, G) = n(«^ - Pi) =<X{ G{ak) = (-)"™< n 



(52) 



k,l 



1 = 1 



depends on symmetric functions of the roots {a} and/or {/?}, thus - according to Vieta formula - depends only 
on the coefficients of F and/or G, and is representable as determinant 



R{F, G) = det 

{rn-\-n) X {m-\-ri) 



/ Wo Vi ... Vn 

Wo Vi ... Vn 

Vo Vi 
Wo Wi ... Wm 

Wo Wi ... Wm 



\ 



V 



Wo Wi 



(53) 



The resultant obviously vanishes when F(x) and G{x) have a common root xq, accordingly the matrix in H53|) 



annihilates the column- vector (l,a;o,Xo, 



n+m 



)• 



Directly from the definitions it follows that discriminant of a product of two polynomials is decomposed in 
the following way 



D{FG) = D{F)D{G)R^{F,G). 



(54) 



If F{x) — F[{x) is divisible by G(x), F — H mod G (for example, H is the residual from F division by G), 
then H52|l implies that 



R{F,G)^ \{ F{(i)^ \{ H{P)^ R{H,G) ^ [] G{j) 

P: G(/3)=0 13: G(/3)=0 ■y: H{-()=Q 



(55) 



where proportionality signs imply the neglect of factors like v„ , w„i and minus signs in (|52|l . in this form the 
relation is applicable to any analytic functions, not obligatory polynomials. (Note that the products at the 
r.h.s. of (|55|l over the roots of H and, say, of _ff + kG will be the same - up to above-mentioned rarely essential 
factors.) 



4.11.2 Discriminant variety in entire A4 

To define the discriminant in analytic (non-polynomial) case one can make use of both definitions (|50|l and H51|l 
and either take the limit of a double-infinite product over the roots or handle determinant of an infinite matrix. 
It is an open question, when exactly this can be done and when the two limits can coincide. 

A possible approach in the case of determinantal representation can make use of the following recursive 
procedure. Discriminant of a polynomial oq -I- aix -I- . . . -|- a^^ix'^^^ may be written as (pre-factor a'^^^ in H51|l 
is important here) 

d 

Dd+i{ao, . . . , Qd+i} = al ■ Dd{ao, . . . , a^} + ^ &d+i,fc{ao, . . . , a^} • a^+i, 

fc=i 

where Dd{ao, . . . , ad} is discriminant of oq -I- aix -|- . . . -I- adx'^ of the previous degree and bd+i,k{ao, ■ ■ ■ , id} are 
polynomials of coefficients ao, . . • , ad- Then for an analytic function locally represented by series f{x) = ^ a^x* 
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with numerically given coefficients {a,} we have a sequence 

Di = 1, D2 = of - 4aoa2, 
D3 = al{al - 4aoa2) + (18aoai02 - 40^)03 - 27ala\, 

Dd,... 

of numbers, equal to the values of the corresponding polynomials. The numeric limit of this sequence may be 
regarded as the value of the infinite expression 

...{al-...{al-{al-{al-l + 62i{ao, ai}a2) 

+hi{ao, ai, 02)03 + &32{ao, ai, a2}al) 
+64i{ao,ai, 02, 03)04 + 642{ao, 01,02,(13)04 + 643(00, oi, 02, 03)04) 

d 

+ •••)+ XI ^d+i.jioo, . . . , ad)ajj+i) . . . 

i=l 

(this computation is an algebraic counterpart of the chain fractions) which may be called the discriminant 
D{f) of the analytic function /. 

Since in the case of polynomial / the coefficients of all functions F„ are polynomials {al^\a)} of the 
coefficients {ok} of the scries /, such procedure can work simultaneously for / and Fn{f). The possibility of 
generalizations from polynomial to analytic case can be formulated in the form of the following 

Hypothesis: {o,) e C„ iff D{Fn) = 0, i.e. some orbits within the radius of convergence of series f = Y1 ^i^^ 
merge iff the above defined discriminant D{Fn) of series Fn = X^a^"^(o)x'' is zero. 

Example: For fc{x) = e^' - c we have Di = 1, D2 = 4(c- 1), D3 = 4(c- 1), D4 = 8(c- l)(2c- 1)^, . . .. 

It is easy to check that all higher Dk will be also proportional to c — 1, so that discriminant D , defined 

as numerical limit of this sequence vanishes at c = 1, where the two real roots of fc{x) merge (and go to complex 
domain) . 

4.12 Discussion 

The above considerations suggest the following point of view. 

Given a (complex) analytic function f{x) we get in the domain U°°{f) C X a set of discrete /-invariant 
indecomposable subsets. These subsets are orbits of a Z-action on U°°{f), generated by /. The union of 
periodic orbits gives a Z-representation Vodf) which may be regarded as an intrinsic characteristic of /. Each 
periodic orbit belongs to a set of roots of some G„(x; /) of minimal degree n, which can be assigned to the orbit 
as its degree. As the shape of / changes, i.e. as / moves in the moduli space M, the /-orbits move in X. The 
dynamics of this motion is worth studying. For /'s of certain shapes some orbits can merge. The corresponding 
subset in A4 can be regarded as inverse image (puUback) V* of the discriminant variety T> C Ai, induced by the 
functions Gn{x; /). We believe that discriminants and resultants can be well-defined not only for polynomials, 
but also for analytic functions, while the variety V* makes sense in even more general - topological - setting. 
We identify the closure of the union ^mn=i'^*nn = Uj^jl?* as the boundary of the Mandelbrot set, which 
separates the moduli space M into disjoint components, what can by used to classify analytic functions. Even 
more structures (Julia sets, secondary and Grand Mandelbrot sets) arise in a similar way from consideration of 
pre-orbits and grand orbits of / and of their change with the variation of /. 

Even in multi-dimensional case for every map f{x) € M. one can take its iterations Fn{x; /) = /°"(x) — x 
and construct the real-codimension-two subsets and 7?.J^„ in M., made out of zeroes of reduced discriminant 
dn{f) and resultant rmnif) functions.^'* A union Um=i ^nm = has real codimension one in M. and can be 
i^The function dn can be defined as irreducible constituent of 

D =1 ^"(^//) = 
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considered as a boundary of a codimension-zero stability set Sn C A^, so that d„ e dSn and dSnf^dSm = T^mn- 
What is non-trivial, the intersection of stabihty domains, not only their boundaries, is just the same (and has 
real codimension two!): dSn U dSm = T^tnu- This makes the structure of universal Mandelbrot set M = yJnSn 
non-trivial, and requires description in terms of projections to multipliers trees and representation theory, as 
suggested in this paper. Sections of the Mandelbrot set by real-dimension-two manifolds (by one-complex- 
parametric map families Hc <Z M) help to visualize some of the structure, but special care should be taken to 
separate intrinsic properties of M from peculiarities of particular section M(/i) = M U /i. For example, the tree 
structure is universal, while the number of vertices and links is not. 



and Tmn = whenever Fm {x; f) = Fnix; f) =0 (though seemingly overdefined, this system defines a real-codimension-two subspace 
in if m|n or n\m). Stability set Sn can be defined as 




where e.v.(A'.) denote eigenvalues of matrix A. 



51 



5 Map f{x) = + c: from standard example to general conclusions 



In this section we use the weU pubhcized example to illustrate our basic claims. Namely: 

1) The pattern of periodic-orbit bifurcations in real case is well described in the language of discriminants. 

2) Discriminants carry more information in two aspects: there are more bifurcations in real case than 
captured by the period-doubling analysis and there are even more bifurcations in complex domain which are 
not seen in real projection. 

3) Fractal Mandelbrot set is in fact a union of well-defined domains with boundaries, which (i) are described 
by algebraic equations, (ii) are smooth almost everywhere (outside zeroes of associated discriminants), (iii) are 
densely populated by zeroes of appropriate resultants. When resultant vanishes, a pair of smooth domains touch 
and the touching point serves as a single "bridge", connecting the two domains. 

4) Julia set is related to inverse limit of almost any grand orbit: its infinitely many branches originate in 
the vicinities of infinitely many unstable orbits, which constitute the boundary of the algebraic Julia set. 

5) Changes of stability (bifurcations of Julia set) always occur when something happens to periodic orbits, 
but these orbits that matter are not necessarily stable. Unstable orbits can be best studied with the help 
of grand orbits, because even pre-images of stable orbits carry information about unstable ones: all grand 
orbits originate in the vicinity of unstable orbits (grand orbit's periodic branches are more-or-less in one-to-one 
correspondence with unstable orbits). At the same time, everything what one can wish to know about grand 
orbits is again describable in terms of discriminants and resultants. 

6) If / S Ml, the fractal Julia set J(/) is in fact a continuous (but not necessarily smooth) deformation of a 
unit disc. The unit circle - the boundary of the disc - is densely covered by grand orbits, associated with unstable 
periodic orbits, and one stable periodic orbit, together with its grand orbit, lies inside the disc. The vanishing 
of the relevant resultant implies that the stable grand orbit from inside comes to the boundary, intersects 
and exchange stability with a particular unstable grand orbit, which becomes stable, quits the boundary and 
immerses into the Julia set. This results in merging the corresponding points of the boundary and changes the 
topology of Julia set - a bifurcation occurs. Topology of the Julia set J(/) is defined by position of the map 
/ in the Mandelbrot set and by the path through " bridges" , which connects / with the central domain of the 
Mandelbrot set (i.e. by the branch of the powerful tree Ti). 

7) If / e Mka with fc > 1 the disc/ball is substituted by fc-dependent collection of disjoint discs/balls, grand 
orbits are jumping between the discs, but return to original one every n times, if the orbit is of order n. All 
bifurcations when / moves between different components inside the given Mka are described in the same way as 
in 6). If / leaves Mandelbrot set, f ^ M, there arc no stable orbits left and Julia set looses its "body" (interior 
of the disc), only the boundary survives, which contains all the bounded grand orbits. 

5.1 Map f(x) = x'^ + c. Roots and orbits, real and complex 

fc{x) = -|- c is a well known example, examined in numerous papers and textbooks. This makes it convenient 
for illustration of general arguments. 

5.1.1 Orbits of order one (fixed points) 

• Si{fc) - is the set of roots of 

Fi(x;/,) = x2-x + c = Gi(x;/,). 
Since the order of polynomial fc is d = 2, the number of roots of Fi — Gi is two, each of the two roots, 

is individual orbit of order one. 

The orbits — roots are real, when c < j. 

The action of /(I) leaves each of the orbits = roots intact. The action of /(I) can be lifted to that of entire 
Abelian group Z, the same for all periodic orbits of arbitrary order (see s l4.1U|l . This action, however, does not 
distinguish between conjugate and self-conjugate orbits, to keep this information an extra Z2 group is needed. If 
we denote by Vq representation of Z when the group acts as cyclic permutations on the sequence of q elements. 



52 



and label the trivial and doublet representations of Z2 by superscripts and ±, then the representation of 
Z (g) Z2 associated with -^1(2;) is 2vi = ® . 

Discriminant d{Gi) — D{Gi) = D{Fi{fc)) = 1 — 4c, its only zero is at c = j, where the two order-one orbits 
intersect. At intersection point the orbits are real, x — ^ and remain real for real c < j. On the real line it 
looks like the two real roots are "born from nothing" a,t c — j. 

Since f'{x) — 2x the plus-orbit is unstable for all real c, the minus-orbit is stable on the segment —3/4 < 
c < 1/4. 

In the domain of complex c stability region 5*1 is bounded by the curve \2x\ = 1 where x is taken from (|56|l . 

i.e. 



. |l±Vl-4c| = 1, or 
^^1- c=l/4 + e^^sinV/2. ^ ^ 

In polar coordinates (r, (p) with the center at j the curve is given by r = sin^ ^ , see Fig ll4l Our expectation 
is that the points on this curve (an everywhere dense countable subset in it) will be zeroes of the resultants 
R{Fi, Fn), or r{Gi, Gn) to be precise, with all possible n. 

The critical point of the map fc = x'^ + c, i.e. solution to f^ix) = is Wc = 0. Equation Fi{wc', fc) — 
-F'i(0;/c) = /c(0) = has a single solution, c = 0, which lies inside the single elementary component of 
5i =a(o)[l]. 



5.1.2 Orbits of order two 

• S2{fc) - the roots of 

F2ix-Jc) = (x^ -x + c)ix^ +x + c+l) = Giix-Jc)G2ix;fc). (58) 
F2 is divisible by Fi since 1 is a divisor of 2, the ratio 

G2{x; fc) = F2/F1 =x^+x + c+l. 

The cP = A roots of F2 form two conjugate order-one orbits and one (^-^-^ = 1) self-conjugate order-two 
orbit. The new orbit - in addition to the two order-one orbits inherited from Si{fc) - is of order 2 and consists 
of two points — ^ ± ^~2~"^'^ ~ the zeroes of 6*2. 

The action of / (2) interchanges these two points. Representation of Z (g) Z2 , associated with i^2 , is2i;i©ti2 = 

Discriminant 

D{F2) = (4c- l)(4c + 3)^ = 
= DiGi)DiG2)R\G2, Gi) = d{Gi)<f{G2yiG2, Gi), 

and 

D{Gi) = d{Gi) = 1 - 4c, 
D{G2) = d^{G2)RiG2,Gi) = -3 - 4c, 

d(G2) = -1, 
i?(G2,Gi) = r(G2,Gi) = 3 + 4c. 

The only new zero is at c = — |. Since it is a zero of R{G2, Gi), the new orbit, associated with the roots of 
G2{x;fc), intersects an old one, related to zeroes Gi{x;fc)- Indeed, at c = —3/4 the stable order-one orbit 
intersects the order-2 orbit and loses stability, while the order-2 orbit becomes stable and stays real for all real 

Since at intersection point c = — | the two roots of G2 merge, D[G2) has a simple zero. Similarly, below 
all zeroes of R{G,G) will be also zeroes of either D{G) or D{G) (in particular, D{G) are usually reducible 
polynomials over any field X), moreover their multiplicities are dictated by the properties of intersecting orbits. 
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4. 



^1 'v. .1 



Fi gure 14: Stability domains in the plane of complex c: a) for orbits of order one and two, b) for orbits of orders up to four. 



Stability criterium for the order-2 orbit of fd^) is 

r \nx)nf{m = \^xUx)\ = \'ix{x^ +c)\<i 

1 G2(a;;/c) =a;2+a; + c+l = 0. 

The boundary of stability region S2 in complex plane c is obtained by changing inequality for equality. Then 
x{x'^ + c) = —x{x + 1) = c + 1, and the boundary is just a circle 

dS^: \c+l\ = \ (59) 

with center at c = — 1 and radius i, see Figll4l It intersects with the curve (|57|) at c = — |, which is the zero 
of i?(G2,Gi) = 3 + 4c. 

Eq. G2{wc\ fc) = G2(0;c) = c + 1 = now implies that c = —1, and this point lies inside the single 
elementary domain of 5*2 = (t'-^''[2|1]. 

5.1.3 Orbits of order three 

• Sz{fc) ~ the roots of 

^3(^; fc) ^{x^-x + c) (x^ + x^ + x\?>c+l)+ x^{2c + 1) + 

+ x^{ic^ + 3c + 1) + x{c^ + 2c + 1) + (c^ + 2c2 + c + 1)^ = 

= Gi(x;/,)G3(x;/e). 

^3(0;) is divisible by Fi{x) ~ Gi[x) since 1 is a divisor of 3, the ratio G3(a::) = F3[x)/ Fi[x). 

The <P = S roots of F3 form two conjugate order-one orbits and two (^-5-^ — 2) conjugate order-3 orbits. 

Representation of Z (g) Z2, associated with F3 is 2vi ® 2u3 = {v^ ® ) ® {v^ ^^3)- 

Discriminant 

D{F:i) = (4c- l)(4c + 7)3(16c2 + 4c+7f = 
= D{Gi)D{G^)R^{G^, Gi) = d{Gi)D\G^)r\G^, Gi) 

and 

D(Gi) =d{Gi) ^ l-4c, 
Z?(G3) = d3(Gi)r2(G3, Gi) = -(4c + lf{l&c^ + 4c + if, 

d(G3) =-(4c + 7), 
i?(G3, Gi) = r(G3, Gi) = 16c2 + 4c + 7. 



54 



The new - as compared to Fi (but not F2, which is not a divisor of ^3) - zeroes of D{F:i) are at c ~ — | and 

at c = — I ± The latter two values are zeroes of R{G3, Gi) and they describe intersection of two order-3 

orbits with the stable order-1 orbit (one of the two roots of Gi) at essentially complex c. No trace of these 

intersections is seen in the plane of real x and c. The two order-3 orbits intersect at the remaining zero of D(F3), 



which is real, c = — |, and they (orbits) remain real at all real c < 



-|. However, since c 



-| is actually a 



root of ^(Ga), it has nothing to do with Gi, and the new real order-3 orbits look - from the real point of view 
- "born from nothing", not at the position of any other "previously existing" orbits. In complex domain the 
orbits always exist, and just intersect (coincide) get real at c = — |. 
Stability criterium for an order-3 orbit of fc{x) is 



\f'{x)f{f{x))f'irHx))\ = |8x/e(a;)/f (x)| 
= |8a;(a;2 -|- c){x'^ + 2cx^ + c)\<l, 
G3(a;;c) =0. 

This system is hard to solve explicitly. Solving instead 

^3(0, c) = c^ + 2c^ + c+l = Q 



we obtain three points, lying inside the three elementary domains of 



S3 = a(°)[3]UaW 



This is the first example when we observe the appearance of a-parameters. 

5.1.4 Orbits of order four 
• Sn{fc) - the roots of 

^4(2;; /c) = {x^ -x + c){x^ + X + C+1) [x^'^ + 6ca;^° + x^ + 

+ a;^(15c^ + 3c) + 4ca;^ -|- x^{2Qc^ + 12c^ + 1) + x^{Qc^ + 2c) -|- 
+ a;^(15c^ -I- 18c^ + 3c^ + 4c) -|- a;^(4c^ + Ac^ + 1) + 
+ x^ (6c^ + 12c* -I- 6c^ -I- 5c^ -I- c) -I- a;(c* -I- 2c^ + c^ 2c) -|- 

{(^ -I- 3c^ -I- 3c* 3c^ 2c^ -h 1)) = 
= F^{x)Gi{x) = G^{x)G^{x)Gi{x). 

F4 is divisible by F\ and F^ since 1 and 2 are divisors of 4, the ratio G4 = Fj^jF^ is a polynomial of order 12. 

The d* ~ 16 roots of F4 form two conjugate order-one orbits (inherited from G\), one self-conjugate order-2 
orbit (from Gi) and three = 3) new order-4 orbits, two conjugate and one self-conjugate. 



Representation of Z (g) Z2, associated with F4, is 

Ivx © W2 © 3?;4 = iyX ® ) I 



1U2 



Discriminant 



D(F4) = (4c- 1) (4c + 3)3. 
•(4c + 5)®(16c2 - 8c + 5)''^(64c3 + 144c2 + 108c + 135)* = 
: I)(Gi)i?(G2)Z?(G4)i?2(G2,Gi)i?2(G4,Gi)i?2(G4,G2) = 
= rf(Gi)rf2(G2)rf*(G4)r3(G2,Gi)r^(G4,Gi)r6(G4,G2) 



where 



D(Gi) 
i3(G2) = 



4c, 
4c, 
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D(G4) = (4c + 5f{16c^ - 8c + 5)^(64:0^ + Uic^ + 108c + 135)'', 

R{G2,Gi) = 3 + 4c, 
RiGi,Gi) = (16c2 - 8c + 5), 
i?(G4,G2) = (4c + 5)2 

and 

d{Gi) l-4c, 
rf(G2) = -1, 
did) = 64c^ + 144c2 + 108c + 135, 
r(C?2,Gi) = 3 + 4c, 
r(G4,Gi) = 16c2 - 8c + 5, 
r(G4,G2) =4c + 5. 

The roots of D{F4) describe the following phenomena: 

c = J and c = — I were aheady examined in connection with Fi and F2. 

c = — I is a zero of D{Gi), thus it could describe either intersection of two order-4 orbits or merging of roots 
of a single order-4 orbit. In the latter case that orbit would degenerate into the one of order 2 or 1. Since this 
c is simultaneously the zero of R{Gi,G2), the right choice is intersection of self-conjugate order-4 orbit with 
the previously existing (inherited from G2) self-conjugate order-2 orbit. Since at intersection point two pairs 
of roots of G4 merge, the order of zero of D{Gi) is two. For lower real c < — | the self-conjugate order-4 orbit 
remains real. 

c = -J ± I are complex zeroes of R{Gi,Gi) and describe intersection of conjugate order-1 orbits with the 
conjugate ordcr-4 orbits in complex domain. No traces of this intersection are seen in the plane of real x and 
c. Since at intersection points all the four roots of G4 merge together, D{Gi) has cubic zeroes. 

The remaining three roots of D{F4) are roots of d{G4) alone and describe entirely the world of order-4 orbits: 
the G| 3 intersections of 3 orbits. Namely, 

c = — 1(1 -|- 2^/"^) - two conjugate order-4 orbits intersect, get real and stay real for lower real values of c. 
Look as "born from nothing" from the real perspective. Since 4 pairs of roots need to merge simultaneously, 
the order of zero of D{G4) is four. 

c = — I -h |2^/'^(1 ± i\/2>) - self-conjugate order-4 orbit intersects with one or another of the two conjugate 
order-4 orbits in complex domain. No traces at real plane. D{G4) has quadruple zero. 

Stability criterium for an order-4 orbit of fc{x) is 

\f{x)f{f{x))f{r^{x))fir^ix))\ = \2^xUx)ff{x)ff{x)\ = 
= \8x{x'^ + c){x'^ + 2ca;2 + c^ + c){x^ -|- . . .)| < 1, 
G4(a;;c) = 0. 

Solving instead of this system 

G4(0, c)=c^ + 3c^ + 3c^ + 3c^ + 2c^ + 1 = 
we obtain six points, lying inside the six elementary domains of 

1 |Ja(2)[2 2|1]. 

5.1.5 Orbits of order five 
• ^5{fc) - the roots of 

F,{x; /,) = {x^-x + c)(x30 + ...) = Gi(x)G5(x) 
There are = 6 order-5 orbits. 



54= u-^°^^"] U^^'^ 
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D(G5) =rf'(G5)r4(Gi,G5), 
i?(Gi, G5) = r(Gi, G5) = 256c'' + 64c3 + IGc^ - 36c + 31, 
d(G5) = -(4194304c" + 32505856ci" + 
+ 109051904c° + 223084544c* + 336658432c' + 
+492464768c^ + 379029504c^ + 299949056c^ + 211327744c^ + 
+120117312c2 + 62799428c + 28629151) = 

= - ((4c)" + 31(4c)i° + 416(4c)9 + . . .) (60) 

The powers 4 and 5 appear here because when an order-5 orbit intersects with the order-1 orbit, the 5 different 

roots of -F5 should merge simultaneously, thus discriminant has a zero of order 4 = 5 — 1, while at intersection 
of any two order-5 orbits, the 5 pairs of roots should coincide pairwise, so that discriminant has a zero of order 
5. 

G5(0; c) is polynomial of order 2 • 8 — 1 = 15 in c, thus equation G5(0; c) = has 15 solutions. These should 
describe the would be G| — 15 pairwise intersections of 6 order-5 orbits. However, since the power of c in d{G5) 
is only 11, wc conclude that only 11 of such intersections really occur. Instead there are 4 intersections between 
order-1 and order-5 orbits (since r(Gi, G5) is of order 4 in c). 

15= 11 4, 

and stability domain S5 consists of 15 elementary domains. 

5.1.6 Orbits of order six 
• ^^{fc) - the roots of 

Fe{x:f,) = Gi{x)G2ix)G3{x)Ge{x) = 
{x^-x + c){x^+x + c+l){x'' + ...) {x'^ + ...) 

Here the degrcc-54 polynomial Gg = F6F1/F2F3 = Fq/ Fi{F'2/ F\){Fi/ F\) since 6 has three divisors 1,2 and 3. 
There are ^'^ "'"'^ = 9 order-6 orbits. 

R{Gi,Ge) = r(Gi, Ge) = lec^ - 12c + 3, 
R{G2,Ge) = r2(G2, Ge) = (lec^ + 36c + 21)^, 
R{G3, Ge) = r^{G3, Gq) = (64c=^ + 12%c^ + 72c + 81)^ 
D{Gq) = -cf{Ge)RHGi,Ge)R^{G2,Ge)R{G3,Ge) = 
= rf6 (Ge , Gi ) (Ge , G2 )r3 (Gg , G3) , 

d{Ge) = 1099511627776c20 + 1044536G463872c^9 + 448738183086080^8 + 
-M21736553037824c^7 -t- 245929827368960c^6 -t- 399107688497152c^5 + 
-F535883874828288c^4 -|- 617743938224128c^3 -|- 631168647036928c^2 + 
-h576952972869632cH + 484537901514752c^0 -|- 376633058918400c^ + 
+263974525796352c^ + 1 735440 17002496c^ + 1049855221882880*^ + 
+589050857041920^^ + 33837528259584c'' + 155559159624960^ + 
+85587727468320^ + 1167105374568c + 3063651608241 = 

= (4c)20 + ... (61) 

G6(0, c) is polynomial of degree 32 — 1 — 1 — 3 = 27 in c: there are 27 solutions to G6(0, c) = 0. On the other 
hand there could be up to G| — 36 pairwise intersections between 9 order-6 orbits, however there are only 20, 
instead there are 2 intersections between order-6 and order-1, 2 - between order-6 and order-2, and 3 - between 
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order-6 and order-3 orbits (all these numbers are read from the powers of c in the reduced discriminant ^(G's) 
and reduced resultants ^(Ge, Gi)). We have 

27 = 20 + 2 + 2 + 3, 
stability domain Sq has 27 elementary components. 

In this subsection we used the example of fc{x) = + c to demonstrate the application of discriminant 
analysis. We explained in what sense the standard period-doubling bifurcation tree in Fig llSI is a part of a 
more general pattern. Even on the plane of real x and c much more is happening. Period doublings do not 
exhaust all possible bifurcations, other orbits of various orders are "born from nothing" - not at the positions 
of previously existing orbits - if we move from higher to lower c or " merge and disappear" if we move in the 
opposite direction. If we go away from real c and x, period doubling gets supplemented by tripling, quadrupling 
and so on. In complex domain no orbits can be born, merge or disappear - they can only intersect in different 
ways. Moreover, such intersections do not need to leave any traces on the real c ^ x plane. Specifics of the 
iterated maps Fn is that whenever they hit the discriminant variety the whole orbits intersect: many different 
roots coincide simultaneously. As n increases, the pattern becomes more and more sophisticated and begins 
looking chaotic as n — > oo. Still the whole story is no more no less than that of the intersecting algebraic 
varieties (infinite-dimensional, if we speak about infinitely-large n). We used the standard example in order to 
formulate the appropriate language for discussion of these - in fact, pure algebraic - phenomena. 

5.2 Mandelbrot set for the family fc{x) = + c 

Fig 1 141 shows stability domains for different orbits (roots of the polynomials Fn and G„) in the complex c plane. 
Fig 1 151 shows the full picture, where one can easily recognize the standard Mandelbrot set PHI- We see that it is 
nothing but the union of stability domains for orbits of different orders, which touch at single points: at zeroes 
of associated resultants. The picture also shows in what sense the boundary of every particular stability domain 
is formed by the zeroes of resultants (the zeroes constitute countable dense subsets in these boundaries). Note 
also that the two domains touch only if the order of one orbit divides another, otherwise the resultants do not 
have zeroes. For example, R{G2tG^) = 1, and it is easy to see that indeed these two r-polynomials can not 
have a common zero, since Gz{x\ fc) — {x^ + 2ca;^ + x ^ (? + c)Gi{x\ fc) = 1. Similarly, R{Gi, Gq) = 1 etc. 

According to this picture every point in the Mandelbrot set (i.e. every map / from a given family - not 
obligatory one-parametric) belongs to one particular stability domain, and there is a uniquely defined sequence 
of "bridges" (tree structure), which should be passed to reach the "central" stability domain - the one with 
a stable fixed point. As we shall see below, this sequence (path) is the data, which defines the structure of 
the Julia set J(/). As will be demonstrated in seclHI from the point of view of this picture there is nothing 
special in the maps fc{x) = x'^ + c, neither in quadraticity of the map, nor in unit dimension of the family: the 
tree-forest-trail structure of Mandelbrot set is always the same. 

However, particular numbers of vertices (elementary domains ti), links (zeroes of Tmn's), cusps (zeroes of 
dns) and trails are not universal - depend on the choice of the family /i of maps (on the section of the universal 
Mandelbrot set). Moreover, cyclic ordering of links in a vertex of a tree can also change, if fj, crosses a singularity 
of the resultant variety (i.e. if R{rki,rmn) = 0): the fat-graph structure of the forest is preserved only locally: 
under small variations of the 1-parametric families fi. 

5.3 Map f(x) = -|- c. Julia sets, stability and preorbits 

In this subsection we describe a few examples of orbits and pre-orbits to illustrate that normally the orbits tend 
to particular stable periodic orbits " in the future" and originate from entire set of unstable periodic orbits " in 
the past". We suggest to call the closure of this latter set 9J^(/) = ©_(/) the boundary of algebraic Julia set 
of /. We demonstrate that there is no difference between the behaviour of preorbits of different points a; G X: 
for X one can take invariant point (i.e. orbit of order one) of /, stable or unstable, a point of any periodic orbit 
or a point with unbounded grand orbit,- despite behaviour of the orbits is absolutely different in these cases. 
Often the boundary of Julia set is defined as pre-image of unstable invariant point, but our examples seem to 
demonstrate that this is unnecessary restriction. 



58 




Fi gurc 15: JMaiidclbrot set = union of all stability domains in the plane of complex c for the family of maps fc{x) = x + c. 
Boundary of the white region is densely filled by collection of all zeroes of all resultants = discriminant variety. The picture below 
shows the real section of the "Julia sheaf, by the plane of real x and real c: for each c on the real line shown are the real orbits of 
different orders, stable and unstable. 
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Now, for given c € M{ii) we can switch to what happens in the complex x plane (provided X = C). In this 
plane we have infinitely many periodic orbits. Their points form the set 0{fc) = 0+{fc) U O- (fc), the union of 
all zeroes of all Fn{x; fc)- There are a few stable periodic orbits and all the rest are unstable. For polynomial 
/c(x), all other stable orbits are actually located at x = oo, unstable orbits separate them from the finite stable 
orbits. Generic grand orbit approaches either one of the stable orbits from 0+{fc) or goes to infinity. The 
orbits which do not tend to infinity form Julia set J(/c) C X. The boundary of Julia set dJ{fc) has at least 
three possible interpretations: 

(i) it is a closure of the set of all unstable orbits, FigtTCl 



5J(/,) = 0_(/,); 

(ii) it is a closure of a bounded grand orbit, associated with any single unstable periodic orbit (for example, 
with an orbit of order one, FigEI^ unstable fixed point of fc, ~ or any other, Figll9land Fig l20|l . 



dJifc) = GO^; 

(iii) it is the set of end-points of all branches (not obligatory unstable, Fig l21l or periodic, Fig l22(l of any 
generic grand orbit (this set is already full, no closure is needed). 

Hypothetically all the three definitions describe the same object, moreover, for / e M, i.e. when the map 
belongs to the Mandelbrot set, this object is indeed a boundary of some domain - the Julia set. If / quits 
the Mandelbrot set, the object is still well defined, but is not a boundary of anything: the "body" of Julia set 
disappears, only the "boundary" survives. 

For the family fc = x^ + cwe know that there is exactly one stable orbit for all c G M(.t^ + c), see FigllSI 
for c ^ M{x'^ + c) all orbits (not just some of them) go to x = oo, so that J{x^ + c) — 9. 

• c € Ml. The set 0+ consists of a single stable invariant (fixed) point (orbit of order one) x^-* = i — ^^^^'^ . 

• c e M2- The set 0+ consists of a single stable orbit of order two formed by the points x'"^'' = — i± . 

We now elaborate a little more on the definition (iii) of the Julia set boundary. Generic grand orbit is a 
tree with d links entering and one link exiting every vertex. Bounded (non-generic) grand orbit of order n ends 
in a loop, with d — 1 trees attached to every of the n vertices in the loop. Degenerate grand orbits (bounded 
or unbounded) have some branches eliminated, see the next subsection 15.41 In the rest of this subsection we 
consider non-degenerate grand orbits. 

Going backwards along the tree we need to chose between d possible ways at every step. Particular path 
(branch) is therefore labeled by an infinite sequence of integers modulo d, i.e. by the sequence of elements of 
the ring Fd {ii d = p the ring is actually a field and the sequences are p-adic numbers). According to definition 
(iii) these sequences form the boundary of Julia set - this can be considered as a map between real line R 
(parameterized by the sequences) and OJa, provided by a given grand orbit. Different grand orbits provide 
different maps. These maps deserve investigation. Here we just mention that among all the sequences the 
countable set of periodic ones is well distinguished, and they presumably are associated with periodic unstable 
orbits, which are dense in OJa according to the definition (i). 

Important example: The map f{x) = (the c = case of -|- c) has a stable fixed point a; = and 
unstable orbits of order n consisting of points on unit circle, which solve equation x'^ — x and do not lie on 
lower-order orbits: 

e^p(^27:i^^y fc = 0,...,n-l. (62) 

The Julia set - the attraction domain of x = - is the interior of the unit circle (the unit disc) and all periodic 
unstable orbits fill densely its boundary, as requested by (i). Unstable fixed point is a; = 1 and its grand orbit 
consists of all the points exp (^^) , Z = 0, 1, . . . , 2™ — 1 of its m-th pre-images, which fill densely the unit circle, 
as requested by (ii). The grand orbit of any point x = re*"^, which does no lie on a unit circle, \x\ = r ^ 1, 
consists of the points x^ , with all integer k. For negative k — —m the choice of one among 2™ phase factors 
exp (^S^) should be made in order to specify one particular branch of the grand orbit. The branch can be also 
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Fi gure 16: The few lowest (of orders 1,2 and 3) periodic orbits. If all higher-order orbits are included, the closure provides the 
boundary dJ of Julia set for the family fc{x) = a;^ + c by definition (i). The same orbits are shown on different fibers of Julia 
sheaf: a) c = 0.2, b) c = —0.1 + 0.75i, c) c = —0.8, d) c = —1.76, e) c = 0.3. For a)-c) values of c e M there is exactly one stable 
orbit inside the Julia set and all others are at the boundary. 
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Fi gure 17; Grand orbit (circles) of unstable fixed point ~ ^ ~^ ' Pi'^ii^^S^s of orders from 1 to 4 are shown. Its closure 

provides the boundary dJ of Julia set for the family fc{x) = a::-^ + c by definition (ii). The same grand orbit is shown for different 
fibers of Julia sheaf: a) c = 0.2, b) c = —0.4, c) c = —0.8. The grand orbit of unstable fixed point consists of all the "extreme" 
points ("spikes") of Julia set for Re c < 0, while these are the "pits" for < Re c < 1/4. The pits are approaching the points of 
grand orbit of a stable fixed point, which lie inside the Julia set, and merge with them at c = 1/4, see Fig llSI Both "spikes" and 
"pits" are dense everywhere in the boundary dJ, if pre-images of all levels are included. 



62 




Figure 18: Decomposition of Julia set at c = 1/4, where the two fixed points, stable and unstable,- and thus their entire grand 
orbits - coincide, a) The pattern is shown as obtained by moving from inside the Mandelbrot set. White lines show peculiar 
circular /-flows in the vicinity of the stable fixed point (similar fiow patterns occur in the vicinities of all other cusps). When 
the boundary of Mandelbrot set is crossed, all the cusps (located at all the points of the coincident grand orbits of two merged 
fixed points, which densely populate dJ) dissociate into jjairs of points, thus dJ acquires holes (almost everywhere) and fails to 
be a boundary of anything: the "body" of Julia set "leaks away" through these holes and disappears, b) The resulting pattern 
as obtained by moving from outside the Mandelbrot set. Shown is unlargcd piece of the Julia set, which now consists ony of the 
everywhere discontinuous boundary. The holes grow with increase of Re c — i > 0. The spiral structure inside the former Julia 
set is formed by separatrices between the flows, leaking through the holes and remaining "inside". Infinitely many separatrices (at 
infinitely many cusps) are formed by (infinitely many) branches of the grand orbit of the fixed point. At the same time each point 
on the grand orbit is a center of its own spiral. Spiral shape is related to the nearly circular shape of the / flows in the vicinity of 
the stable fixed point and its pre-images just before the formation of holes (for Re c slightly smaller than i). If the boundary of 

Mandelbrot set is crossed at another point, which is a zero of some other dn (and belongs to some it£'''), the same role is played by 
grand orbits of the intersected orbits of order n: holes are formed at all points of these intersecting grand orbits. 
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— ^ it ^ ^ ^'^ . Again, the same grand orbit is shown for three different fibers of Julia sheaf: a) c = 0.2, b) c = —0.4, c) c = —0.8. 
This grand orbit also fills the boundary dJ densely, in accordance with the definition (ii). 




Figure 20: Non-degenerate bounded grand orbit of order 4 for cubic /(x), d = 3, c = 0.3 — one more illustration of the definition 
(ii). The first two preimages are shown. The dashed lines enclose each of the 9 second preimages of the whole orbit. 
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Fi gurG 21: Grand orbit of a stable fixed point = 0.2. Its end- (or, better, starting-) points approach the boundary dJ of 

Julia set, if one goes against the / flow and fill it densely, according to the definition (iii). The first four preimages are shown. 
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labeled by the soqiioncc of angles {v?,,, }. such that Lpm+i = ^(y'm + 27rsm), where {sm} is a sequence with entries 
or 1. If the sequence is periodic with the period ^ {^Pm + 27rp^-' , where 

Pm) = Xl"=o 2*STO+i- Then the origin of the corresponding branch is a periodic orbit of the order n, containing 
the point on the unit circle with the phase ip^tn = ^^z^, which solves the equation ipm^ = ^{(j)^^'> + 27r/)m^)- 
Different m. E (toq + 1, . . . , mo + n) describe n different elements of the periodic orbit, and mo labels the place 
where the sequence becomes periodic. For pure periodic sequences one can put mo = 0. 
For n = 1 we have two options: 

si = 0, pi = 0, (^(^) = and 

= 1, = 1, (^(1) = 27r (63) 

which have equivalent limits: the two branches originate at the single unstable orbit of order one, i.e. at the 
point X = 1. 

For n = 2 there are two options for the s-sequences: 

s = {0,l}, P2 = 2, 'P^'^ = y, 

s = {l,0}, P2 = l, = y, (64) 

which describe the single existing orbit of order 2 as the origin for two different branches. The sequences 
s = {0, 0} and s = {1, 1} are actually of period one, not two 
For n = 3 there are six options for the s-sequences: 

s = {0,0,1}, PS = 4, <p(3) = y, 

S = {0,1,0}, P3 = 2, <^(') = y, 
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S = {0,1,1}, P3 = 


6, 


^(3) = 


1 o 
IZTT 


s = {1,0,0}, P3 = 


= 1, 




Zi /I 

" T' 


S = {1,0,1}, P3 = 


5, 




IOtt 
~' 


S = {1,1,0}, P3 


= 3, 







which describe the two orbits, (2, 4, 8) x y and (6, 12, 10) x y , of order 3 as the origins for six different branches. 
Continuation to higher n is straightforward. 
For more general map f{x) = x'^ we have 

,d-l) = Gd- 
ji-i 

1=0 

In order to understand what happens when c 7^ we need to study the deformation of grand orbits with 
the variation of c. The deformation deforms the disc, but preserves the mutual positions of grand orbits on its 
boundary, i.e. the structure, partly shown in Fig l23l As soon as we reach the boundary of the domain A/i(/i) 
in Mandelbrot set, the attractive fixed point loses stability and this means that its grand orbit approaches the 
boundary of the disc - and intersects there with an orbit of order n (if / approaches the point A/i n A/„ S dMi). 
Every point of GOi merges with n points of GOn, which are located at different points of the disc boundary, 
so that the disc gets separated into sectors - the Julia set decomposes into a fractal-like structure, see Figsl^ 
and 123 When / enters the domain M„, the orbit 0„ gets stable and its grand orbit fills densely the interior of 
the disc, while Oi remains unstable and lies on the boundary, keeping it contracted (glued) at infinitely many 
points. Moreover, the angles a between the n components merging at contraction points increase smoothly from 
a = at Ml n M„ to a = — when the next bifurcation point M„ n Mi is reached. At Af„ n Mi new contractions 
are added - at the points, belonging to grand orbit GOi , which become merging points oi l/n sectors of our 
disc. And so on. 



(/fm+l = -(</?m + 27rs„i), S„ie(0,l,. 



(n) 



rf" - 1 ' 



5.4 Map f{x) = + c. Bifurcations of Julia set and Mandelbrot sets, primary and 
secondary 

The real slices of the lowest orbits and pre-orbits for the family fc{x) = a;^ + c are shown in FigEH In order to 
understand what happens beyond real section, one can study FigEHI 

In this case z{c) — c - this is the point with degenerate pre- image, which consists of one rather than two 
points; this is ramification point of the Riemann surface of functional inverse of the analytic function -f c. 
Actually, z{c) = c for all the families fc{x) = x'^ + c, but then it has degeneracy d — 1. This z(c) is /-image of 
the critical point w{c) =0 of multiplicity d — 1. 

For this family of maps one can straightforwardly study the pre-orbits and compare with the results of the 
resultant analysis - and this is actually done in Figs l28l2??l Shown in these pictures are: 

• the two fixed points (orbits of order one), i.e. zeroes of Gi{x) = x^ — x + c: 

X = — ± : 

2 2 

• their first pre-images, i.e. zeroes of Gi^i{x) — x'^ -\- x + c: 

1 VT^^ 
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Fi gure 23; Unit disc (Julia set at c = 0) with location of the first preimages (orbit and the first preorbit) of unstable orbits of 
orders one (circles), two (diamonds), three (boxes) and four (crosses). Locations of the points are shown for c = 0, but mutual 
positions of grand orbits remain the same for all complex c (though particular points can change order: a dot can pass through a 
cross — when c = Mi D M2, - but another dot will pass through the same cross in the opposite direction. Shown are also the lines, 
connecting the two points of the orbit of order two and the two points of its first pre-image - they will be among the contracted 
separatrices in Fig l24l — and analogous lines for an order-three orbit - to be contracted in Fig l25l 




Figure 24: Julia set at c = - |, i.e. when c = o-p n . See also FigEU 
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Figure 27: Julia set at c = — | + i-^, i-e. when c = o-j'^' D o-g"*'. This time the pairs of points from the grand orbit of order two 
should be strapped, as well as the sextets of points from a grand orbit of order six. 



^"li^^-S^:***"-": 

■■■■ 

^^^^^^^^^^ 


\ 

* 


— ■ 


* 



Figure 28: One more kind of section of Julia sheaf. The lowest two orbits of fc{x) = x'^ + c: of orders one (Gi(x;c) = 0) and 
two {G2{x;c) = 0) together with their first two pre-orbits {Gi^i{x;c) = 0, Gi^2{x;c) = and G2,i(a;;c) = 0, G2,2{x;c) = 0). The 
section Im c = 0, Im x = of entire 2-complex-dimensional pattern is shown, real c is on horizontal line, real a; — on the vertical 
one. Distinguished points on the horizontal line, where the roots of various polynomials Gn,s coincide, are: c = 1/4 - the zero of 
di{c) = 1 — 4c and thus also of -D(Gi^i) ^ di ; c = — the zero of •u'i(c) = c and thus also a zero of Wi^2ic) = wi{c)wi^2{c) and 
of D{Gi/2) ~ Wi/2', c = —3/4 — the zero of ri2(c) = 3 + 4c, thus vanishing at this point arc also discriminants D(G2), -D(G2.i), 
D{G2,2) and resultants x|i,2: ''i,l|2,l) ''i,2|2,2i proportional to r\2\ c = —1 - the zero of t02(c) = c + 1 and of the resultants 
^2\2,l} i'2\2,2^ ''2,i|2,2i proportional to 102; c = —2 — the second zero of 'uji,2(c) = c + 2 and thus of D(Gi,2) H^i,2. The zeroes 
of W2,2 = c^ + 1, where the two branches of G2,2 = intersect, are pure imaginary and not seen in this real section. w\^i and W2,i 
are identically unit for d = 2. 
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Figure 29: Instead of the section of the two-complex-dimensional pattern by the plane Im c = 0, Im 1 = 0, shown in Fig l28l 
this picture represents projection onto this plane, with Re c and Re x on the horizontal and vertical axes respectively. If it was 
projection of the same section, the only difference from Fig l28l would be attaching of horizontal lines to the right of the bifurcation 

points, corresponding, say, to Re ^§ ^ ^^"^" ^ = i for c > 1/4 (while there are no traces of these complex roots in the section, 

they are well seen in projection). If we now slightly change the section, from arg(c) = modTr to a small, but non-vanishing value 
(i.e. consider a section of Mandelbrot set not by a real axis, but by a ray, going under a small angle), then in projection we get the 

curves, shown in these pictures, a) The fixed points x = i it ^^2'^'^ • '-') '^^'^ same fixed points together with the orbit of order two, 

2; = — i it ^"1"^" . c) Fixed points together with their first pre-images, x = — ^ ^ ^^2'^'^ ' '^^ T'iie same with the second pre-orbit 
(zeroes of Gi,2) added, e) The same with the orbit of order two added. 
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(the lower branch of Gi.i = is mapped by / onto the upper branch of Gi — 0, (c, x) — (0, 0) is a fixed point 
of a family fc{x)); 



• second pre-images, i.e. zeroes of Gi^2{x) 



+ (2c + 1)2:2 ^ c(c+ 2): 




l + 2c yi-Ac 
2 ^ 2 



• the orbit of order two, i.e. zeroes of G2(x) = a:^ + a; + c + 1: 



X 



1 V-3 - 4c 

2 2 



• its first pre-image, i.e. zeroes of G2.i(a;) = a;^ — a; + c + 1: 



1 V-3 - 4c 

2 2 



. (this time the lower branch of G2,i = is mapped by / onto the lower branch of Gi = 0, (c, a:) = (—1,0) is 
not a fixed point of a family fdx)); 

• its second pre-images, i.e. zeroes of G2^2(a;) = + (2c — l)x^ +0^ + 1. 



From these pictures it is clear that along with the zeroes of ordinary discriminants and resultants (points 
c = 1/4 and c = —3/4), well refiected in the structure of universal Mandelbrot set, new characteristic points 
emerge: the zeroes of w„ and Wn^s (like points c = 0, c = — 1 and c = —2). Some of these points are shown in 
Figl3Ultogether with the Mandelbrot set. By definition these additional points are associated with the boundary 
of Grand Mandelbrot set. It seems clear from FiglJOlthat there is also an intimate relation between some of 
these points with the trail structure of the Mandelbrot set itself. In Fig l31l we show what happens to Julia set at 
zeroes of Wn and Wn,s- Clearly, zeroes of Wn.s mark the places where Julia sets change the number of connected 
components. 

Discriminants for lowest pre-orbits are given by the following formulas (note thatc?2 = 1 and all Wn,i — 1 for 




d= 2): 
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Fi gurc 30l Mandelbrot set, shown together with the zeroes of some Wn (crosses) and Wn,! (circles), associated with the more 
general structure of Grand Mandelbrot set. 
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5.5 Conclusions about the structure of the "sheaf of JuHa sets over moduH space 
(of JuHa sets and their dependence on the map /) 

1. Julia set J(/) is defined for a map fs E fi C A4. If / G o'n(A*), J{f) is an attraction domain of the 
stable /-orbit of order n. Unstable periodic orbits and all bounded grand orbits (including pre-orbit of the 
stable orbit) lie at the boundary dJ{f). Every bounded grand orbit is dense in dJ{f). This property allows 
to define dJ{f) even for / ^ M, i.e. outside the Mandelbrot set: as a closure of any bounded grand orbit (or, 
alternatively, as that of a union of all periodic orbits). Just outside Mandelbrot set this closure fails to be a 
boundary of anything ~ there are no stable periodic orbits and Julia set J(/) itself does not exist. 

2. If / G cTnil^) C Afi(/i), the Julia set J(/) is actually a deformation of unit disc. Its exact shape is 
defined by grand orbits of unstable periodic orbits, in particular all its "exterior points" /spikes are preimages 
of unstable fixed point (s). Its topology is even more transparent: dictated by position of / on the tree-powerful 
of the Mandelbrot set - as explained at the end of sec 15. 31 

In more detail: J(/) is a union of real domains in X (sectors of the unit disc): 

J(/)=U^J,(/). (67) 

In variance with the Mandelbrot set decomposition into elementary domains (T„'s, the constituents Jy (decom- 
position of the unit disc) and even their labelings v are not uniquely defined. 

2.1. At particular point (map) / G cr^^^ n cri^+^^ = ^4^' n 9cr^^'*'^\ where i?(Gfc,G„) = and two orbits 
Ok and On of orders n = mk intersect, every point of Ok coincides with exactly m = n/k merged points of 0„. 
Then every point of grand orbit GOk is a singular point of dJ{f), where exactly m components from the set 
{Jy[0„,Ofe]} "touch" together at zero angles, a = 0, see examples of {k,n) = (1,2) in Figl^ {k,n) = (1,3) 
in FigEHl {k,n) = (2,4) in FigE()l and (fc,n) = (2,6) in Fig[2I| There are no other singularities of dJ{f) at 

/G4^)nai^+^\ 

2.2. If we start moving / inside <j^^^ (k < n), every rn-ple from {Ju[0n,0k]}, which was touching at a 
single point for / G fife H cr„, acquires common boundary segments and the singularities get partly resolved: a 
point turns into m points and the full non-degenerate orbit 0„ is restored, see Figsl^andl^ On this side of 
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Fi gure 32: The deformation of Fig |24l Julia set in the vicinity of intersection point c = '^^^ ^1 cr!^^ with c CE ^[^^ • If is obtained 
from the unit disk by strapping the pairs of points on the boundary, which belong to grand orbit of the orbit of order 2. Domains Jv 
overlap and soparatrices connect the points of grand orbit of order 2, which are strapped at the intersection (phase transition) point. 

The length of separatrices can serve as an order parameter in this phase (ctj^^). Actually, there are infinitely many different order 

parameters inside cr^i^\ associated with transitions to different cr^^'. Separatrices are closures of peculiar (not-bounded) /-orbit 
and of its pre-images. 




It is 



Fi gure 34: Deformation of Fig l24l in another direction: Julia set in the vicinity of intersection point c = with c G 

In this phase the order parameter is the angle a > between the merging domains. It increases from a = at transition point 
c = cTj^' n cTj^' to the maximal value a = 27r/m = tt when c reaches other points of the boundary dcr^'^ , i.e. when the next phase 
transition occurs, e.g. when c = CTj^' P| fFig l26l or c = p| cr^^ (Fig l27i . 
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the (fc, rt)-" phase transition" the "order parameter" is the length of the common segments (separatrices) : the 
deeper in crj. the longer the segments. These segments can be also considered as peculiar unbounded grand 
orbit which densely fills "separatrices" between the merging orbits 0„ (on dJ{f)) and Ok (inside J{f)) when 
we approach the transition point from below - from the side of (k < n). 

2.3. If we start moving / inside ui^^^'' (n > fc), every m-ple from {Ji,[Om Ok]}, which was touching at a 
single point for / G crfc n (t„, continues to intersect at a single point - an orbit Ok is now unstable and should 
stay at the boundary dJ{f),- but angle a (the "order parameter" above the (A:, n)-"phase transition") is no 
longer vanishing, but increases as we go deeper into an, see FigsI^andOHl This is how the orbits Ok and On 
exchange stability: 0„ leaves dJ{f) while Ok emerges at dJ{f). Since n > k this is not so trivial to achieve 
and the problem is resolved by creating highly singular merging-points of m varieties at all the points of the 
grand orbit GOk, when it emerges at dJ{f). 

2.4. If, while moving inside cr^^^^'' {n > fc), the map / reaches a new intersection point an'^^^ H crj^^^\ 
N > n, our (A;, n)-structure (n/fc-merged points at GOk) feels this! Namely, when it happens, the angles ak,n 
reach their maximal values of our m — n/k merging varieties turn into needles in the vicinity of the merging 
points (i.e. at points of GOk), see FigsEHlandETI 

2.5. All this describes the behaviour of Ji, sets, associated with the subvariety a^^^ ncri^^^^' C dM{^) at the 
boundary between smooth constituents of Mandelbrot set. However, if we are at some point inside ak, there 
are many different structures of this type, coming from different boundary points (intersections with different 
Ckm with all integer m), and no one is distinguished. Thus the decomposition J(/) = U^Ji, is ambiguous for / 
inside ak and gets really well defined only exactly at dak ~ though is an everywhere discontinuous function on 
this boundary, changes abruptly as / moves along dak- 

2.6. If / crosses the boundary 9Af (/i) and goes outside the Mandelbrot set then there are no stable orbits left 
and nothing can fill the disc: only the boundary remains, the Julia set has no "body" - it becomes a collection 
of disconnected curves, which are no longer boundaries of anything (but are still densely filled by bounded 
grand orbits and/or unstable periodic orbits). This dissociation of the boundary occurs at discriminants points, 
where the two orbits of the same order, one stable and another unstable on the Mandelbrot-set side cross and 
become unstable beyond the Mandelbrot set. When they diverge after crossing, the orbits (and their entire 
grand orbits) pull away their pieces of the boundary of Julia set, and the "body" of Julia set "leaks away" 
through emerging holes, see FiglTHl 
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Despite it is not exhaustive, this decomposition structure carries significant information, in fact it is in- 
formation about topology of Juha sets and singularity structure of their boundaries. A remote region ai^^ in 
the main component of Mandelbrot set can be reached from ct^^^ by a path which lies entirely inside 

M(/i) and passes through a sequence of points ak H ct; connecting different components. Passage through every 
such separation point changes Julia set in above-described way, drastically changes its topology and angles a. 
However, the change of the shape (i.e. geometry rather than just topology) of J(/) when / is varied inside 
every particular component cr„ is not described in equally exhaustive manner. 

In other words, the suggested theory gives a complete description of the following characteristics of Julia 

set: 

- " Extreme" points (like " spikes" and " pits" in Fig ll7|l - they belong to grand orbits of unstable fixed points 
(orbits of order one); 

~ Merging points of sectors and needles - they belong to grand orbits, which became unstable on the path 
of / from inside the region a[^\iJ,) in Mandelbrot set Mi(/i), where the order-one orbit was stable. Every 
merging point is associated to particular "bridge crossing" at the point (T„ Hak in Mi{^), k\n, where the stable 
orbit Ok "exchanges stability" with 0„. The number of merging components is equal to m = n/k. The angles 
between merging components of Julia set depend on how far we are from the corresponding " bridge crossing" , 
and change from a = (touching) at the bridge to a = ^ (needle) when the next bridge is reached. 

3. The shape of dJ{f) is defined pure algebraically, because all the bounded grand orbits lie densely inside 
dJ{f). However, in variance with the boundary dM of Mandelbrot sets, dJ{f) can not be decomposed into 
smooth constituents (like 9tT„ which are smooth almost everywhere): no countable set of points from bounded 
grand orbits fill densely any smooth variety (like the zeroes of R{Gk, Gkm) did for given k and arbitrary m in 
the case of dak)- 

4. Unbounded (generic) grand orbits lie entirely either inside J(/) or beyond it. Each unbounded grand 
orbit originates in the vicinity of entire dJ(f): for every point of dJ{f) and for every grand orbit there is a 
branch which comes (with period d) close enough to this point. If grand orbit lies in J(/), it tends to a stable 
orbit - of order n if / G (7„. If grand orbit lies outside of J(/), it tends to infinity. One can say that infinity 
is a "reservoir" of stable orbits: if the degree d of polynomial map / is increased, the new stable orbits come 
from infinity. 

5. For / € Mna lying in other components of Mandelbrot set, disconnected from Mi (but linked to it by 
densely populated trails), the structure of Julia set is similar, with the only exception: the starting point is not 
a single disc, but a collection of discs, formed when the map / leaves Mi and steps onto one or another trail 
leading to the other components Mna- Both periodic orbits and pre-orbits jump between the boundaries (if 
unstable) or interiors (if stable) of these discs. With this modification, the same description of Julia sheaf is 
valid, with " monodromies" described by the same procedure of interchanging stable grand orbits inside with 
unstable grand orbits on the boundary, which causes fractalization of the boundary. 

6. If / ^ M{iji) gets outside the Mandelbrot set, the Julia set is no longer a disc or collection of discs: it looses 
its "body" (which would be occupied by the stable bounded grand orbit) and consists only of the boundary 
(where unstable periodic orbits and their grand orbits live). Its bifurcations can still be described algebraically, 
but now the intersections of unstable orbits will be the only ones to play the role. These bifurcations are 
controlled by the Grand Mandelbrot set. 
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6 Other examples 



In this section we present some more examples of the few-parametric famiHes of maps fc{x) in order to show 
that our considerations and conclusions are in no way restricted to the peculiar case of fc{x) = x'^ + c. For every 
family /i C we describe the Mandelbrot set M{fj,) as the union of stability domains M„(/i), its boundary 
dM(fi) as the closure of appropriate discriminant variety, the structure of grand orbits and of algebraic Julia 
sets. 

We are going to illustrate and support the following claims: 

• Consideration of various 1-parametric families shows universality of the tree structure of Mandelbrot sets 
and non-universality of particular numbers: shows the difference between n, fc, m- and a-parameters. 

• Particular 1-parametric families (like {fc{x) — x'^ + c} or {/(a;) = cx"^ + x'^}) are not in generic position - 
even among the maps of given degree - and the corresponding sections of Mandelbrot /discriminant variety are 
not fully representative: have additional special features, symmetry properties and accidental degeneracies. 

• 2-parametric families can be considered as families of 1-parametric families. Their analysis helps to 
understand how the non-universal components of the description change and helps to understand their nature 
and relation to next-level algebraic structures like singularities of resultant varieties, controlled by the resultants 
of the higher order. Clearly, non-universality is no more but the corollary of our restriction to ordinary resultants 
and further work will provide unique interpretation to all properties of particular sections/families ^. 

In variance with the previous section [S] we do not give exhaustive treatment of every particular example 
below. Each example serves to illustrate one or another particular aspect of the problem and a particular line 
of further development. 

6.1 Equivalent maps 

Different maps / and / can be equivalent from the point of view of our considerations if they are related by a 
diffeomorphism 0ofX, f o <p — <j> o f, or 



Then, obviously, /°" o cj) — (p o /°", Fn o (j> — 4> o Fn, G„ o (j) = (j) o Gn and so on. If / and / are polynomials 
of the same degree - what will be the most important application of equivalence below - than cj) should be a 
linear transformation, (f>{x) = ax + b. 

6.2 Linear maps 

6.2.1 The family of maps fa/i — a + /3x 



f{c^{x))^c^{f{x)). 



(68) 



In this case 




m 



l)[x + 




and decomposition formula H20|l reads: 



r(n) 



Fn — Gfe, 



k\n 



Gi{x\ fap) = (/3 - l)x + a, 
Gk{x; fap) = gkiP), k> 1. 



(69) 



The functions gkiP) are circular polynomials: 



r(n) 



k\n 
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with 



For simple p > 2 



.91 =/?-!, 
52(/3) + ! = -.91 (-/?), 
ff3(/3) +1, 
34 (/?) + 
55(/3) = /3' + /3' + /3' + /3+l, 
56(/3) = /3'-/3+l = .93(-/3), 



(70) 



52p(/3) = = P"-' - P"'' + . . . - /3 ■ 

9p^m = ^^lr_\ = /^''^''-'^ + . . . + /3P + 1 

and so on. 

All Gk{x; fais) in H()9(l except for Gi do not depend on x. This means that there are no orbits of order 
fc > 1 unless P is appropriate root of unity - when there are infinitely many such orbits. For example, when 
P = ~1 any pair (x, a — x) form a second-order orbit of the map x a ~ x] and when /3 = ±j any quadruple 
(x, ±ix + a, —X + (1 ± i)a, —ix ± ia) is an orbit of order four etc. 

Of course, the family of linear maps is a highly degenerate example, still this peculiar section of Mandel- 
brot/discriminant variety should be reproduced in particular limits of more general families of maps below. 

6.2.2 Multidimensional case 

We use the chance to illustrate the new ingredients of the theory which arise in multidimensional situation. 

Take X = C and consider the linear map / : x ~f Bx where i? is an r x r matrix. The iterated map 
remains linear, f°^{x) = B"-x and F„(a;) = (B" — I)x. Irreducible components Gk = gk{B) are independent 
of X for fc > 1, just as in the one-dimensional situation. However, for r > 1 degeneration pattern is richer: 
of interest are all situations when a pair of eigenvalues coincide, matrices acquire Jordan form and some two 
eigenvectors become coUinear. 

The case of r = 2. The eigenvector x = (xi, X2) of _B = ( l*^^ ) is a quadric in CP"*" in homogeneous 

\ 021 022 J 

coordinates, i.e. satisfies 

]^X^Qij{B)x.j = 612X2 + (611 - 622)xiX2 - 1^21x1 = 0. (71) 

The linear map B is degenerate when determinant of the matrix Q{B) (or, what is the same, discriminant of 

lijXj I 



the polynomial XiQijXj / x\ in xxjxi vanishes). 



det Q{B) ^ (611 - &22)' + 4612&21 = 0. 

We call this expression discriminant of the linear map B. For r > 2 consideration is a little more sophisticated 
and we do not discuss it here, see J^l- The resultant of two maps A and B vanishes when some eigenvector of 
A gets collinear to some eigenvector of B. Of course, different iterations of the same linear map have the same 
eigenvectors, and resultant analysis for linear maps remain trivial. Still, it is important to remember about 
these structures in analysis of non-linear multidimensional maps. 

6.3 Quadratic maps 

This class of examples helps to demonstrate, that even equivalent maps (related by diffeomorphism in M) can 
have differently looking Mandelbrot sets. Of course, this is not a big surprise: non-linear change of variables 
(say, c ^ c^) in Mn fi changes the shape and even the number of domains in FigtTsl 
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6.3.1 DifFeomorphic maps 



Consideration of generic 3-parametric family of quadratic maps fap-fix) = a + Px + jx'^ can actually be reduced 
to that of fc{x) = + c by the rule (jHEjlj because 



provided 



ifap^ix) + I = /c(72; + |) 



«7+f-^ (72) 



Then ir^l^(x) + I = /c(7/a/37(a^) + f ) = hifdix + f )) = ffi^fx + f ) and, in general, 



7a,(x) + | = /r(7a:+f). 



In other words, /a/j-y : M. ^ Ai\s difTeomorphic to /c : Al. Discriminants and resultants for /a/37 

obtained from those for /c by a substitution H72|l . However, since H72I) is non- linear (quadratic) transformation, 
it can change the numbers of zeroes, and the numbers of elementary domains in sections of V* can be different 
for different 1-paramctric sub-families of {fafi-y}- 

For two 2-parametric families to be analyzed below we have equivalences: 

2 n Aq + 2p-p^ 
X +px ^ q: p = p, q ~ a-f, c — 

and 

1 - 6^ 

72;^ + (6 + l)a; : p = b+l,q = 0, c= — - — . 



6.3.2 Map / = + c 

First of all, we summarize the discussion of Mandelbrot set in the previous section |S1 in the form of a table: 



/ = a;2 + c 


n 


1 2 


3 


4 


5 


6 


# of orbits 
= <ieg^[Gn]/n 


2 1 


2 


3 


6 


9 


# of el. domains 
= deg^[G„(iOc,c)] 


1 1 


3 


6 


15 


27 


totcil ^ of cusps 
= deg^ [d„] 


l-4c i 
1 


-7-4c 
1 


64c3 + 144c2 + 
3 


see eo.lKOI 

11 


see eg. 1611 

20 


total #: of {n,n/m) 

touching points 
= deg^[T-„ for 

m=2 


3+4c 

1 
1 


- 


T'24 
4c+5 

1 
1 


- 


64c3 + 128c2 + 
+72C+81 

Q 


m=3 




16c2 + 4c + 7 
2 




- 


16c2 + 36c + 21 
2 


m=4 


- 


- 


ri4 

16c2 - 8c + 5 
2 


- 


- 


m=5 








256c* + 64c^ + 
+16c2 - 36c + 31 
4 




m=6 


- 


- 


- 




ri6 

16c2 - 12c + 3 
2 


■^n'^ '■ # of el. domains 
with q(p) cusps 












p = 0, q=l 


1 


1 


3 


11 


20 


p=l, q = 


1 


2 


2 


4 


3 


p = 2, q = 








1 





4 


p = 3, q = 


















= (4c)6 + {4cf + (4c)4 + (4c)3 + 15 • {icf - 17 • (4c) + 127, 
ri8 = 2560"* + 32c2 - 64c + 17, 



and so on. Stars in the table stand for too long expressions. Afn is the number of domains (T„ of given level 
p and order n, which still differ by the values of parameters rrii and at that are ignored in this table. For 
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one-parametric families the number q of cusps often depends only on the level p. Taking this into account, we 
can write down - and check in this and other examples - the set of sum rules for the numbers 



^AAif)=degJG„K,c)], 

oo 

^g(p)AAif)=degJd„], 

oo 

=J2degArn,n/m]- (73) 

p—1 m 

Actually, for every given n only finitely many values of p contribute: p is the number of non-unit links in the 
branch of the multipliers tree, thus p < log2 n. From H73|l one can deduce an identity which does not include 
sums over p: 

] = degc[G„(wc,c)], (74) 

m 

and in many cases (including the families x'^ + c) the r.h.s. is just Mn{d). Also in many cases (again including 
x'^ + c) (7(0) = d — 1 and q{p) = d — 2 for p > 0. In such situations 

[d - 1)AA(0) + (d - 2) E degJr„,„/„] = degjd„], (75) 

m 

and we also have a consistency condition between H74|l and (|75|l : 

degjd„] + Edegc[r„^„/™] = {d- l)degJG„('u;c, c)]. (76) 

m 

Coming back to our particular family -I- c, one can easily check that consistency condition is satisfied by the 
data in the table, and also Mn^ = deg^[dn], as requested by (|75|) . 

For analysis of Julia sheaf above table should be supplemented by the data concerning pre-orbits and 
associated resultants, but this remains beyond the scope of the present paper. 

6.3.3 Map /-y^o = 72;^ = -I- (6+ l)x 
According to H72|) we have 

c=^-^ = i-(^ = i(l-0, (77) 
2 4 4 4 4^^' ^ ' 

and one can expect that quadratic cusp for domain cti at c = -j disappears after such transformation. In fact, 
transformation 177|) makes from cardioid cti in c-plane a pair of unit discs in 6-plane with centers at 6 = and 
6 = 2, which touch each other at the point 6 = 1, S": = cr^'^-' [l-|-]Ucr*^°^ this is a good example, demonstrating 
that links with m = 1 in multiples tree sometime contribute to the structure of Mandelbrot set. 
The boundary dSi is described by a system of equations H38|l : 



dSi : 



\2-ix + l3\ = 1 
72;^ -\- I3x — X 



i.e. X = OT X = \f3\ = 1 or 1/3 — 2| = 1 and |1 ± &| = 1- Accordingly the resultants ri„(&) should have 



1-/3 
7 ' ^ ^ 

roots lying on these two unit circles. Actually, they are all made from products of circular polynomials H7U|) : 

ri„(6) = 7^"g„(/3)g„(2 - f3) = 7^"5«(1 + b)gn{l - 6), 

r2,2„ = 7'^'"5n(5 - b^), 
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in particular (for 7 = 1) 



l-^n2{b) - 3 + 4c = (/3 + 1)(3 ~ /?) = (2 + b){2 - 6), 
7-6^13(6) = 16c2 + 4c + 7 = + /? + l)(/32 - 5/3 + 7) = (3 + 36 + b^){3 - 36 + 6^), 
7-1^14(6) = 16c2 - 8c + 5 = + - 4/? + 5) = (2 + 26 + 62)(2 - 26 + 6^), 

7"^°fi5(fc) 256c4 + 6Ac^ + IBc^ - 36c + 31 ^ 
= {p'^ + (3^ + (3^ + 13 + - 9(3^ + 31/32 _ ^ 
= (5 + 106+ 106^ + 56^ + 6'')(5 - 106+ 106^ - 56^ + 6^), 
7-54^16(6) = 16c2 - 12c + 3 = (/32 - /3 + l)(/32 - 3/3 + 3) = (1 + 6 + 62)(1 - 6 + 6^), 

7"^^?-24(6) =4c + 5==6-62 = l + (5- 6^), 
7~^'*r26(6) = 16c2 + 36c + 21 = 6^ - 116^ + 21 = (5 - b^f + (5 - 6^) + 1, 

7"^°*?-36(&) = 64c3 + 128c2 + 72c + 81 = -6^ + 116"* - 376^ + 108, 



The last expression, for r^Q, has nothing to do with circular polynomials: the map H77|l transforms crj , which 
is quadratic cardioid in c-plane, into a bouquet of two unit discs in 6-plane; it converts Uj^' , which is a disc in 
c-plane, into domain 16^ — 51 < 1 in 6-plane; but it does not simplify higher ct„, which are complicated in c-plane 
and remain complicated in 6-plane. 

What is interesting in this example, dependence on 7 is trivial in 6-plane, and the pattern of Mandelbrot 
set (its section by quadratic maps) remain the same in the limit of 7 — > 0. In this limit the biggest part of 
Mandelbrot set is associated with orbits, lying at very big x ~ 7^^- Naive limit 7 = 0, considered in above 
s l6.2l ignores such orbits and only the unit disc |/3| < 1 (domain a^"^) is seen by examination of linear maps. 
However, the other parts are actually present (nothing happens to them in the limit of 7 = 0), just not revealed 
by consideration of strictly linear maps. Similarly, consideration of quadratic maps alone can ignore (overlook) 
other pieces of Mandelbrot set, and so does every restriction to maps of a given degree. 

Powers of 7 in above formulas for resultants characterize intersection of orbits at infinity. Since in this 
particular case all periodic orbits except for a single fixed point tend to a; = 00 as 7 —> 0, the asymptotic of 
resultants is defined by the following rule (see ea. H80|l below): 

i?(Gfe,G0-rL~7^^(^'-^'-), k>l, 

where Ni — Ni{d ^ 2) is the degree of the polynomial Gi{x; f) for quadratic map /, and Ni/l is the number of 
periodic orbits of order /. 

We omit the table for this family, because all the numbers in it are obtained by doubling the corresponding 
numbers in the table for + c (provided degrees of all polynomials are counted in terms of 6 rather than c). 

6.3.4 Generic quadratic map and f = x'^ + px + q 

Expressions for generic quadratic map jx'^ + l3x + a can be obtained from formulas below by a trivial change 
of variables: x — *■ jx, p = P, q = aj. 

Gi ^ Fi ^ x^ + ij) ~ l)x + g, 
G2^^^x^ + {p+l)x+{p + q+l), 

Gi,i = = G2 - 1 = a;2 + (p + l)a; + (p + g + 1), 

Gi,2 = ^x'^ + 2px^ + {2q + p^+p+ l)x^+ 

GiGi^i 
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+ {2pq + p2 + p)x + {q^ + pq + 2q+p), 
p„ 1 

G2,i = „ „ „ ^x^ + {p- l)x + {q + 1), 

Lrl{_z2Lrl,l 

^ ^2^2 ^ 

2,2 

CTiLr2(-Ti_i(jri_2(-T2,l 

= a;'* + 2pa;3 + (p^ + 2^ + p - \)x^ + (p^ + 2pq - p)x + (g^ + + 1) 

The critical point is w = —p/2^ the "irreversible point" z = f{w) = ^^j^- Accordingly, for generic map 
w = -/3/27 and z = ^22_£!. 



UI3 



wi = 9 


-I 


4 ' 






W2 = 


1 + 


Wl, 






= 1 + 2p(7 + g + 4 




p' 


2,p^q^ 


ipq^ 




4 


4 ^ 


2 








3p'*g 




2 4 16 




~ 64 


+ 16 " 


^ 32 ' 



Wll = 1, 

u;i2 = 2 + Wl, 



6.3.5 Families as sections 

We now use these examples to briefly discuss the geometric interpretation of the theory. The Mandelbrot set in 
Fig 1 151 is obtained as a union of infinitely many elementary domains a in the plane of complex c, glued together 
in a very special way along the discriminant variety (which has real codimension one and is a closure of a union 
of complex varieties 7?,^ „ of complex codimension one). Each domain a can be described as a real domain, 
defined by some equation of the form E(c) < 1 with a real function S(c) of a complex variable c. 

Switching to the two-parametric family / — x^ + px + q, we get the Mandelbrot set, which is a similar union 
of domains a in the space of two complex variables p and g, with each a defined by the equation g) < 1, 
where essentially q) = i^c — q + \ — j{p — 1)^) . 

Since we can not draw pictures in C^, we restrict the whole pattern to real p and q: this badly spoils the nice 
picture, but preserves the main property which we are going to discuss. After restriction to real c, Fig llSl turns 
into collection of segments (implicit in the lower part of that figure), see Fig|2H Instead we can now extend the 
picture in another direction: to two-dimensional plane {p,q). Segments turn into parabolic strips, see Fig l37l 
where different families, like x'^ + c and x^ + [b + l)x, are represented as different sections. It is clear from the 
picture why one and the same Mandelbrot set looks very different if restricted to different families. Another 
version of the same pattern is presented in Fig l38l where Mandelbrot consists of straight rather than parabolic 
strips (after complexification strips become "cylindrical" or, better, toric domains). However, beyond quadratic 
family such ideal "torization" is not possible, different resultant varieties can not be made exactly "parallel", 
because sometime they intersect, see eas. (|79f) for the simplest example - at level of cubic families. The best one 
can try to achieve is the representation in terms of somewhere-strapped tori, represented by chains of sausages, 
see Fig0Blat the end of the paper. It is an interesting question, whether Mandelbrot set can be represented as 
a toric variety. 
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Figure 36: The Mandelbrot set for the family a;^ + c shown together with its section by a real line Iin c = 0, where it turns into 
coUeetion of segments. These segments form disconnected groups — sections of different components M„ (actually, only the big 
group, associated with Mi, and a tiny one, associated with M3, are seen in the picture; if the section was not by a real line, a single 
Ml would provide many disconnected groups of segments). Within each group the segments, corresponding to different touch 
by their ends, real zeroes of resultants, and free ends are real zeroes of discriminants (perhaps, those of pre-orbit polynomials, i.e. 
zeroes of dn or Wn,s)- 
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Fi gure 38; The same picture in the real (c, p)-plane: p plotted on absciss axis, c — on ordinate. The Mandelbrot set is represented 
by collection of straight horizontal strips. Resultants and discriminants depend only on c and not on p thus their zeroes form vertical 
lines, separating the strips. In general situation (beyond quadratic families) the Mandelbrot set can not be made "cylindrical": 
the obstacle is provided by intersections of different resultant and discriminant varieties, controlled by higher resultants and 
discriminants, see 1791 for examples. The best one can hope to achieve is a "sausage-chain" representation, see Fig l4(il at the end 
of this paper. 



6.4 Cubic maps 

Generic cubic map - the 4-paranietric family faf3-ys{x) — a + ^x + ^x^ ~\- 5x^ - is difFeomorphic to a 2-parametric 
family fp.q{x) ^ x^ + px + q, 

27(5 

In addition to some initial information about this 2-parametric family we provide a little more details about 
several 1-parametric sub-families of cubic maps, 

fc[x) = .T^ + c with (p, q) = (0, c), 



3c' 27c Vcj ' 

3 ^ 2 ■ 1 . x (1 - a)^ 2(1 - a)^ -f- 27a^ 

fdx) = ax^ + 1 - a)x^ + c with {p, q) = 

\ 3a nay/a 

respectively (in the last case a is considered as an additional parameter). 

Our main purpose will be to demonstrate four (inter-related) new phenomena, not seen at the level of 
quadratic maps: 

- Keeping degree of the maps fixed, we ignore (hide) the biggest part of Julia set, which is detached from 
the "visible" part and is located at x = 00. Julia set decays at peculiar "decay" bifurcation points, which are 
presumably are zeroes of Wn,s[f)- See Fig|2] 

~ For essentially multiparametric families p, of maps (like the 2-parametric fp,q{x) — x^ + px + q), different 
components Vkn and dn of the resultant/discriminant variety intersect and higher order singularities of this 
variety can be revealed in this way. At these singularities particular sections of universal Mandelbrot set, 
associated with 1-parametric families, reshuffle, and new components Af„Q, can split from previously existing. 
Thus these reshufflings are responsible for formation of the trail structure of particular Mandelbrot sets. 

- If the Zd symmetry of the map f{x) = x'' -I- c is slightly, but fully broken (e.g. by addition of a term 
px or with small non- vanishing beta or 7), then the Zj-i symmetry of the section of Mandelbrot set gets 
broken in peculiar way: the central domain Mi and the nearest (in a relevant topology, actually, those with 
largest sizes p„) components M„q, continue to possess the Zd-i symmetric shape, while remote components 
Mna acquire the non-symmetric form (characteristic for the Mandelbrot set for the family x^ + c). See Figl39l 
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Figure 39: Different, close and remote, components of the Mandelbrot set in the complex c plane for the family ax^ + (1 — a)x + c 
for different values of a (considered as an additional parameter): a) a = 4/5, b) a = 2/3. One can see how the symmetry of different 
domains is changing with the variation of a: the smaller a the larger is the low-symmetry domain. See Fig |4Ul for description of 
transition (reshuffling) processs and s lb.4.5l for comments. Halo around the domain in the right pictuires (at high resolution) is an 
artefact and should be ignored. 



Detailed analysis of reshuffling process, shown in Figsl40l and 1411 reveals how the non-central components and 
trail structure of Mandelbrot set is formed. 

- In generic situation a given map (at a given point in the Mandelbrot space) can possess several stable 
periodic orbits (not just one as it happens in the case of x'^ + c - and this brings us closer to continuous 
situation, where several stable points, limiting cycles or other attractors can co-exist). Bifurcations of Julia set 
occur whenever any of these orbits exchange stability with some unstable one. See Figl42l 

6.4.1 Map fp^q{x) ^ + px + q 

Gi{x) ^x^ + {p-l)x + q, 
G2{x) ^x^ + {2p + l)x^ + 2qx^ + {p^+p+ l)a;2 + ((2p + 1)^ - l)x +p + q^ + l, 

G,,jix) = G2-l, 

The two critical points are w = i^/— p/3, or p — — 3w^. Here are some irreducible components of pre-orbit 
resultants and discriminants (they are all products of two factors, associated with +w, and —w, such products 
are polynomials in integer powers of p; we list expressions which are not too long): 

wi = (q- 2w^ - w){q + 2W^ -|- to) = ^(27cj^ + 4p^ - 12p^ + 9p), 
W2 = (4«)^ - 210** - 4w'-^q - 2w^ + wq + q'^ + l)(4to^ - 2w''' + 4w'-^q - 2to^ - wq + q^ + 1), 
= (g - 2w)^ + 2w){q + 2w^ - 2w) = ^(27g^ + 4p^ + 24^^ + 36p), 
«,2,i = (g^ + 1 - 2'wq + w'^ - iw^q + Aw^ + Aw'^){q'^ + 1 + 2wq + w'^ + 4w^q + Aw'^ + 4u)^), 

Some discriminants and resultants: 

di = -{27 q^ + 4p3 _ I2p2 + i2p - 4), 
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Figure 40: Shown is a sequence of views of Mandelbrot sets for different values of a in the 2-parametric family {ax^ + {I — a)x^ + c} , 
interpolating between the 1-parametric families {x^ + c} at a = 1 and {x^ + c} at a = 0. See s l6.4.5l for additional comments and 
footnote 1161 in that section for precaution concerning this particular Figure: only the right parts of the pictures can be trusted, 
the left parts should be reflections of the right. Arrows show positions of some zeroes of di, d2 and r\2- The full sets of these zeroes 
are shown in fully reliable (though less picturesque) Fig |41l 
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Figure 41: A sequence of pictures, showing positions in the complex c-plane of zeroes of di{c) (big cirlces), d2{c) (big crosses), 
d3{c) (big boxes), ri2(c) (small crosses), ria (c) (small boxes) for the 2-parametric family {ax'^ + (1 — ci)x^ -\-c}, interpolating between 
the 1-parametric families {x^ + c} at a = 1 and {x^ + c} at a = 0: a) a = 2/3, b) a = 1/6, c) a = 1/15. These are particular 
special points in the patterns, presented in Fig |40l Three arrows point to a zero of di, a zero of ri2 and a zero of d2. After the two 



roots of ri2 merge at the real-c line, this happens at a ■ 



5-V2T 



1/5, they start moving in opposite directions along this line. 



We arbitrarily choose to point at the right of the two zeroes (because the left one is not seen in the partly erroneous Fig |4Ui . 
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Figure 42: Coexisting stable orbits for the family cx^ + x^: a) c = —0.22, b) c = —188, c) c = —0.178. Shown are Julia sets, and 
stable orbits lie inside J(/). In the left hand part of the Julia set there is a stable fixed point. The stable orbit inside the right 
hand part is different: shown is the exchange of stability between another fixed point (stable in the case (a)) and an orbit of order 
two (stable in the case (c)). In (b) the transition point is shown, where the fixed point and the orbit of order two intersect and 
Julia set changes topology. 
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= 27g2 + V - 1 V + 16, 
d2 = 27 + Ap^ + 24p2 + 48p + 32, 
mil = 27(g + w + 2w^)(gf -w~ 2w^) = 27wi, 
Z?(Fi,i) _ I?(Gi,i) 



dl 



= -27(g - 2w + 2w^)((7 + 2w - 2w^) ^ -27u;i,i, 
''111, 2 = 27(g + w + 2w^)(9 - w - 2w^) = 27u;i, 
^i,i|i,2 = R{GiA,GiaY'^ ^ 27{q + W + 2w^){q -w~ 2w^) = 27u;i, 



Zeroes of wi,i(/) form a submanifold, where the Julia set can decay into two disconnected components 
{{p,q) = (—3/4,3/4) is a particular point on this submanifold, another point is {p,q) = (0,0), where Julia set 
turns into an ideal circle.) As cubic / degenerates into quadratic map, namely, when q with p fixed, one 
of these two components travels to x = oo. See Figl31l 

Since in this subsection we are dealing with a two-parametric family, we can observe intersections of various 
components rk^n and d„ of the universal discriminant variety (these intersections lie in complex codimension 
two). They are seen in the section by complex-dimension-two family f = +px + q as particular points: zeroes 
of the resultants of d{p, q) and r(p, g), considered as polynomials of q for a fixed p or vice versa. For example: 

Riri2,di\q) = 16•36(3p-5)^ 
R{ri2,d2\q) = 16 -3^(3^ + 2)2, 
R{dl,d2\q) = 16-3i°(p2+p+l)2; 

(78) 

i?(ri2,di|p) = -2i°(72V + 32), 
i?(ri2,d2|p)- 28(72^ + 256), 
R{di,d2\p) = 2*3^(27^^ + 16). 

6.4.2 Map = + c 

Gi ^ — X + c, 
G2 = x*^ + x"^ + 2cx^ + x^ + cx + + 1, 
Gi^i = G2 — 1, 

The double stable point is w = 0, and 

W2 = {l + c'r, 
W3 = (l+c2 + 3c4 + 3c6 + c*)2, 
W4 = c*, 

Wll = c^, 
«;i2 =c2(c'^ + 3c^ + 3)^ 

W2I = (l+c2)2, 
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W22 = l+ic^ + '2c^+ic}'+ic}^+c}\ 

W31 = (1 + + 3c^ + 3c^ + c^f , 

These quantities are obviously factorizable - in seeming contradiction with our claim that they are all irreducible. 
In fact, factorization is accidental and is lifted by infinitesimal deformation, for example, by taking p ^ 0, see 
the previous subsection 16.4.11 Such deformation also breaks the accidental coincidences, like that of W2 and 

W21- 



f = + c 


n 


1 


2 




3 




4 


5 


# of orbits 
= deg^[G„]/n 


3 


3 




8 




18 


48 


# of el. domains 
= deg^[G„(tOc, c)] 


1 


2 




8 




24 


80 


total # of cusps 
= degjd„] 


4 - 27c2 
2 


(27c2 + 
2 


32)i 


* 

12 




* 

40 


? 
152 


total #: of {n,n/m) 

touching points 
= de£ for 
















m=2 


- 


7^12 

27c2 + 
2 


16 


- 


729c* 4 


- 1620c2 + 1000 
4 


- 


m=3 


- 






729c'' + 27c2 + 169 
4 






- 


m=4 










729c'* 


ri4 

- 324c2 + 100 
4 




m=5 














312^8 + 394^6+ 
+3<^46c* - 7 • 27 ■ 263c2 + 11* 
8 


■^ri^ '■ # of cl. domains 

a<f ' c s„ 

with g(p) cusps 
















p=0, q=2 


1 







4 




16 


72 


p=l, q=l 





2 




4 




4 


8 


p=2, q=l 















4 





P=3, q=l 






















92 



da = 3^^c^^ + 2 ■ 3^^c^° + 3"3571c* + 2''3^11 ■ 19 ■ 107c^ + 
+25 3**13^137c'' + 2''3''l3''l7c2 + 2*13^^. 

d4 = (150094635296999121c2'' + 1200757082375992968c2^ + 
+4203267461712259335c2° + 8399740516065395253ci* + 
+10909964351274746418c^'^ + 10526401881511556976c" + 
+8522156414444085612c^2 +55446117194182680000^° + 
+2750472027922567500C** + 1314354779366400000c« + 
+459901255680000000C'' + 167772160000000000) • 
■(282429536481c^^ + 1757339338104c" + 4642459719687c^2 + 
+6806074010589ci° + 6891783220746c® + 5994132959232c® + 
+41 18269132800c'' + 17394617548800^ + 1073741824000) = 
= 33«c24 + ... = (27c2)i2 + ... 

c?4 appears factorizable, but this is an accidental factorization, lifted by any deformation of the symmetric 
family + c. 

The sum rules H73|l now imply: 

■^n^ + degc[''«,n/m] = degJG„('u;c, c)], 

m 

,n/7n 

5]AA('')=^deg,[r„,„/„] 

p>l m 

and the consistency condition 

degcK] +^degc[r 

n,n/m, 

m 

all obviously satisfied by the data in the table (question marks near some numbers in the table indicate, that 
they were found not independently, but with the help of the sum rules). 

With these explicit expressions we can illustrate the small general theorem, formulated and proved after 
ea. H41(l . Namely, solutions to the system 

where Hn = {-F^, Gn}, are given by zeroes of irreducible discriminant dn- 

For n = 1 we have: G'l = — 1, the system has solutions, i.e. G[ and Gi possess a common zero, when 
their resultant D{Gi) = 4 — 27c^ vanishes, and Hi = —Qx never has common zeroes with G'^. Thus solutions 
of the system are zeroes of the irreducible discriminant di, which in this case coincides with D{Gi). 

For n = 2 G'2 = 6x^ + 4a;^ + 6cx^ + 2x + c and H2 — — 36x(a:^ + x + c){x^ + c) {x^ — x + c) both depend on c and 
have a common root whenever their resultant i?(iJ2, Gj) = -2'^'^3'^^c'^{27c^ - 16)(27c2 + 16){27c'^ + 8) vanishes. 
At the same time, the system can have solution only when G2 and G2 possess common zeroes, i.e. when 
D{G2) = {21(? + 16)(27c^ + 32)^ vanishes. We see that a pair of the roots, 27c^ + 16 is actually eliminated, 
because they are also zeroes of R{H2,G'2) and what remains are the roots of the irreducible discriminant 
d2 = 27c2 + 32. 
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Fi gurc 43) A fragment (fibers shown at several points only) of the Julia sheaf for the family -\- c — the analogue of Figs ll6l 1171 
for a;2 + c and|^for + c, a) c = + 0.7i, b) c = + 0.9i, c) c = 0.2 + l.OSi, d) c = 1.14i, e) c = 0.45. 



The first few pre-orbit discriminants for tfie family x'^ + c are: 



ns D{G ns) 

11 -i^dlwii 

12 -i'^'^ d\wiwl^wi2 

13 -i^^^df-wiwl^wl^wiz 

14 -3567^54^13^,27^,9^^3^^^^ 

21 36d|r?2^«2i 

22 ■i^^dfrl^wl^W22 

31 324^6^4^^^^ 



The Jiiha sheaf for this family is shown in Figl43l 
6.4.3 Map fc{x) = cx^ + 

Gi = cx"^ + — x, 
G2 = e'x^ + 2c2a:5 + c(c + l)a;'* + 2ca;^ + (c + l)j;2 + x + 1, 
Gi_i = G2 — 1, 

Critical points w — {0, ^^}- Contribution of critical point w = can be ignored in most applications, and 
we define in this case: 



Wi — Gi{w) — Fi{w) — 
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^^•^= G^(^ = ~3^(,^"9j■ 
c = 4/9 means that q) = (-3/4, 3/4). 



/ = cx^ + 


n 


1 


2 3 


4 


^ of orbits 
= degx [Gn]/n 


3 


3 8 


18 


# of el. domains 
= degc[G„(£,c)] 


1 


3 12 


36 


total # of cusps 


4c + 1 

1 


ic5/2(4- 13c + 32c2) * 
2 10 


* 

.■^1 


total # of (n,n/m) 

touching points 
= degc[r„,„/^] for 








m=2 


- 


c5(6c + 3) 
1 


''■24, 

ci°5(i00c3 - 220c2 + 47c + 20) 
3 


m=3 




c22(169c2 + 68C + 7) 
2 




m=4 






rn 

c™(100c2 +44c + 5) 
2 


A/'n''' : # of cl. domains 

(Jji l— On 
with q(p) cusps 








p=0, q=l 


1 


2? 10? 


31? 


p=l, q=0 





? ? 


? 


p=2, q=0 





? ? 


? 


p=3, q=0 





? ? 


? 




d3 = 


c''''(1235663104ci0 - 765891776c^ + 315356704c** 





-107832976c'' + 33146817c^ - 9493768c'' + 2585040c'* + 
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+233040c3 - 123072c2 + 42752c + 7168) = 
= c66(2»136ciO + ...) 

d4 = c'^°°(+180143985094819840000000000000c^i - 39631676720860364800000000000c=*° + 
+8693002811987722240000000000c2^ - 13044663126539108352000000000c28 + 
+1173757802691870851072000000c2^ - 251521821595369943859200000c26 + 
+571799816500523292950528000c2^ + 31829281120908656668016640c24 + 
+4498587759541291509802744c2^ - 223935602294557922308802760^2 - 
-4844558197092071004917119c2^ - 391359614494908199135640c2° + 
+847227219438292501416048c^^ + 319824545156577533637200ci* + 
+43548434837473605814592c^^ + 33835232982782581502464ci'' + 
+6794824688429177434112ci5 - 295875264228762906624ci* - 183114267457753161728ci3 + 
+165660797805589659648ci2 _,_ 93048634733415038976c" + 19011408301729447936ci° + 
+25910347327455363072c^ + 16153160956010037248c* + 4359984836267474944c^ + 
+933600163325280256c'' + 3141248453967872000^ + 94391968968212480c'' + 
+17132514753642496c^ + 17820882315509760^ + 99381248262144c + 2319282339840) = 

= c500 (26^51=^031 + ...), 
ri5 = c^^^*^ (146410* + 125060=^ + 40210^ + 576c + 31) = o2="^((llc)'' + . . .) 



It is instructive to reproduce the powers of c in these expressions by direct analysis of roots of the polynomials 
Gn{x) to demonstrate the strong correlation between different roots of iterated polynomials. We consider here 
only the simplest case: the periodic orbits of orders one and two in the limit c 0. In this limit all polynomials 
Gn{x; /) for cubic map / should turn into the same polynomials, but for quadratic map /, thus the powers of the 
polynomials should decrease appropriately: from Nn{d — 3) to Nn{d ~ 2). In particular, the degrees of Gi{x) 
and G2{x) should change from 3 to 2 and from 6 to 2 respectively. For cubic map Gi{x-, /) has three roots, one 
of them grows as i, another two remain finite, so that Gi{x) ^ c[x + 2) x{x— 1) + 0(c). (In this case the roots 
can be found exactly, they are: x = and x = ^(—1 ± Vl + 4c).) G2{x) has six roots, four grow as c 0, 

and two remain finite. However, the four roots can not all grow as because then G2{x) — Y^i=i{x — pf^) 
would grow as c^~^ = instead of turning into finite quadratic polynomial + a; + 1 + 0{c). Actually, only 
two of the four routs grow as while the other two - only as 

Still, this is only the beginning. Unless the coefficients in front of the singular terms are carefully adjusted, 
we will not reproduce correct c-asymptotics of di, c?2 and ri2. Actually, the first discriminant is simple: 

2-3-2 M \ n 



di=D{Gi)^c''-'^-j ^c" + 0(c), 

where the first factor is the standard a'^f' '^ from discriminant definition, and the second comes from squared 
differences between the singular root of Gi{x) and two finite ones. 
Similarly we can naively estimate: 

n.~c»-(---)V---^)'(--,»)'(.<..--V. 



4 



4 



•(^p^i^^-^j [p^'^-p^'^) ^ci^-^-^-^-^-^(l + 0(c))~c^(l + 0(c)), 

where the first factor is ^^^^ i)N2+NiN2) ^^^^ s ifj.yj) , the three next come from differences between the singular 
root of Gi and the roots of G2, while the remaining three terms come from the differences of two finite roots of 
Gi with the roots of G2 (here p^^^^ and p^^^ denote finite roots of Gi and G2, there are two and two). However, 
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here we run into a problem: the asymptotics is wrong, actually ri2 ~ c^(l + 0(c)). The reason for this is 
that the factor (^ — ^) is actually finite, ~ c", rather than ^ as we naively assumed, because the singular 
roots of Gi and G2 are strongly correlated: the root of Gi is — i + 0(1), while the two most singular roots of 
G2 have exactly the same asymptotics: — -i + 0{^/c)^. coefficients in front of ^ are the same for all these three 
roots! 

For D{G2) we have: 

D(G2) = dlru ~ 




(t - t) ' (;;3TT7W) ' (1 + ^ (t - t) ' + 

We know already, that the remaining difference is not since both most singular roots of G2 are —- with 
identical coefficients (—1). However, there is more that: they are actually —-±iy/c+ . . ., so that the difference is 
actually ~ ^Jc, and D{G2) ^ c^*'. There is no similar cancellation between the two less singular roots, because 
they are actually ~ with opposite rather than equal coefficients. Thus, since we know asymptotics of 

D{G2) and ri2, dj - 0^(1 + 0(c)). 

In this case the Mandelbrot set has no Z2 symmetry (like for the family + c), therefore we expect that 
all the cusps belong only to the zero-level components cri^^' , i.e. q{p) = 1 for p = and q{p) = for p > (like 
for the family + c). The sum rules ()7l-i(l imply in this case: 

AA(°)+^degJ 

' n,n/m ] = degJG„(£,c)], 

771 

AA^o) =deg,[d„], 
J^AAt^^^^deg.K,^/™] 

p>l m 

and the consistency condition 

degc[d„] -|-^deg^[r 

] = degJG„(e,c)], 

m 

all obviously satisfied by the data in the table, provided the number of elementary domains is counted as the 
c-degree of G„(x, c) with small c-independent x = e (exactly at critical point a; = all Gn{x, c) vanish, and at 
another critical point x = — 2/3c they have negative powers in c). 

6.4.4 — x^ + ^x^ 

The family cx^ + x^ is diffeomorphic (i.e. equivalent) to x^ + 72;^ with 7 = Xj ^fc^ which is much simpler from 
the point of view of resultant analysis. 

This time the map has two different critical points, {w/} = {0,-27/3}. Therefore one could expect that 
the number of elementary domains is counted by a sum of two terms: 

# of elementary domains = deg..y[G„(0, 7)] + deg. 

However, G„(a; = 0,7) vanishes identically because of the high degeneracy of the map / = a;"^ + 70;^ at a; = 
and does not contribute. Also two out of three roots of G\ (—3-, 7) for n = 1 coincide (and equal zero), and 
label one and the same elementary domain. Thus, actually, 

# of elementary domains = deg. 





97 



f = x^ + ■yx^ 


n 


1 2 3 4 ... 


# of orbits 
= dega, [Gn]/n 

# of el. domains 
= deg^[G„(-%^,7)]-5„,i 

dn 

total # of cusps 
= deg..y[dn] 


3 3 8 18 

2 6 24 72 

-y2 _,_ 4 j(4^4 _ ;^3^2 + 32) * * 

2 4 20 62 


'''n,n/m ^^^^ 
total # of (n,n/m) 

touching points 
= deg^[r„_„/^] for 

m=2 

m=3 

m=4 


ri2 ''24 

37^ + 16 - 2O7® + 477* - 2207^ + 1000 
2 - 6 

77"* + 687- + 169 
4 

ri4 

57'' + 4472 + 100 
4 


J^-n'^' # of sl. domains 
with cusps 
p=0, q=l 
p=l, q=0 
p=2, q=0 
p=3, q=0 


2 4 ? 20 ? 62 ? 
? ? ? ? 
? ? ? ? 
? ? ? ? 


ds = 7168a^° + 427520^* - 1230720^^ + 2330400^" + 
+25850400^2 - 9493768ai° + 33146817a* - 1078329760*^ + 
+315356704a* - 7658917760^ + 1235663104 

di = 23192823398407''2 + 993812482621447^° + 17820882315509767^^* + 
+171325147536424967^5'^ + 943919689682124807^'' + 3141248453967872007^^ + 



+9336001633252802567^" + 43599848362674749447''** + 161531609560100372487*" + 
+259103473274553630727'*-* + 190114083017294479367"^ +930486347334150389767*" + 
+1656607978055896596487^* - 1831142674577531617287^** - 2958752642287629066247^" + 
+6794824688429177434112732 + 3383523298278258150246473° + 4354843483747360581459272^ + 
+31982454515657753363720072*' + 84722721943829250141604872* - 391359614494908199135640722 - 
-484455819709207100491711972° -223935602294557922308802767^* +449858775954129150980274471** + 
+318292811209086566680166407" + 5717998165005232929505280007^2 _ 2515218215953699438592000007^° + 
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+11737578026918708510720000007* - 130446631265391083520000000007® +86930028119877222400000000007* - 
-396316767208603648000000000007^ + 180143985094819840000000000000 



As already mentioned in the previous subsection 16. 4. Hi in this case the Mandelbrot set has no Z2 symmetry 
(like for the family a;^ + c) , therefore we expect that all the cusps belong only to the zero-level components crif'' , 
i.e. q{p) — 1 for p = and q{p) = for p > (like for the family x'^ + c). The sum rules (|73(l imply in this case: 

m 

A/'i") =deg^[d„], 

=^deg^[r„,„/„] 

p>l m 

and the consistency condition 

deg^[d„] + ^deg^[r„^„/„] = deg. 

m 

all obviously satisfied by the data in the table. 



27 
3 ' 



27 
3 ' 



Comparing the data in tables in this section and in the previous s l6.4.3l we see, that the numbers of elementary 
domains differ by two - despite the two maps are equivalent (diffeomorphic) in the sense of eg . (1681) . This is the 
phenomenon which we already studied at the level of quadratic maps (which are all diffeomorphic) in s 16. 3. 51 
different families of maps (even diffeomorphic) are different sections of the same entity - universal discriminant 
variety, but what is seen in the section depends on its (section's) particular shape. Relation between cx^ + 
and + jx'^ is such that c = 1/7^, so it is not a big surprise that in the 7 plane we see each domain from the 
c plane twice. 

It is instructive to see explicitly how it works, at least for the simplest case of af\ This (these) domain(s) 
is (are) the (parts of) stability domain 6*1, defined by relations H38(l . 

J cx^ + x^ - X = 
\ |3cx2 + 2x1 < 1 

for the family cx^ + x'^ and 

( x"^ + jx^ — X — 
{ \3x^ + 2jx\ < 1 

for the family x^ + 70:^. The first of these systems is easily transformed to 

/ |3 - x| < 1, or a; = 3 - e'"^ 

while the second one - to 

13-7x1 < 1, orx= 

The critical point (single) of the first system, where dc/dcj) = 0, is at Xcr = 2, thus Ccr = —j ~ the zero of 
discriminant di{c) = 4c+ 1. The critical points of the second system, where d'y/d(f> — 0, are at Xcr — ±«, thus 
Ccr — —2 =F * ~ the zeroes of discriminant ^1(7) =7^ + 4. These points are the positions of the single cusp of 
a single elementary domain a^f^ in the first case and the two cusps of two elementary domains <jf± (one cusp 
per domain) in the second case. 
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6.4.5 Map fa-c = ax^ + (1 - a)x'^ + c 



We now use another possibility, provided by the study of 2-parametric families: we look at the situation with 
broken Zd symmetries and at the interpolation between Zd and Zd_i. In the family {ax^ + (1 — a)x^ + c} the 
Z3 symmetry of Julia set and Z2 symmetry of Mandelbrot set occurs at a = 1 , while at a = it is reduced to Z2 
for Julia sets and to nothing for Mandelbrot set. Accordingly, the number of cusps of elementary components 
of Mandelbrot set should decrease as a changes from 1 to 0. It is instructive to see, how continuous change of 
parameter a causes change of discrete characteristic, like the number of cusps (i.e. that of zeroes of irreducible 
discriminants (i„(c)). Part of the answer is given by Fig|2ni at a = 1 all the zero-level components cri"-' have the 
2-cusp shape; when a is slightly smaller than 1, only remote components ali'^ change their shape from 2-cusp to 
1-cusp type, while central domains stay in the 2-cusp shape; and the smaller a the closer comes the boundary 
between 2-cusp and 1-cusp shapes. Analytically, this means that positions of zeroes of c?„(c) depend on a in a 
special way. 

Transition reshuffling process ("perestroika" ) is a separate story, it is shown in Figl40land it repeats itself with 
all the components on Mandelbrot set, though the moment of transition is different for different components, as 
clear from Fig|2ni Fig^Hl demonstrates that what was a single central component Mi with the multipliers-tree 
internal structure for the family {x^ -|-c} at a = 1, deforms and splits into many such components as a decreases 
from 1 to 0. Most of these detached components disappear at c = 00 as a ^ 0, but some remain, including 
the central domain Mi and, say, the Ma-domain for the family {x^ + c}. This picture clearly shows how the 
trail structure between Mi and M3 (and its continuation further, towards the end-point c = —2) is formed for 
{x^ +c} from a single domain Mi for {x^ +c}. Our analytic method allows to describe and analyze this process 
by tracking the motion of zeroes of resultants and discriminants with changing a. For our 2-parametric family^^ 

rfi = (1 + af + 2{a - l)(a + 2)(2a + l)c - 27a^c^, 
= [(a + 3)(3a + 1) - 2(a - l)(a + 2)(2a + l)c + 27a'^c^] , 
d2 = a'^/^i [{ia* - 29a^ + 82a^ - 29a + 4) - 2a{a - l)(a + 2)(2a + l)c + 27a^c^] , 
ri3 = a22(7 + 40a + 4c - ia^c + iSa^c^ + ISa^c^ -I- 2160^0^ -I- 16c^+ 
+40a3+48ac2-t-16a^c2-324a'*c^+2ca-216a^c^+75a2+7a''+729a'*c'*+-131a^c2-2a'*c-8a2c-t-8a^c-hl5a2c2+324a3 



27g2 + 4p3 + 24p2 + 36p 4(20^* - 8a^ + Sa^ + 10a - 7) - 12a2(a - l)^c + Sla^c^ 
^^■^ = 27 = 81^^ 

and positions of zeroes are shown in Fig^] As a ^ 0, of two zeroes of di{c) one tends to i -f 0(a) (position of 
the cusp for the family {x^ + c}), another grows as ~ ^ ~ M ^i'^)' two zeroes of ^12(0) one tends 

to — I + 0(a) (position of the intersection o-j"^ n a^^ for the family {x^ -|- c}), another grows as ~ ^ ~ 

^ -I- 1 -f 0{a); both zeroes of 0^2(0) grow - as — i -1-2-1- 0(a) and ^ ~ if ^i^)- On their way, the 

pairs of complex (at a = 1) roots of ri2 and ^2 reach the real line Im c = and merge on it when the higher 
discriminants 

D(ri2|c) - (a - if (a + 2f{2a + if - 27a^{a + 3)(3a + 1) = 
= 4a^ + 12a^ - 84a'^ - 296a^ - 840^ + 12a + 4 = 4(0^ - 5a + l){a^ + 4a -I- 1)^ = ((2a - 5)^ - 21)((a + 2f - 3f 

^^Alternative technical approach (especially effective for a close to 1), can substitute our family by a diffeomorphic one, {x^ + 
'yx^ + a} with 7 = 01 = C\fa. For this family expressions for discriminants and resultants are somewhat simpler, for example 

di = -21c? - (47^ + 187)0 + (7^ + 4), 

ri2 = TJo? + (47^ + 187)a + (87^ + 16), 
d2 = i(21c? + (47^ -I- 187)0 -I- {47'' - 137^ + 32)), 
ri3 = 7290"' + (2167^ + 9727)q^ + (I67'' -I- I447'' + 3517^ -I- 27)a^-|- 
-|-(475 + 227^187)0 4- (77'' -1- 687^ + 169), 



Instead interpretation in the region of a near is less transparent. 
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and 



D{d2\c) - a^ia - if {a + 2f(2a + if - 27a^{'ia'^ - 2Qa^ + 82a^ - 29a + 4) == 40^(0^ - 



1)^ 



vanish. Actually, this happens at a 



0.208712152 (i.e. 7 = \/3) and 



15 w 0, 127016654 



respectively. After that the roots diverge again and move separately along the real line towards their final 
values at c = — I and c = 00. Similar is the behaviour of zeroes of ri3, ^3 and of higher resultants. The "main" 
splitting, shown in the central picture (with a — 1/7) of Fig l40l takes place at a « 0.141 - it corresponds to 
the splitting of "end-point" of the component Mi of Mandelbrot set (which is a limiting point at n — > 00 of 
the sequence of zeroes of resultants r„_2n)- The "final splitting", shown in the next picture (with a = 1/8) of 
Fig l40l takes place at a « 0.121, which is the zero of the discriminant D(wi_2|c) - this is obvious from the fact 
that the "end-point" c = — 2 of Mandelbrot set for {x^ -f c} is a zero of i«i,2(c), and two such end-points can 
be seen to touch at this cadre of Fig^UI Of course, there is nothing special in 'Wi^2 - except for that 1 and 2 
are small numbers and that wi,2(c) has real zeroes, - other stages of reshuffling (similar catastrophes, involving 
smaller components of Mandelbrot set) take place when other D{wn,s\c) = 0. 

From Fig^] it is clear, that for small a an exact - only mirror-reflected - copy of the Mandelbrot set for 
the family {x^ -f c} (which is located near c = 0) is formed in the vicinity of c = 00 (far to the left of Fig l41|l . 
More evidence to this statement is provided by analysis of relation H38|) . For example, for stability domain 5*1: 

|3aa;2 + 2(1 - a)x\ < 1, 
; = a; — (1 — a)x'^ — ax^ . 



At small a this domain decays into two identical parts: 



2\x\ < l + 0(a), 
X — x^ + 0{a), 



and, for x 



2(l-a) 
3a 



X, 



c = X — (1 — a)x^ 



2(l_a)i(l + ^i) 



^3 _ 



< 1, 



= (-2^-lir + l)-(i-S') + 0(a). 



In general, the map x — > ax -f (1 — a)x -I- c is equivalent to a; — s- ax + (1 — a)x + c with 



4(1 -a)3 2(l-a) _ 2(a- l)(a + 2)(2a + l) 



27a2 



3a 



27a2 



This is the peculiar Z2 symmetry of the Mandelbrot set for our family. 

To avoid possible confusion, note that Z2 symmetry, discussed in this paragraph, is not respected in Fig l4()l obtained with 
the help of Fractal Explorer program 1111 . The program obviously misinterprets some properties of Mandelbrot sets in the case of 
non-trivial families, like {ax^ -|- (1 — a)x^ + c}. Still it seems to adequately describe some qualitative features even for such families 
and we include Fig |4Ul for illustrative purposes, despite it is not fully correct. Actually, the right parts of the pictures seem rather 
reliable, while the left parts should be obtained by reflection. Fully reliable is Fig l41l but it is less detailed; it contains information 
only about several points of the Mandelbrot sets. Some of these points (located in the right parts of the pictures) are shown by 
arrows in both Figs l40l and 1411 to help comparing these pictures. 
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J. -1 , /1 \ 2 , small a 
f = ax^ + (1 ~ a)x + c, , 

' \a near 1 


n 


1 


2 


3 


4 


# of orbits 
= deg^ [G„]/n 


3 


3 


8 


18 


# of el. domains 
= 2deg,[G„(0;c)l 


2 

[i; 


4 
[2] 


16 

[8] 


48 
[24] 


total # of cusps 
= deg^[d„] 


2 


2 


12 


40 


total # of (n,n/m) 
touching points 

= deg^[r„,„/„] for 










m=2 




2 




4 


m=3 


- 


- 


4 


- 


m=4 








4 


J^n''^ '■ # of cl. domains 
with cusps 










p = 0, g = 1 

[q=2] 


2 

[1] 


2 

[0] 


12 
[4] 


40 

[16] 


p = 1, q = 

[q=i] 




[0] 


2 

[2] 


4 
[4] 


7 

[4] 


p = 2, g = 

[q=i] 




[0] 




[0] 




[0] 


7 

[4] 


p = 3, g = 

[q=i] 




[0] 




[0] 




[0] 




[0] 



Every map of our family has two critical points, x — and x — 0, i.e. x — ^^3^7^, which are both independent 
of c. Therefore the number of elementary domains, deg^[Gn{wc;c)] — 2deg^[Gn{0',c)]. In this table most 
data without square brackets corresponds to the case of small a, i.e. to the small deformation of the family 
{x^ + c}. However, even a minor deformation causes immediate switch of degrees of (i„(c) and r„ fc(c) to their 
values, characteristic for the cubic map family {x^ + c}. At the same time, the numbers q{p) of cusps remain 
the same as they were for the quadratic map family {x^ + c}. The table shows what this means for the numbers 
of elementary domains. Numbers in square brackets count elementary domains in the vicinity of the {x'^ + c} 
family, when some of components got merged (see Fig l4()ll and the numbers of cusps increased. Thus the table 
shows transition between cubic and quadratic families in terms of discrete numbers. Of course, at least one of 
the sum rules and consistency condition also change under the transition, since so do the numbers q{p): 

degjd„] = g(l)degJG„K,c)] + A^i^), 

and 

degc[rf„] + ^ degc[r„^„/„j] = g(0)deg^[G'„(wc, c)]. 
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while two other relations remain the same for all values of a: 

m 

^AA(^)=5]degJr„,„/„]. 

One more warning should be made: since transition occurs at different values of a for different components 
of Mandelbrot set (see Fig l39(l . the meaning of "close" (when we mention a close to 1 or to 0) depends on 
considered component (actually, on the value of n in the table). 

6.5 Quartic maps 
6.5.1 Map + c 

This section of Julia sheaf was used as an illustration in Figs|7|and|S| 
In this case 

Gi = X* — X + c, 

Ga = x^^ + x^ + 3cx^ +x^ + 2cx^ + 3c^x^ + x^ + cx^ + c^x + c^ + 1, 

Gi.i ^ G2 - I ^ (x'^ + X + c)ix^ + 2cx^ + + c^). 

Factorization of the last polynomial is accidental and is lifted by infinitesimal variation of the family a;^ + c, e.g. 
provided by additional term pa;. The same happens with accidental factorizations in the following formulas. 
The triple critical point is w = 0. 

wi = (?, 
U.12 =c«(c3 + 2)3(c'^ + 2c' + 2)^ 
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n 


1 


2 


3 


4 


# of orbits 
= degj,[G„l/n 


4 


6 


20 


60 


# of el. domains 
= deg^[G„(uic, c)] 


1 


3 


15 


60 


total # of cusps 
= degj,[d„] 


-27 + 256c3 
3 


65536c6 + 
+152064c3 +91125 

6 


* 

39 


? 
165 


'"n.n/m 3,nd 

total # of (n,n/m) 

touching points 
= degcK.n/m] for 










m=2 


- 


256c3 + 125 
3 




r24 

167772 16cS + 53411840c« + 
59113216c3 + 24137569 
9 


m=3 






ri3 

65536c8 - 2304c3 + 9261 
6 




m=4 


- 


- 


- 


ri4 

65536c« - 24064c3 + 4913 
6 


J^n'^' # of el. domains 

with (j(p) cusps 










p=0, q=3 


1 





9 


45 


p=l, q=2 





3 


6 


6 


p=2, q=2 











9 


p=3, q=2 















ds = (1099511627776c^^ +43679817400320^^ + 6678573481984c^ + 

+4811867160576c*^ + 1590250910976c3 + 1312993546389) ■ 
•(18446744073709551616c^'' + 1495915652227383951360^^ + 
+525955697232230481920ci* + 1038531305524659486720c^^ + 
+12353982755577863864320^2 + 8949781635344105472000^ + 
+419572901219568058368o'5 + 1819801492452326799360^ + 79779353642425058769) 
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The sum rules H73|l now imply: 

+ X!^^Sc[^«,n/m] = deg^[G„(u;c,c)], 

m 

+ 2 ^ clegJr„,„/„] = degjdj, 
J^AAt^^'^^degJ^^/J 

p>l rn 

and the consistency condition 

degjd„] +^degc[r 

n.n/m 

] = 3deg^[G'„(wc,c)], 

m 

all obviously satisfied by the data in above table. 
6.6 Maps fd;cix) = x'^ + c 

Julia sets for these maps have symmetry Zd, because pre-image /°(~^^(a;) of any point x is Zd-invariant set 
le^^a;*-"^-* , A: = 0, l|, and thus entire grand orbits, their limits and closures are Zd invariant. Sections 

of Mandelbrot set, associated with these families, posses Zd-i symmetry, because of invariance of the equations 
Gn{x] c) = under the transformations x — > e'^^x, c — > e'^^c. 

The critical point for all these maps is (multiple) Wf = Q. The values of polynomials G„ with n > 1 at this 
critical point, G„(0; c), depend on d^^^ only, therefore the systems of roots G„(0; c) — Q possess Zd-i symmetry 
(and the root c = of Gi(0, c) = c is a stable point of this symmetry). This Zd-i is the symmetry of entire 
Mandelbrot set. 

Gi(0,c) = c, 
G2(0,c) = l + c-^-i, 

G3(0,c) = l + c'^-i(l + c''-i)', 



G4(0,c) = 1 + c' 



,d~i 



1 + c'^-^ (l + c'^-'^Y ] -1 



d 



d 



= 1 + V (1 + c'^-Y^' 

s— 1 ^ ' 



Degrees of these polynomials satisfy recurrent relation 



r(n)-\ 

deg,[G„(0,c)]=d"-i- degJGfc(0, 



which is similar to relation for powers of G„(a;), 



T(n)-1 



fei» 



Therefore degJG„(0,c)] = iV„(d)/d. 

Note, that f„(/(0);c) = /(F„(0;c)) - /(O), and if F„(0;c) = 0, then F„(/(0);c) = 0. For our maps 
f{x) = x'' + c, /(O) = c, and F„(0; c) implies F„(c; c) = 0, actually F„(c; c) = F^(0; c). 
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The sum rules lf7S|) imply: 

■^n^ + X!^^Sc[^«,n/m] = deg JG„ , c)] , 

m 

degjd„] = (d- 2)deg,[G„K,c)] +AAio), 
5]AA(f)=^degJr„,„/„]. 

Consistency condition 

degjdn] + ^ degJr„^„/„] = (d - l)deg^[G'„(u;c 

m 

All these relations are satisfied by the data in the following table: 
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f = x<^ + c 


n 


1 


2 


3 


4 


5 


6 


Nn{d) 


d 


did - 1) 


d(d2 - 1) 


d2(d2-l) 


d(d4 - 1) 


d(d2 -l)(d3+d-l) 


# of orbits 

= dcg JGn]/n 
= N„id)/n 


d 


d(d-l) 


d(d^-l) 


d^Cd^-l) 


d(d*-i) 


d(d2-l)(d='+d-l) 




2 


3 


4 


6 


6 


# of cl. domains 
= dcg^ [G„('UJc, c)] 
= (d) Id 


1 


d-1 


d^-l 


d(d2 - 1) 


d^-l 


(d2-l)(d3+d-l) 


tot&l ^ of cusps 
= deg^[(i„] 


d — I 


(d — l){d — 2) 


(d - l)(d2 - 3) 


(d^ - l)(d^ - d-1) 


(d - l)(d'' - 5) 


(d — l)(d^ — 2d2 — 3d + 2) 


total 4t at (ri n/m") 

touching points 
= deg^[r-„_„/^] for 














m=2 




d-1 




(d-l)2 




(d-l)3 + 2(d- if 
= (d-l)2(d+l) 


m=3 


- 




2(d- 1) 




- 


2(d-l)2 


m=4 








2(d- 1) 






m=5 










4(d-l) 




m=6 


- 


- 


- 


- 


- 


2(d - 1) 


J^n^ '■ # of cl. domains 

a'*'' r 9 

with cusps 














p = 0, q = d — I 


1 





(d-l)2 


(d-l)2(d+l) 


(d- l)2(d2 +2d + 3) 


(d - l)d(d3 + d^ - 2) 


p = l, q = d-2 





d-1 


2(d - 1) 


2(d- 1) 


4(d-l) 


2(d - 1) + (d - 1)3 
= (d- l)(d2 -2d + 3) 


p = 2, q = d-2 
















2(d-l)2 + 2(d-l)2 
= 4(d-l)2 


p = 3, q = d — 2 





















Of the total of {d — l){d + + d — 1) domains ae there are: 

• (d— + rf^ — 2) zero-level components cTg*^^ - the centers of new isolated domains Msc,; 

• 2{d — 1) first-level components attached to the single a'^^ at 2{d — 1) zeroes of the resultant rig, 
located on the boundary of one between every of — 1 cusps and d — 1 zeroes of ri2 (where the ctj^^ 
components are attached to cr^°^; 

• (d— 1)^ first-level components with d—1 attached to each of the {d— 1)^ zero-level components 

at {d— 1)^ of {d— l)^{d+ 1) zeroes of rse in exactly the same way as d — 1 components a^^ are attached to the 
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Fi gure 44: A fragment (fibers shown at several points only) of the Julia sheaf for the family + c. Julia sets have symmetry 
Zs, Mandelbrot set has Z4 symmetry: a) c = -0.6, b) c = 0.4 + 0.62i, c) c = 0.269 + 0.72j, d) c = 5827 + 0.8834i, e) c = 0.6 + 0.6i. 

central ctJ''^ at zeroes of ri2] 

• 2(d— 1)^ second-level components <J^q \ with 2((i— 1) attached to each of the d—\ first-level components 
CT2^'' at zeroes of in exactly the same way as 2{d — 1) components are attached to the central af^ at 
zeroes of ria; 

• 2((i — 1)^ second-level components (J^^\ with d — 1 attached to each of the 2{d— 1) first-level components 
CTg^'' at remaining 2(d — 1)^ zeroes of rgg in exactly the same way as 2{d — 1) components cr^^ are attached to 
the central (tJ°^ at zeroes of ri2. 

For fragments of Julia sheaf for various families see FigsEI and El {d = 2), 031 ("^ = 3), 15.51 (d — 4), ^] 
(d= 5) and|331(d = 8). 

6.7 Generic maps of degree d, f{x) = Z^^Lo ^a^^ 
It is useful to introduce the following gradation: 

X — ^ A X, 
fix) ^ X-'f{x), 
on — > X~^ai, 
F„(x) ^ A-iF„(x), 
Gi(x)^A-1Gi(.t), 
Gfe(a;) ^ Gk{x) for fc > 1 

Then, since Gk(x) is a polynomial of degree Nk{d) in x, it is clear that 

Gu{x)=a7^=^ n (^-^^'^) 
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Fi gure 45; A fragment (fibers shown at several points only) of the Julia sheaf for the family + c. Julia sets have symmetry Zg , 
Mandelbrot set has Z7 symmetry: a) c = 0.7, b) c = 0.73 + 0.25i, c) c = 0.068 + 0.9775i, d) c = 0622 + 0.543i, e) c = 0.786 + 0.504i 



where the roots pi ^ \ ^ pi are sophisticated algebraic functions of the coefRcients a^. Consequently, the 
resultant 

(w„(d)-^„^i)Wfc(d) + jv„(d)(Wfc(d)-afc 1) N„{d) Nk{d) 

R{G,,,Gu)^a, n -/^i') (79) 

1=1 i'=i 

Assuming that n > k we can simplify the expression for degree of ad to Sk,iN^{d) ^ 

Example. As a simple application of these formulas, take d — 2. Then, as 02 0, all the roots grow as 
a^^, p;'^^ ~ 0^2^- The only exception is one of the two fixed points which remains finite, but its differences 
with other roots are still growing as a^^. Then from H79() we conclude that in this case 

R{Gn Gk) ~ Q,(2JV„(2)iVfc(2)-5,,iJV„(2))^-W„(2)7Vfe(2) ^ ^^„(2)(Arfc(2)-5,,i) ^.g^^ 

We used this result in s 16. 3. 31 above. 

When Qfd for (i > 2, then Nn{d — 1) out of the Nn{d) roots of G„(x) (points on the orbits of order n) 
remain finite, while the remaining Nn{d) — Nn{d — 1) grow as negative powers of ad- However, different roots 
grow differently, moreover, some asymptotics can coincide and then differences of roots can grow slower than 
the roots themselves: all this makes analysis rather sophisticated (see subsection 16. 4. 31 for a simplest example). 
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7 Conclusion 



In this paper we suggest to begin investigation of Julia and Mandelbrot sets as algebraic varieties in the phase 
space X and the moduli space of maps M respectively. We suggest to identify the pure topological and/or 
pure algebraic entity - the universal discriminant variety^^ "D* = TZ* of real codimension one in - as the 
boundary dM of Universal Mandelbrot set M. The discriminant variety TZ* naturally decomposes into strata 
of real codimension one, 

oo 
k=l 

which, in turn, are made out of comp/ea;-codimcnsion-one algebraic varieties 7?,^ , 



n = 



U 

m— 1 



The intersection 

is non-empty only when k is divisor of n or vice versa. This internal structure of the universal discriminant 
variety is exhaustively represented by the basic graph: a directed graph, obtained by identification of all vertices 
of the multiples tree (see s l2.3|l with identical numbers at the vertices (this is nothing but the graph, discussed 
in S I4.10I in relation to representation theory of Z-action on periodic orbits). Irreducible resultants 7?.^ are 
associated with the links of this graph and TZ*f,'s - with its vertices. Each 7?,* separates the entire moduli space 
A4 into disjoint strata, with one or another set of periodic orbits of order n being stable within each stratum. 
An intuitive image can be provided by a chain of sausages: every sausage is stability domain for particular set 
of orbits of given order, every chain is bounded by particular TZ'^, two chains, labeled by n and k touch along a 
spiraUing Hue, represented by 7?.* j., see FigQB] 

Ordinary Mandelbrot sets, familiar from the literature, are obtained as sections of this universal discriminant 
variety by surfaces of low dimension, usually of complex dimension one. In such section a single stratum and 
a single irreducible resultant 7?.* ,„„ can look as a set of circles and points respectively (labeled by additional 
non-universal parameters a in s.3). Irreducibility of resultants guarantee that all these points with given n and 
TO belong to one and the same entity. 

While such consideration seems to be potentially exhaustive for Mandelbrot sets (though their trail structure, 
shapes, fractal dimensions and Feigenbaum parameters still need to be better described in algebraic terms), it 
is not so clear about Julia sets: we did not identify them as (infinite) unions of well defined domains, and we 
did not establish full control over their geometry. Instead, see s l5.5l we described Julia sets as strapped discs in 
X, with identified points at the boundary. Identification of every group of points separates some sectors from 
the disc, and since at every bifurcation infinitely many groups of points are identified, the emerging structure 
of the Julia set looks fractal, see Fig03 Identified arc points of various grand orbits, and orbits involved in 
this procedure depend on the place of the map in discriminant (Mandelbrot) variety: Julia sets form a kind of a 
sheaf over universal discriminant variety. We made just a few steps towards description of the structure of this 
sheaf, its monodromy and singularity properties. We came close to description of the underlying combinatorial 
structure, which is a sheaf of rooted trees (skeletons of Julia sets) over the basic graph (the skeleton of the 
discriminant /Mandelbrot variety) . 

Even more subtle is the issue of effective action and adequate description of cyclic and chaotic RG fiows. 
Exact relation between r-functions, phase transitions and Mandelbrot sets remains to be found. Generalizations 
to arbitrary fields X and possible relations to multi-dimensional flows are not explicitly worked out. Last, but 
not the least, the obvious relations to landscape theory, which studies the distributions of algebro-geometric 
quantities on moduli spaces and their interplay with rcnormalization-group flows, are left beyond the scope of 
the present paper. We are just at the beginning of an interesting story. 

^^Since all resultants of iterated maps appear in decompositions of their discriminants, see eQ. I22l . it is a matter of taste, which 
name, resultant or discriminant variety to use. However, to avoid confusion one should remember that elementary irreducible 
components are resultants 7?.* ^ , and irreducible discriminants are rare (they can be identified as regularized diagonal 7?.* ^ . 
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Ill 



Figure 47: Getting Julia set by strapping a ball. Actually there are double-infinity of strappings: going in the direction of 
particular phase transition we strap the points at the boundary, associated with different orbits; and even if a direction is fixed, 
infinitely many points, belonging to entire grand orbit, should be strapped. 
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